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Abstract 

We consider 2- and 3- dimensional Schrodinger or generalized 
Schrodinger-Pauli operators with the non-degenerating magnetic field 
in the open domain under certain non-degeneracy assumptions we 
derive pointwise spectral asymptotics. 

We also consider asymptotics of some related expressions (see 
below). For all asymptotics loops rather than periodic trajectories 
play important role. 



Introduction 

In this paper we consider 2D- and 3D-magnetic Schrodinger operator (13.1.1) 
satisfying assumptions (13.1.2)— (13.1.5) and consider pointwise asymptotics 
of e(x, x, 0) and also asymptotics of expression (6.3.4): 



(0.1) I : 

as well as of 
(0.2) J 



J J e(x, y, 0)e(y, x, 0)u(x, y)^(x)V»i(y) dx dy 
' IS e ( x ' x ' ) e (- )/ '- )/ ' ) a; ( x '- )/ )^ 2 ( x )^ 1 (- )/ ) a,xa, ^ 



with a function cj(x, y) satisfying assumption (6.3.5): 
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(0.3) cu(x,y) := Q(x,y;x — y) where function Q is smooth in 6(0, 1) x 
6(0, 1) x 6(IR 2 \ 0) and homogeneous of degree — k (0 < k < 2) with respect 
to its third argument, 

and with smooth cut-off functions ipi,ip2- 

These two expressions play a role in the applications to the multiparticle 
quantum theory of Part VIII. Actually, instead of asymptotics of (0.2) we 
consider related estimates of 

(0.4) K : = 

(e(x, x, r) - h- d N x {T)) (e(y, y, r) - h~ d N x {t))u{x, y) dxdy 

where M x {r) denotes some approximation to e(x, x, r). 

We assume that V/F satisfies non- degeneracy conditions (13.3.45) and 
(13.3.54) i.e. 

(0.5) V < -e in 6(0,1), 

(0.6) \F\ > e in 6(0, 1), 

(0.7) |V^| x 1 in 6(0,1); 

in asymptotics of expressions (0.1)-(0.2) we will be able to drop (0.5) and 
replace (0.7) by 

V V 

(0.8) |V — | < e =^ |detHess— | > e 

or 

, V V 

(0.8) + |V — | < e =^> detHess— > e . 

As we consider Schrodinger-Pauli operator (13.5.3) and [ih > 1 conditions 
(0.7)-(0.8) + will be modified to 

(0.7)' |-^ + (2m + l-?)| + |V^| xl in 6(0, 1) Vm G Z+; 

V V V 

(0.8) |- + (2m + l-0l + |V-| <e | det Hess -| > e , 

(0.8) + I— + (2m + 1-^)1 + |V-| < e => detHess- >e 
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(a) prolate cycloid (b) perturbed prolate cycloid 

Figure 1: Classical trajectories 



respectively. 

Non-degeneracy assumptions (0.7)-(0.8) + eliminate periodicity at least as 
|jK<1 but there are plenty of loops as on figures la and lb and we know that 
loops play the same important role in the pointwise asymptotics (of e(x, x, 0)) 
as periodic trajectories play in the local asymptotics (of J e(x,x, 0)^(x) dx). 
There also will be many paths from x to y and back and these paths play 
important role in the asymptotics of expression (0.1). 

Fortunately, even if there are plenty of looping points they are rather 
exceptional on the classical trajectory. Surely, one needs to remember that if 
we study e(x, x, 0) point x is fixed but direction £ varies and thus trajectory 1 ) 
moves so it passes through x and remains on the energy level (see figures 
2a-2b below). However for the geometrical simplicity we will sometimes talk 
about x moving along trajectory. As we study expression (0.1) x,y both 
move. 

In fact, situation is more complicated than this due to the uncertainty 
principle: as fi increases not cyclotrons becomes smaller and tighter but the 
lines retain their width and become wider in comparison with the cyclotron, 
and symbolically we can consider figures 3a-3c below. 

We also generalize section 6.3 and we consider expression (0.1). It follows 
from section 6.3 that I = X w + 0(h 1 - d - K ) as fi = 0(1) with X w defined 
by the same formula but with e(x, y, r) replaced by 

(0.9) e^{x,y,r) := {2nh)- d [ e ih x y - </£. 

Jg{y,i)<v{y)+r 

^ Even if one calls them trajectories they are only x-projections of actual trajectories 
which live in (x, £)-space. 
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Then the standard rescaling technique implies the same asymptotics but 
with the remainder estimate 0(fih 1 ~ d ~ K ) provided 1 < /i = o(/? _1 ). 

Let d = 2. Then in the general case it is the best remainder estimate 
possible while 0(fih~ 1 ) is the best possible remainder estimate for 

(0.10) J e(x,x,T)ip{x)dx; 

recall that this estimate cannot be improved in the case of constant magnetic 
field F and g Jk = const, V = const. 

On the other hand, in the generic case the remainder estimate for (0.10) 
is o(fi~ 1 h~ 1 ) and the principal part is J e MW (x, r)ip(x) dx if F does not 
vanish and \i < hr 1 . 

Meanwhile for d = 3 the remainder estimate for I' is 0(/?~ 2 + / u/j~ 1 ~' 5 ) 
and under very mild non-degeneracy assumptions it is 0(/?~ 2 + /i/?^ 1 ). 

Therefore our purpose is to get the sharper remainder estimate for (0.1) 
under the same non-degeneracy assumptions. This is a very daunting task 
since for (0.10) periodic trajectories were the main source of trouble and 
they were broken in the generic case; for (0.1) loops are also the source of 
trouble, and in the generic case former periodic trajectories generate a lot 
of loops. 

It is rather easy to derive asymptotics with the sharp remainder estimate 
0(fi~ 1 h~ 1 ~ K ) as d = 2 and 0(h~ 2 ) as d = 3 but with the principal part 
given by very implicit Tauberian formula when l T is defined by (0.1) but 
with e(x,y, r) replaced by its Tauberian approximation 

(0.11) e T (x,y,T) :=h- 1 I F^ h - w (x T {t)u(x, y, t)) dr; 

J — oo 

such remainder estimate is a rather easy corollary of the results of Chapter 13. 
However deriving of an asymptotics with the sharp remainder estimate and 
rather explicit principal part is much more difficult. 

Plan of the paper 

In sections 1 and 5 we consider respectively 2D- and 3D-pilot-models with 
the constant g' k and F and linear V and derive representations of solutions 
and pointwise asymptotics as magnetic field is weak enough. We also derive 
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asymptotics of expression (0.10) with ip{x) replaced by ip 1 (x) which is ip(x) 
scaled with a spatial scaling parameter 7; these asymptotics are instrumental 
in the asymptotics of (0.2) and estimates of (0.4). 

In sections 2 and 6 we prove the same asymptotics respectively for the 
general 2D- and 3D-operators; as magnetic field is strong we use explicit 
expressions for the pilot-model operators as approximations. 

In sections 3 and 7 we prove important preliminary results and derive 
sharp asymptotics for expression (0.1) but with the Tauberian principal 
part. 

In sections 4 and 8 we pass from the Tauberian approximations to Weyl, 
magnetic Weyl or pilot-model approximations. 

Finally, in section 9 we consider estimates of (0.4) with different ap- 
proximations, leaving to the reader very similar asymptotics of expression 
(0.2). 

The results of this paper are not always sharp or very explicit, but could 
be made either sharp or completely explicit. 



1 Pointwise asymptotics: 2D-pilot-model 
1.1 Pilot-model in M 2 : propagator 

Consider the pilot-model operator 

(1.1) A = A := h 2 D\ + (hD 2 - /ix^ 2 + 2ax v 

We are interested in the Schwartz kernel U(x,y, t) of the propagator 
e lh tA . Making h- Fourier transform with respect to x 2 >— > £2 and rescaling 
xi 1 — y fixi, t 1 — y fit we arrive to 

(1.2) U(x,y,t) = (2tt/7)-V f u{x l ,y l ;r ] ,t)e ih - 1 ^-^dr ] 

with u(x 1 ,y 1 ;?], t) the Schwartz kernel of e' h lfa with ID-Harmonic oscillator 

(1.3) a = h 2 Dl + (xi - t]) 2 + 2aiJL~ 1 x 1 = 

h 2 D 2 + {x 1 -r] + aii' 1 ) 2 +n~ 1 a (2r] - aij," 1 ) . 
^ v ' ^ v ' 
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Recall that for the Harmonic oscillator b = D 2 + x 2 the Schwartz kernel of 

e itb 2) ig 

(1.4) (27r)-5 e f-( t )| sin (2t)|-| exp(-^ (cot(2t)(x 2 + y 2 ) - 2xy csc(2t))) 

where a(t) = 1, 2, 3, 4 at (0, §),(§, vr), (vr, ^ ), , 2tt) respectively and it 
is 27r-periodic and therefore for Harmonic oscillator b' = h 2 D 2 + x 2 the 
Schwartz kernel of e lth b is 

(1.5) (27rft)-*e^ CT(t) |sin(2t)|-?x 

exp^-^ _1 (cot(2t)(x 2 +y 2 ) - 2xy csc(2t)) j ; 

one can prove it easily by rescaling x i— > xfrk. 

For operator (1.3) the Schwartz kernel of e lth a is obtained by plug- 
ging x := xi — r] + cuff 1 , y ■= yi — + a^ 1 and multiplication by 
exp^Ti -1 // -1 *^ 7 ?))- Therefore, multiplying (1.5) by 

(1.6) (27rn) _ V 2ex p('^ 1 (A i_1 ttt(2^ - a/i -1 ) + (x 2 - y 2 )v)) 
and integrating by r] we arrive to 3 - 1 

(1.7) U{x,y,t) = {Inh)-^ 2 J e^^lsmilt)]- 1 ^ 1 ^ 1 ^^ dr] 
with 

(1.8) £ := ~ cot(2t)((xi - 77 + a^ 1 ) 2 + ( yi - 77 + a^ 1 ) 2 ) + 

csc(2t)(xi - r/ + ayu _1 )(yi - rj + a/i" 1 ) + (x 2 - y 2 )v+ 

tfi~ 1 a(2r] — ttyU -1 ). 

We can rewrite (1.7)-(1.8) after integration by rj as 

(1.9) U(x, y, t) = /(47r/i)-y csc(t) e' 7 ^*^ 

2 ) This is a metaplectic operator, see 1.2.3 or H. ter Morsche and P. J. Oonincx [MO]. 

3 ) After rescaling x i-» [ix, y fiy, t n- /it, fi n- 1, /? n- h — fih; however we treat 
U(x, y, t) and e(x, y, r) as functions, not as densities with respect to y and it leads to the 
factor fi 2 . 
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with 

(1.10) 0:=-Jcot(t)(xi-yi) 2 

+ ^{x 1 +y 1 +2afi~ 1 ){x2-y 2 +2tafi- 1 )-^cot(t)(x2-y 2 +2tafi' 1 ) 2 -ta 2 fi- 2 ; 
here the critical point with respect to rj is 

(1.11) r] = ^(xi + yi + 2a//- 1 ) - ^ cot(t)(x 2 - y 2 + 2to /U - 1 ). 
So, we have proven 

Proposition 1.1.1. For the pilot-model operator (1.1) 3 - ) mM 2 the Schwartz 
kernel U(x,y, t) of e' h tA is given by (1.9)— (1.11) and the Schwartz kernel 
e(x, y, t) of the spectral projector is given by 

(1.12) d T e{x,y,r) = {2nh)- 1 F t ^ h -i T U{x, y, .) = 

pirh)- 1 J e- mrltT 'U{x,y,t) dt 

and 

(1.13) e(x,y,T) = (27rh)- 1 J T (J e' 1 ^' U(x, y, t) dt) dr = 

(27r)- i y"(-/t)- 1 e-' 7i " ltT (y(x,y,t) t/t 

mt/i £/ie iosi integral taken in the sense of the essential value at 0. 

In the case of the strong magnetic field we will need an alternative 
representation. Starting from formula for ID-harmonic oscillator 

(1.14) e(x 1 ,y 1 ,r)= ^ v m { Xl )v m { yi )d(r - {2m + 1)) 

meZ+ 

with Hermite functions -u m we after the same rescaling and transition to the 
pilot-model as before arrive to 
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Proposition 1.1.2. For the pilot-model operator (1.1) 3 - ) in~R 2 the Schwartz 
kernel of the spectral projector is defined by 

(1.15) e(x,y,r) = 

(27r)"V 2 ^ -1 XI / ^( 7 ^ + ^^( x i-yi)) t; '"(^-^k x i-yi)) x 

e(r - a/i-^Xi+yj) - la^hkr] - a 2 ^ 2 - (2m + l)h}e m ~ l ^- y2)7) cfy 
and i/ius 

(1.16) <9 T e(x,y,r) = 

(47ra) _ V 3 ^ - ^ Yl J ' v ^( r i + h ~H x i - yi)) v m(v - fr~H x i - y0) x 



meZ+ ' 



g/ft 2(x 2 -y 2 )r? 



with 

(1.17) ?? m := (2a) _ VH( r - a/x _1 (xi + y x ) - a 2 ^ 2 - (2m + l)ft). 

1.2 Tauberian estimate 

Consider now pilot-model (1.1) in 6(0, l)cXc M 2 . As dynamics (classical 
or microlocal) starts in 6(0, e) it is confined to 6(0, 1) for \ t\ < T* = e/i and 
thus in this time interval we can use formulae of the previous subsection 
(modulo negligible term). 

1.2.1 Preparatory estimate 

Let us rescale 3 ) and set x = y = in (1.9)— (1.11): 

(1.18) 1/(0, 0, t) = /(47rft)-y csc(t) e''*" 1 ^ 
with 

(1.19) 4>(t) := t 2 a 2 fi~ 2 cot(t) + a 2 ^ 2 t\ 

(1.20) 7] = -ta/ji' 1 cot(t) + a/ji~ 1 . 
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We are interested in 

(1.21) F t ^ h -wXT(t)U(0,0,t) 
with 7=7*. Plugging (1.18) we arrive to 

(1.22) F t ^ T x T (t)U(0,0,t) = 

/(47rft)-y J Xr(t) csc(t)e /fi " 1( * (t) - T) dt = 

(Anhy 1 ^ 2 J xr(t)csc(t)exp^/^ 1 (-t 2 a 2 /i _2cot (f) + tp~ 2 c? - tr)) eft. 
Applying stationary phase to 

(1.23) <p(t) := -t 2 p~ 2 a 2 cot(t) + tp~ 2 a 2 - tr 
we get 

(1.24) p~ 2 a 2 (t 2 - tsin(2t)) = (r - //~ 2 a 2 ) sin 2 (t). 



Remark 1.2.1. Consider corresponding classical trajectory: 
(1.25) X! = pcos(2s) + rj — ctp~ x , £i = psin(2s) + ctp^s 

and there is a self-intersection iff 
| cos(2 Sl ) = cos(2s 2 ), 

| — ps\n(2s 1 ) + 2a/i _1 s 1 = — psin(2s 2 ) + 2ap~ 1 s 2 ; 

then S2 + Si = 7r/c with /c G Z (because we cannot fulfill the second equation 
as Si 7^ s 2 and s 2 — Si = nk); then 2si = irk — t, 2s 2 = irk + t and the second 
equation is ps\n(nk — t) = —ap _1 t or, equivalently, p = — ap _1 t csc(nk — t). 

Further, we need to satisfy pcos(jik — t) + rj — a/i -1 = (as x 1 = is a 
level of intersection) and (77 — ap^ 1 ) 2 + 2ap~ 1 (i] — ap" 1 ) + a 2 p~ 2 = r (to be 
on the energy level r) and then p 2 cos 2 t — 2ap~ 1 p cos(nk — t) = r — a 2 p~ 2 . 

Plugging p = — ap~ x t csc(nk— t) we conclude that this equation becomes 
(1.24). We can find rj to satisfy condition pcos(2si) + rj — ctp~ l = (as (0, .) 
is a point of intersection, or equivalently, rj = — ap~ 1 t cot(t) + ap" 1 . 
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We need to justify the stationary phase method. Note first that the 
spacing between successive stationary points is ~ ir. Therefore 

(1.27) t k = -t-k, t k ~ nk, s\n(tk) ~ a// _1 r _ ^7r/f 
and 

(1.28) The number of stationary points is ~ 2ir~ 1 \a\~ 1 fi,T~i . 
Recall that 

(1.29) (p'(t) = a 2 fi~ 2 (t 2 csc 2 (t) - 2t cot(t) + l) - r 
and 

(1.30) ip"(t) = -2a 2 fi- 2 csc(t)(t 2 csc 2 (t)cos(t) - 2t csc(t) + cos(t) ) ; 
therefore 

(1.31) Let k 7^ 0; then (f"(t k ) x — 2r cot(t k ) as | cos(t^)| > C e and moreover 
(p"(tk) ~ — 2rcot(t k ) as | cos(t/c)| 3> £. 

Therefore as |sin(t k )| > e, |cos(t k )| > e stationary points are non- 
degenerate and we can use the stationary phase method. 

In the near-pole zone {|sin(t)| < e} we need to remember about sin- 
gularity as sin(t) = 0; so we need to introduce a scale i x |sin(t)| and 
the stationary phase method is expected to work only as |y?"|i? 2 > CqH or 
equivalently | sin( t fc ) j > C h or, finally 

(1.32) \k\ >k := C £ _1 £ 
with 

(1.33) e := a^ l {r - a 2 ^ 2 )^. 
Summarizing, we arrive to 

(1.34) Condition (1.32) is fulfilled for all k ^ as e > C h; otherwise 
stationary phase method fails for \k\ < k. In particular, k = C fi 2 h as 
|a| x 1, r x 1 and condition (1.32) is fulfilled for all k ^ as /i < e h~2. 

Meanwhile, there is no singularity in the near-equator zone {| cos(t)| < e} 
but there is a degeneration as cos(t) = and therefore we need to introduce 
a a scale t x | cos(t)| and the stationary phase method is expected to work 
only as |<^"|£ 2 > CqH or equivalently | cos(t)| > C max(e, h^). 
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Definition 1.2.2. (i) {|sin(t)| < max(e e, C h)} is a near-pole singular 

zone; 

(%%) {| cos(t)| < C max(e, h^)} is a near-equator singular zone; 

(Hi) {| sin(t)| < max(e £, CqK), | cos(t)| > C max(e, hk)} is a regular zone. 

We apply the stationary phase method in the regular zone and justify it 
and we estimate contributions of the singular zones. 

Note that the simple integration by parts brings factors h/((p' sin(t)) and 
h<f"/(ip') 2 and contribution of zones where these expressions are less than 
h s are negligible. However at this moment we want just to estimate. 

Remark 1.2.3. As long as fih < 1 we will assume that r x 1; while we can 
always achieve a x 1 by rescaling x h-> |a|x, /i >->■ h (->■ h we will 

need many intermediate results and it is simpler to deal with the general 
case | a | < 1 from the very beginning. Without any loss of the generality 
one can assume that 

(1.35) < a < 1 ( e < fi' 1 ). 
First, consider zone 

(1.36) {t : | sin tj < (1 -e )e|t|, \t\ > e }. 

Then \(p'\ x e 2 t 2 \ sin(t)| -2 and \tp"\ < e 2 t 2 \ sin(t)|~ 3 and both factors 
\h/(p' sin(t) | and \h(f" /((p') 2 \ do not exceed h\s\n(t)\/e 2 t 2 < h/e\t\ and are 
less than 1 as h\ sin(t)| < e 2 t 2 . This always holds in zone (1.36) if H < e. 

Then multiple integration by parts brings factor C(h\ sin(t)|/e 2 t 2 )' with 
an arbitrarily large exponent / and contribution of k-th "tick" 4 ) to (1.22) 
does not exceed 

C / u 2 r 1 (n|sin(t)|/e 2 t 2 ) , |sin(t)r 1 . 

Therefore summation with respect to t along zone (1.36) intersected with 
k-th tick returns Cfih~ 1 (h/e\k\) 1 as \k\ > k with k defined by (1.32); then 
summation by k results in Cfih^^h/e) 1 as K < e (1 < \k\) and Ciibr x k x 
Cil 2 e- x &sh>e(k< \k\). 

In the latter case as 1 < \k\ < k we integrate by parts only as 
|sin(t)| < hr l e 2 t 2 and contribution of k-th. tick to (1.22) does not exceed 

4 ) I.e. interval {t : \t — t k \ < ir — e }- 
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Cfih^ 1 \og(k/\k\) where logarithmic factor is an integral of |csc(t)| along 
k-th tick intersected with {t : h~ l e 2 t 2 < |sin(t)| < e\t\}. Then summation 
with respect to 1 < \k\ < k returns 0(/i/? -1 ^) again. 
Therefore we arrive to 

(1.37) Contribution of zone (1.36) to expression (1.22) does not exceed 
C^,h~ 1 {h/e) 1 as h< e and Qi 2 ^ 1 as h > e. 

In the same time in the zone 

(1.38) {t : |sin(t)| > (l + e )e|t|, \t\ > e } 

\tp'\ x 1 and \tp"\ < |sin(f)| _1 and both factors do not exceed fo/|sin(t)|. 
Again summation with respect to t along k-th tick intersected with zone 

(1.38) returns C/i/?~ 1 (/i/e|/c|) / as \k\ > k and summation by k results in 
Cfih-^h/e) 1 as h < e (1 < \k\) and /i 2 ^ 1 as h > e (k < \k\). 

In the latter case as \k\ < k we do not integrate by parts and contribution 
of k-th tick to (1.22) does not exceed Cjih^ 1 ; summation with respect to 
1 < |^| < k returns 0(/i 2 e _1 ). Therefore 

(1.39) Contribution of zone (1.38) to expression (1.22) does not exceed 
C^h^ih/e) 1 as h< e and Qi 2 ^ 1 as h > e. 

Note that as H < ee these considered zones (1.36) and (1.38) cover 
completely the singular near-pole zone. On the other hand, as h > ee we 
did not integrate by parts in the second zone as \k\ < k and therefore it can 
be extended to {t : | si n ( t) | > e £|t|} and we arrive to 

(1.40) Contribution of singular near-pole zone to expression (1.22) does not 
exceed Cfih^fi/e) 1 as h < e and C[i 2 e~ l as h > e. 

Claims (1.37), (1.39), (1.39) imply that as |cos(t/f)| > e we can apply 
the stationary phase method. 

Then contribution of k-th stationary point to (1.22) does not exceed 
Cfih~ 1 (h/e\k\)2 and summation with respect to k : max(l, k) < \k\ < C e~ l 
returns C/i/7 _1 (/i/e|/c|)2 x \k\ calculated as \k\ = e" y i.e. C//2£ _1 /7 _ 2 which 
is larger than what we got in (1.37)-(1.40) thus resulting in 

(1.41) Contribution of zone {t : | cos(t)| > e , \t\ > e } to expression (1.22) 
does not exceed C^e^h^i as 1 < ji < h -1 . 
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Consider now the part of the near-equator zone {| cos(t)| < eo, \t\ > eo} 
not covered by (1.37), (1.39); there (1 - e)t* < \ t\ < (1 + e)t* where 

(1.42) t* = (k* + |)7r is the closest to e' 1 number of the form (k + |)7r. 
Plugging t = t* — s with |s| < et*, | sin(s) | < e we find that 

(1.43) r-V(t) = -e 2 s{2t* - s) + sin 2 (s) + 0(^ 2 ) + 0( £ | sin(s)|). 
Consider first zone 

{s: C <|s|<e£ _1 , | sin(s)| 2 < 2(1 — ei )^|s|} 

where \ip'\ > e\s\. Integrating by parts we acquire factor (h/£\f'\) x 
(h/(e\s\)2) because £ = (e\s\)2 is the scale; thus we integrate multiple 
times by parts as \s\ > e -1 /i 2 / 3 ; so contribution of k-th tick to (1.22) does 
not exceed 

Cfxh-^hW/elk - k*\) 3 2'\k - k*\i. 

2 

Then, as e > hi summation with respect to k : \k — k*\ > C returns 
C\iE2 /? _1 (/i3 /e) 1 . On the other hand, as e < h$ summation with respect to 
k : \k — k*\ > £ _1 ^3 returns Cfie^h^ 1 x hfie^ 1 = CjjLie~2h~i. 
Further, the contribution of the zone 

{s : C < |s| < ee" 1 , | sin(s)| 2 > 2(1 + ei)£|s|} 

is estimated in the same way. Finally, as s < both terms in the right-hand 
expression of (1.43) have the same sign and we can extend zones in question 
to cover {s : sin 2 (s) x e\s\} as well and we arrive to 

(1.44) Contribution to (1.22) of the zone {s : |sin(s)| < e, |s| < ee' 1 } 
without subzone 

(1.45) {s : s > C max(l, £ _1 /p), (2 - e)es < sin 2 (s) < (2 + e)es} 
and the near-equator singular zone 

(1.46) {s : |s| < C max(l, e' 1 ^), \ sin(s)| < C max(e, h^)} 
does not exceed C jie^ h' 1 {hk / e) 1 as e > fih and CjAe'^h'^ as e < frh . 
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Obviously, 

(1.47) Contribution to expression (1.22) of zone (1.46) does not exceed 
Cfih^ 1 x max(l, e^ 1 ^) x max(£,h3) = Cfxh^ 1 x max(£:, e^h, h*). 

Finally, in zone (1.45) we can use the stationary phase method. Recall 
that <//' x sin(s) x (e\k — k*\)^', therefore the contribution of k-th tick does 
not exceed C/z/7 _1 (/i/| sin(s) | ) 2 x Cfih~ 1 (h/(e\k — A - * | ) 2 ) 2 and summation 
with respect to k : \k — k*\ < e^ 1 returns C^e^hr^ . 

Combining all these estimates we arrive to estimate (1.49) below; adding 

(1.48) \F t ^ h -i T (xt(t)U(0,0,t)) \ < C^h- 1 

which is the standard semiclassical estimate 0(/i _1 ) with an extra factor fi 2 
due to rescaling we arrive to estimate (1.50): 

Proposition 1.2.4. For the pilot-model operator (1.1) with r x 1, e < /x -1 , 
H < h' 1 m R 2 

(1.49) \F t ^r{{x T (t) -Xt)U(0,0,t))\ < C^e-'h- 1 * 
as T = e < T < 00 and 

(1.50) |F t _*-i T (xr(tM0,0,t))| < Cnh^+Cuh-H-l 
1.2.2 Tauberian estimate 

Remark 1.2.5. (i) Recall that for the pilot-model operator (1.1) in domain 
X, 8(0, £) C X C K 2 with i > Q/i -1 , fundamental solution U coincides 
modulo negligible with one in IR 2 as T < T* = to[is~ l £. 

(ii) We have T* = e fie~ 1 £ rather than T* = €qE~ x I as in Chapter 13 only 
because we rescaled t y-¥ fit; without this rescaling we would have F t ^ h -i T 
rather than F t ^ h -i T and the right-hand expressions (1.49), (1.50) should be 
multiplied by ji^ 1 . 

Due to above we arrive to the statement (i) of proposition 1.2.6 below; 
then the standard Tauberian arguments yield that the Tauberian error does 
not exceed the right-hand expression of (1.50) divided by 7"*: 
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Proposition 1.2.6. For a self-adjoint operator in J?f 2 (X) in the domain 
X, B(0,£) C X C M. 2 , I > Cofi^ 1 , coinciding in B(0,£) with the pilot-model 
operator (1.1) with r x 1, e < \i , < h^ 1 , 

(i) Estimates (1.49), (1.50) hold as T < T* = eofie^i; 

(ii) Tauberian estimate 

(1.51) |e(0, 0, r) - e T (0, 0, t)\ < C^/T 1 + ^ e' 1 \F^) yT 1 et 1 = 

C{eh- 1 + fih-^r 1 

holds with the Tauberian expression 

(1.52) e T (x,y,r) := fir 1 [ (F t ^ h - W x T {t)U(x,y,t)) dr 

J —oo 

with any T G (C £ _1 , tojie^t). 
1.2.3 Micro-averaging 

Note that the second terms in the right-hand expressions of (1.50), (1.51) are 
due to the loops and according to Chapter 13 these expressions would not 
appear if instead of l~ x e(., ., r) we looked at T(e(., ., T)ip) with the uniformly 
smooth function ip. Now we want to consider l~(e(., .,r)-^ 7 ) where is 
7-admissible function and 7 < I. 

So, let us consider l~(e(., .,r)^ 7 ). First of all note that all the above 
estimates will acquire factor 7 2 due to the integration, but there is more: 
integrating by parts with respect to xi and remembering that due to rescaling 
x i—)- fix we also rescale 7 1— > /ry we conclude that the contribution of k-th 
tick acquires factor {h/e\k\jj,^y = (h/e r y\k\) 1 as 

(1.53) \k\ > £(7) := h/e-f. 

This factor "forces" \k\ = max(/c(7), 1) in the estimate of the Fourier 
transform 

(1.54) F t ^r(U^) 
and there are four cases: 
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(a) 7 > h/e; then £(7) = 1 and expression (1.54) does not exceed 
C^/r 1 m in ((ft/ £ ) 5, l) x (/?M)V; 

(b) h/e > 7 > n~ x . Then 1 < £(7) < k = h/e and (1.54) does not exceed 

(c) min(/)/e, / W 1 ) > 7 > /). Then /c < £(7) < e' 1 and (1.54) does not 
exceed Ciih~ 1 (h/e)2 \k{ci)\ 5 7 2 = C^e~ x ^\ 

(d) 7 < /). Then ^(7) > e" 1 and micro-averaging brings no improvement 
and estimates (1.55)-(1.58) below are due to estimates (1.49)-(1.51) 
without micro-averaging. 

Therefore we arrive to 

Proposition 1.2.7. For the pilot-model operator (1.1) with r x 1, e < 
/i < h~ y in IR 2 estimates 

(1.55) f 2 |^M T ((xr(t)-Xr)r(^))| < Qie" 1 /?^) 
and 

(1.56) 7~ 2 |/ r t ^- 1 .(xr(t)r(^ 7 ))| < Qi/T 1 + C^- 1 /? T ( 7 ) 
ao/d as T = e < T < 00 where 

(1.57) /? T ( 7 ) = 
' fi2h~2 as 7 < /?, 

/x^7 - ^ as min(/?/£, yu" 1 ) > 7 > /?, 

7 _1 as > 7 > /i -1 , 

k e/T 1 min((/i/e)5, l) x (h/e^y) 1 as 7 > te -1 . 

Then the standard Tauberian arguments imply 
Corollary 1.2.8. In frames of proposition 1.2.6 

(1.58) 7" 2 |r((e(., ., r) - e T (., ., r))^ 7 ) | < Cfc/T 1 + /? T ( 7 ))r 1 . 
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1.3 Calculations 



1.3.1 Correction term 

We want however more explicit expression for (1.52) which we rewrite as 
(1.59) h- 1 f (F t ^ T Xf(t)U(x,x,t))dr+ 

J — CO 

f (/vW-ft) -1 (xHO - xt (0) *. 0) ^ 

J — CO 

and as /i < /J 5-1 the first term modulo negligible is a Weyl expression 
h- 2 Ny := (47r)- 1 /)- 2 (r - \/(x)) (provided V(x) < r - e) and the second 
expression is a correction term h~ 2 Af x , CO n and our purpose is either to 
estimate or calculate it. 

Let us apply the stationary phase method whenever in makes sense; 
otherwise we just temporarily skip the corresponding zone. Then the main 
part of the answer is 



exp^sign (p"(t k ) + /7i~V(tfc)) 

with 

(1.61) <p"{tk) ~ 2rcot(t fc ) as | cos(t*)| > Ce* 
and 

(1.62) 3:={/c^0: | sin(t k )| > ft, | cos(t fc )| > Cmax(e5, H?)} 
where dimmed restriction is just temporary. 



1.3.2 Estimating correction term 

First, let us estimate the correction term. As | sin(t fc )| > h, \ cos(t k )\ > e we 
have tk x k, sin(t^) x ek, ip" x e _1 /c _1 and the contribution of /c-th tick 
does not exceed Cfj,^e~^h~^ \k\~l. 

In comparison with the previous subsection we acquired factor \k\~ l 
due to the factor (— it)" 1 in (1.59) and it is a game-changer: summation 
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3 1 1 | | 1 

with respect to k returns C[i^e~^h~^ with the smallest possible k: as 

3 11 

h < e we sum from \k\ = 1 and get C^e^^bT^\ as h > £ we sum from 
|/c| = /c := fe -1 and get C/i/? -1 which is the trivial estimate. 

Further, contribution of /c-th tick to the error is C/x/j -1 x [h/ e\k\) r+1 i \k\~ l 
if we use r-term stationary phase approximation. Then summation with 
respect to k : \k\ > 1 returns C fih^ 1 ^/ 'e) r+ 2 as h < e. 

On the other hand, as H > e summation with respect to k : \k\ > k returns 
C\ih~ x and thus this approximation does not make any sense. Therefore 

(1.63) We assume until the end of this subsection that H < e. 

This assumption takes care of many problems: it eliminates the near- 
pole singular zone and away from the stationary points we can many times 
integrate by parts; one can see easily that contribution of zones (1.37) and 
(1.39) does not exceed Cfxh^ft/e) 1 with arbitrarily large /. 

Therefore we need an additional analysis only in the the near-equator 
zone {t : | cos(t)| < e, \t — 1*| < ee^ 1 }. It is easy to estimate its contribution 
to the correction term: in comparison with the previous subsection we gain 
factor e (from factor (— it) -1 ) and we get C/i 2 e~2h~2 x e = C[i 2 e*h~2 which 
is less than the contribution of the near-polar zone. Thus we arrive to 

Proposition 1.3.1. For a pilot-model operator (1.1) with r x 1, e < iT 1 , 
h<e 

(i) Correction term h~ 2 Af x corr does not exceed Cphe^h^ and therefore in 
frames of proposition 1.2.6 

(1.64) := |e T (x,x,0) - h~ 2 M^\ < C^h' 1 + C>MH. 

(ii) In particular, < C^i^h^ 1 as e > fi 5 h. 

1.3.3 Calculating correction term 

Situation in the near-equator zone becomes more delicate if we want to 
estimate an error in r-term stationary phase approximation. Then the 
contribution of k-th. tick does not exceed 

Cneh~\h/\ cos(t k )\ 3 ) r+1 2 x C^ l {h/{e\k - k*\)i) r+ ^ 
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which as r > 1 sums to Cfieh l (k/(e\k — /c*|)z) r+2 |/c — /c*| calculated at the 

3 r+- 

least possible value of \k — k*\; so we get Cfieh^^h/e^y 2 as we sum for 
\k — k*\ > 1, i.e. as e > f&. 

On the other hand, we get C/ieh^ 1 x h^e~ l = C^h~3 as we sum for 
\k — k*\ > foe -1 , i.e. as e < h*. This takes care of zone (1.45); analysis of 
zones covered by (1.44) is easy and its contribution to the error is lesser. 

Finally, the contribution to the error of zone (1.46) does not exceed 
Cfxeh' 1 max(e, e^h, hi) where in comparison with the analysis of the previ- 
ous subsection we acquired factor e. 

Therefore, contribution of the near-equator zone to the error does not 
exceed 

(1.65) + ( C »f\(^Y H + C ^ h ~ l as , > ni, 

[ C/J/? 3 as £ < ft 3 

and we arrive to 

Proposition 1.3.2. For a pilot-model operator (1.1) with r x 1, e < ^r 1 , 
h < e correction term h~ 2 N XiCon is delivered with r-term stationary phase 
approximation h~ 2 Af XiCorr (r) with an error not exceeding 

(1.66) C///7- 1 (ft/ £ ) r+ ^ + C i xeh- 1 {h/ei) r+ ^ + C^eh~l + C/^/T 1 

2 

as e > fts 

and 

(1.67) Cnh- 1 (h/e) r+2 ' + f ih- 1 i + Cfi^h- 1 as e<hl 



1.3.4 Micro-averaging 

In the current framework (including assumption ft < e) micro- averaging does 
not improve estimate of the correction term or the error as | cos(t)| > eo in 
it unless 7 > he^ 1 (in comparison with the remainder estimate which was 
improved as 7 > h). Then, as 7 > he^ 1 in estimates (1.66), (1.67) the first 
term gains factor (h/e r y) 1 , the last term is preserved and all other terms gain 
factor (h/'-f) 1 < e which effectively kills them and we arrive to estimates 
(1.68), (1.71) below. 
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As h < 7 < he^ 1 the first and the last terms do not improve but all 
further terms gain factor (Vt) an d we arrive to estimates (1.69) and (1.70) 
below. Finally, for 7 < /j,h micro- averaging brings no improvement. 

Proposition 1.3.3. (i) For a pilot-model operator (1.1) with r x 1, e < 
e > h, 

(1.68) R^ )7 := 7 - 2 | J (e T (x,x,0)- h- 2 (M™ +M x , con(r) )^dx\ < 

Cnh-^h/eY+^h/e-f) 1 + C^h' 1 as 7 > h/e, 

(1.69) R^ r) , 7 < Cfih-'ih/sY^ + Cfie^h-'ih/e 3 *)^ + fih-^ih/^) 1 

+ C// -1 /? -1 as £ > fti, h<-f < h/e 

and 

(1.70) R^ r) , 7 < Cfih^ih/sY^ + c /U l/)-3(/,/7) , + ck 1 /)- 1 

2 

as e < h$, h < 7 < 
fn) For a pilot-model operator (1.1) r x 1, s < fi^ 1 , e < h 

(1.71) R™ := 7 ~ 2 | j (e T {x,x,Q)-h- 2 M™)vb 1 dx\ < 

fih-^h/e-f) 1 + C/i" 1 /?- 1 as 7 > /i/e- 

1.4 Strong and superstrong magnetic field 
1.4.1 Tauberian estimate 

Consider strong /ih x 1 and superstrong fih > 1 magnetic field for a 
Schrodinger-Pauli pilot-model operator 

(1.72) A := h 2 D\ + (hD 2 - fi Xl f + 2ax x - 3/1/? 

which alternatively means that we assume that |r — 3yu/)| < C in the 
framework of the standard pilot-model (1.1). Formulae (1.15)— (1.16) imply 
that 
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Proposition 1.4.1. For a pilot-model operator (1.1) with \ih > 1 m M 2 as 

1*1 < Co, |y| < Q 

(i) e(x,y, r) = mod 0(/i~°°) as t < fih — e (the lower spectral gap); 

(ii) e(x, y, t) = e(x, y, r') and d T e(x, y, r) = mod 0( / u~°°) as 

(2n — 1) + e < t' < r < (2n + l)fih — e (other spectral gaps); 

(Hi) e(x, x, t) — e(x, x, r') x [ibr 1 as r' < (2n — l)fih — e and 
r>(2n- l)fih + e ; 

(iv) As e > the following estimates hold 

(1.73) \d T e(x,y,r)\< C^e-H-l 
and 

(1.74) \F„ h -i T U{x,y,t)\ < Cfih-H- 1 ! 

and these estimate are sharp as x = y and r is close to Landau levels 
(2n + l)fih. 

Further, as before 7"* = e fie~ 1 £. Then we immediately arrive to 

Corollary 1.4.2. For a self-adjoint operator in domain X, B(0,£) cXc 
M. 2 , £ > Co/i _5 /j2 ; coinciding in 6(0, £) with the pilot-model operator (1.72) 
with t x 1, e < 0, /j, > h^ 1 

(i) Statements (i)-(iv) of proposition 1.4-1 remain true; 

(ii) Formula (1.15) holds modulo 0(n^h~^£~ l ); 

(Hi) Moreover, formula (1.15) holds modulo 0(/i -o °) 
as (2n - 1) + e < t < (2n + l)fih - e ; 

1.4.2 Micro-averaging 

Let us consider micro-averaging. First let us estimate Fourier transform 
where as before we rescaled t h-> fit; formulae (1.15)— (1.16) imply immedi- 
ately 

Proposition 1.4.3. For a pilot-model operator (1.72) in M 2 with [ih > 1 
and 3 = (2n + 1) 

(1.75) 7 " 2 |/WirX7-(r)r(^ 7 )| < C> 2 + C^- 1 /? T ( 7 ) 
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with 

{11 _ i i 

u 2 h 2 as 7 < u 2 /)2, 

7 as 7 > // 2 /j2 . 

Corollary 1.4.4. In frames of corollary 1.4-2 
(1.77) 7' 2 |r(e(., ., r) - e T (., ., r))V 7 | < C(/ie + /? T ( 7 )> 



1-1 



1.5 Weyl and magnetic Weyl approximation 

We can try to use a host of the different approximations but restrict ourselves 
now to Weyl approximation, which makes sense only as [ih < 1. Recall that 

(1.78) |e T (x,x,r) - /r 2 A/" x w (r)| < C/? w 

(1.79) 7~ 2 | J (e T (x,x,r) - h~ 2 M^ (r))V> 7 dx\ < C/? w ( 7 ) 
with 

;i.80) R ^ = {^-h-\ a se>h, 

as e < H, 



;i.8i) /? w ( 7 ) 




as 7 < 

as 7 > h/e. 



On the other hand, let us apply instead the magnetic Weyl approximation. 
Obviously, without micro- averaging Af^ w (r) cannot produce uniform with 
respect to r remainder estimate better than fih -1 due to its own jumps. 
Meanwhile with an averaging (i.e. a micro-averaging with 7 = 1) A/" x mw (t) 
provides a better remainder estimate O^ji^h^ 1 ) (see Chapter 13; for a 
pilot-model we need only assume that e > fJT ). So, as ///? < 1 the magnetic 
Weyl approximation may provide a better remainder estimate as 7 is not 
too small enough but for small 7 Weyl approximation is better. Let us 
investigate this. 

Proposition 1.5.1. (i) For a pilot-model operator (1.1) as fih < 1 

(1.82) R MW ( 7 ) := 7 - 2 | j (e T (x, x, 0) - /T 2 J A/" MW )^ 7 dx\ < 

Cfi-'h- 1 + C/? MW ( 7 ) 
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with 



(1.83) 



iMW 



(7) 



( ^h-'j- 2 (h/e 7 )' 

IMh-^h/ej) 1 
U/7- 1 



as 7 > max(/i _1 , /?/£:), 

as /i" 1 < 7 < h/e, 

as h/e < 7 < 

as 7 < min(/z _1 , /?/e); 



fw) For a pilot-model operator (1.72) as iih > 1 and 3 = 2n + 1 estimate 
(1.84) R MW ( 7 ) < C + 0? MW ( 7 ) 



(1.85) 



(7) 



I' 2 
iih' 1 



as 7 > /i 2/ 7 2 j 
_i 1 

as 7 < /i 2/72. 



Proof. We need to investigate only case /i -1 < 7 < te -1 as in all other cases 
we can easily pass from results for R w (7). 

Note that /7~ 2 A/^ /lw (r) delivers e(y,y,r) for a pilot-model operator A : = 
A y = A y (x, D x ) defined by (1.1) or (1.72) with ax-y frozen at the point y 
(with respect to which we integrate later); so we consider our pilot- model 
operator as a perturbation of this one. 

Let us plug into Tauberian formula the Schwartz kernel of 



;i.86) 



e iftr 1 tA _ e ihr 1 tA 



+ ih' 1 f e inrlt '~ A (A - A) e inrl ^ A dt' 
Jo 



or its iteration for ±t > 
(1.86) 



e ih-hA = e 



iH-hA x+ X 
!</<m-l 



(/n- i y J e ih ~ ltlA (A - A)e m ~ lt2 ' A (A -A) ■■■(A- ~A)e ih ~^' A dhdt 2 ■ ■ ■ dtj+ 



{±ti>0,±t 2 >0 ±tj>0, 

t! + t2 + ... + tj = t} 



{ih- 1 ) m / e ih ~^' A {A - A)e iR ~ lt2 ~ A (A -A) ■■■(A -A) 



Jh^tmA 



dt x dt 2 ■■■dt n 



{±ti>0,±t 2 >0 ±t m >0, 

t 1 +t 2 +...+tj=t} 
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So let us plug (1.86) m into our expression with Xr{t) instead of Xr(t)- 
Then our standard methods imply that j-th term (J = 0, ... , m) of what we 
got does not exceed 

(1.87) C^ih- l {T\\A - KWti^y x (1 + + 

'2 '1 

with Ti = C fie~ 1; y where 7 := max(/i _1 , and 7~ 2 = /i/ £ 7i we l eave 

easy details to the reader. 

Then we can replace the norm by the "effective norm" which is a norm 
restricted to /i7-vicinity of x (in rescaled coordinates) and equal to Ccry. 

Note that T2 < T\ as 7 > 7; then summation with respect to T returns 

Cfih-^T^A-AxWh-'Y x C/i/)- 1 ^) 7 

which provides required estimate (1.83) or (1.84) for all terms with j > 2. 

Finally note that while A y (x — y, D x ) is even with respect to (x — y, D x ) 
(we can always achieve it by the gauge transformation), perturbation 
(A — A y ){x — y, D x ) is odd and therefore only terms with even j survive 
as we plug x = y. This takes care of term with 7 = 1. □ 

Remark 1.5.2. (i) For fih < 1 magnetic Weyl formula provides better ap- 
proximation as < 7 < h/e; otherwise Weyl formula provides either 
better or equally bad (/u/)" 1 ) approximation. 

(ii) For fih > 1 magnetic Weyl formula provides approximation rather than 
just mam term estimate as 7 3> \x 2/72. 

1.6 Geometric interpretation 

Remark 1.6.1. (i) As we mentioned tk are length of the classical loops at 
point (or x if we consider e(x,x, r); then t k = t k (x)). While picture of 
the fixed trajectory and a point x on it (Figure 2a) is more geometrically 
appealing, the correct picture is of the fixed point and different trajectories 
passing through it (Figure 2b); due to assumption a^O "radii" of these 
trajectories differ by 0(/x _1 ). 

Therefore, when we talk here about loop we mean a trajectory looping 
at the fixed point x rather the self-intersections of the fixed trajectory. 
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(a) Intuitive (b) Correct picture 



Figure 2: On the left more geometrically appealing picture, when x moves 
along given trajectory; on the right the correct one when we consider different 
trajectories passing through a fixed point. 

Similarly, when we talk about pole or equator we mean not x near a pole 
or an equator of the fixed trajectory but rather trajectory such that x is 
near its pole or equator. 

So, all further remarks should be interpreted correctly. 

(ii) How to interpret our results from heuristic uncertainty principle? Obvi- 
ously these results mean that ase>li then the majority of trajectories on 
Figure 2b do not return to the original point x = after k = 1 tick 5 ) while 
otherwise this happens only after k ticks with ke ^> h. As the spatial shift 
is ek (after rescaling) this means that the thickness of the trajectory is x H 
after rescaling (and therefore x h before rescaling). 

It also means that if | sin(0)| x 1 where 9 is the polar angle (between 
direction of the trajectory at and the drift direction (0, 1)) then we can 
take an interval in 9 of the magnitude 1, so 9 is a dual variable to x 2 . 
Conversely, with the exception of the interval of the length K/e\k\ around 
9k (where 9k corresponds to the classical trajectory returning to after k 
ticks) trajectories do not return at 0; we can rewrite this condition as 

(1.88) |sin(0-0 k )| > h/(\k\e). 

5 ) In the sense that the difference between two points is observable. 
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This characteristic length matches to the stationary phase estimates we 
derived rigorously. 




(a) (b) (c) 

Figure 3: Looping near pole from microlocal point of view 

(iii) Near equator situation changes drastically: spacing between self-inter- 
sections is x e\ cos(^)| _1 as | cos(0)| > eoe? and there are exceptional 0, 1, 2 
self-intersections with | cos(#)| < eo£2. Uncertainty principle may prevent 
us to know how many of them are no matter how "not large" /i is. 



(c) 



(d) 



Figure 4: Self-intersections near equator. The first winding shown by a 
solid line, the closest to it near equator by a dashed one. 



Consequent self-intersections are distinguishable as | cos(#)| > h\ accord- 
ing to our calculations. 
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Problem 1.6.2. We know that estimates (1.73), (1-74) are sharp as jih > 1. 
Are these estimates sharp as jih <C 1? 

2 Pointwise asymptotics: general 
2D-operators 

Let us recreate the approximation similar to the one of the previous section 
for general operators assuming that K < e i.e. 

(2.1) i2<eh-\ 

We are interested in \ t\ < Coji~ 1 £~ 1 for the original (not rescaled) oper- 
ator as for C ar 1 e^ 1 < \t\ < // _1 7~* := e^e^i dynamics leaves 8(0, c^jT 1 ) 
but remains in 6(0, £). Here 

(2.2) e = a' 1 or 1 as we assume that \W / F\ x a in 6(0, £) with 

(2.3) dji' 1 <a<l ( C fi- 2 < e < fi' 1 ) 
and £ > Coji" 1 with sufficiently large constants C\, c Q . 

2.1 Classical dynamics 

First we consider classical dynamics starting from point x = 0. So far we 
use not-rescaled x, t. Let us freeze g Jk and F there and replace Vj and V by 
their linear germs: 

(2.4) g* = g^(0), V(x) = V(0) + (W(0),x>, 

Vj(x) = Vj{0) + <W y (0),x> ( =► F = F(0)), 

denote corresponding Hamiltonian by a(x,£) and all the object associated 
with it will have bar. We consider original our dynamic system as a pertur- 
bation. 

Without any loss of the generality we can assume that 

(2.5) g jk = u 2 S jk ( a(x, = w 2 (p2 + pi) + V), 

(2.6) w(0) = F(0) = 1, (Vw)(0) = 0, 

(2.7) \/(0) = 0, (W/F)(0) = (-«,()) 
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as we can achieve it by an appropriate change of variables (see footnote 9)) 
and 

(2.8) Vi(0) = 0, V 2 = 

as we can achieve it by the gauge transformation. 

Proposition 2.1.1. Under conditions (2.1)-(2.8) as \t\ < c u,~ 1 e~ 1 

(i) For drift flows z t = <D t (0), z t = 4> t (0) 

(2.9) z^Zt + O^et 2 ); 

(ii) If F = 1 then for Hamiltonian flows (x tl £ t ) = W t (0, 77), (x tl £ t ) = W t (0, 77) 
with magnetic parameters \i and fx = + 0{fi~ 2 )) defined by (2.18) on 
the energy level t < c 

(2.10) (xt,rtt) = {*t,pr%) + o(fi- 2 t). 

Remark 2.1.2. Assumption (2.1) means that in 6(0, ^T l ) the constant part 
of V V dominates over its variable part and it implies that as 

(2.11) yr l et 2 < c Q iT 2 \t\ < Ke\t\ as \t\ < Cq^ 1 ^ 1 
with 

(2.12) k = co^V 1 < Q 1 ^ 

and thus the perturbation of the drift is respectively small. In particular 
drift line deviates from the straight one by no more than cqK. 

Also equator is defined as a point where x t intersect its its first winding 
the last time deviates from (/i , 0) by no more than Co/! -1 k. 

Proof of proposition 2.1.1. (i) Proof of assertion (i) is trivial and left to 
the reader. Just recall that the drift flow is defined by (13.6.5) with an 
extra factor a~ l in the right-hand expression and in our assumptions 4> t : 
(xi, x 2 ) I—?- (xx + et,x 2 ) after rescaling. 

(ii) Recall that pj = £y — /iVy(x), {p^x^} = 5^ and in virtue of assumption 
F = 1 

(2.13) { Pl ,p 2 } = flu- 2 . 
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According to subsection 13.2.1 one can correct 

(2.14) yi := Xi - {pi,p 2 }~ 1 P2 = *i ~ prVp 2 , 
y 2 ■= x 2 + {pi, P2} _1 Pi = x 2 + /i _1 w 2 Pi 

modulo 0(p~ 2 ) so that corrected expressions satisfy drift equation modulo 
0(/U~ 2 ). More precisely, in the current setup 

{a, yi } = pT V ({p 2 , u 2 }{pj - pi) - 2{pi, u 2 } Pl p 2 + {p 2 , V}) 

and for y[ = y 1 + p~ 2 /3iPip 2 + p~ 2 /3 2 (p 2 - P 2 ) with 

(3t = -^u 2 {p 2 ,u 2 }, /3 2 = ^u 2 {p 1 ,uj 2 } 

we have 

modulo terms which are 0(p~ 2 ) and also homogeneous polynomials of 
degrees 1 or 3 with respect to (pi, p 2 ) and thus these terms could be corrected 
by adding to y[ terms which are 0(p~ 3 ) and also homogeneous polynomials 
of degrees 2 or 4; then we arrive to 

{a, ya = (-1) V^ 2 {P3-/o V} mod 0(^- 3 ) /c = 1, 2 

(/c = 2 is considered in the same way as k = 1). Note that the right-hand 
expressions are calculated at point x = x t rather than y" = y"; so we rewrite 
them as 

mod 0(p~ 3 ) /c = 1,2 

and we can correct y£' by adding to y[' terms which are 0(p~ 3 ) and also 
homogeneous polynomials of degrees 2 with respect to (pi, p 2 ). 
Therefore 

(2.15) y' t " ee z t mod 0(// _3 t) 
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where z t = <t> t (0). 

Meanwhile (2.13) implies 

(2.16) j t pj = 2(-iy m _j - (u;%(p 2 + pi) - v x . = 

2(-l)Vp 3 i - 2u X] u- l W + W Xj , 7 = 1.2 

where we replaced p 2 + p\ by u~ 2 W, W = r — V (no error as r is an energy 
level). Let us plug xi(t) = zi(t) + u,- 1 u 2 (z t )p 2 , x 2 {t) = z 2 {t) - // _1 u; 2 (z t )pi 
modulo 0(/i -2 ) in the right-hand expressions; we get 

< 2i7 » sfeHfeV* 



with 



-fi-Wi^co-'W + l/l/ Xl ) X2 -(1 - ^a; 2 ^- 1 !/!/ + H/ X1 ) X1 )^ 
(l-^ 2 W 2 (a; X2 a;- 1 l^+ l/l/ X2 ) X2 ) ^(u^W + W X2 ) Xl 



J := 2// 
and 

K 



-Iuj^W + W X2 



and here we can calculate elements of J at any point of 6(0, /i -1 ) we choose 
(and we choose 0) while elements of K are calculated at z t . Then J becomes 
a constant coefficient matrix and since Vw(0) = 0, u(0) = 1 



J = 2/i 



- / u~ 2 (w XlX2 1/1/ + W X1X2 ) -(1 - ^- 2 (u; XlXl 1/1/ + W XlXl )) 
(1 - ^- 2 (u; X2X2 1/1/ + W X2X2 )) fi- 2 (u X2Xl W + W X2X1 ) 



and its eigenvalues are ±2/// with 
(2.18) /2 = 

/i ( (l-^- 2 (a; XlXl H/+ H/ X1X1 )) (l-^- 2 (u; X2X2 W+ W X2X2 )-^ 4 {lo XiX2 W+W XiX2 ) 2 

= u, + fi- 1 {A{uW)) + 0(/i" 3 ). 

Then J = Q-^C? with J = Jx ^ and (? = / + 0(/i- 2 ). 
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Due to (2.17) 

(2.19) p(t) = e tJ p(0) + f e^'V K{z t ,) dt' = 

Jo 

e tJ p(0) + J* e^'^J-^Kizt,) dt' - r l K(z t ) + e tJ J- 1 K(z ). 

The second term in the right-hand expression is 0(/j,~ 1 et) as J -1 = 0(fi^ 1 ) 
and ^ 7 K(z t /) = 0(e). Note that with 0(fi~ 2 t) error one can replace in the 
two last terms J by J calculated for a pilot-model. Finally, with an error 
0(fi~ 1 et) one can replace K(z t ) by K(z ); however K(z ) and K(z Q ) coincide 
as Vw(0) = 0. 

Replacing J by J we reduced evolution of p to those of the pilot-model al- 
beit with Jl instead of \i. Therefore fi~ 1 p t = [J>~ 1 Pt m °d 0(fi~ 2 t). Replacing 
/i by p, does not affect drift (in frames of the indicated precision). 

But then x t could be found from z t and p t and their drift also is described 
by a pilot-model in frames of the same error and then it is true for fi^ 1 ^ as 
well. □ 

Consider now the general case, i.e. F different from 1 and variable. 
Then differential equations describing (x(t),£(t)) for a(x, £) coincide with 
equations for F~ 1 a(x, £) with F i — >■ 1 and r — V i— > F _1 (r — V) but with the 
"time" 9 satisfying 

df) 

— = F{x t ) = F(z t ) + (u 2 F- 1 )(z t )(F Xl (z t )p 2 - F X2 (z t ) Pl ) mod 0(^ 2 ) 

where z t = <t> t (0). We can correct 9 by 0(u~ l ) eliminating linear terms in 
the right hand expression and therefore 

(2.20) 9 = j F(z t )dt^F(Q)t + \{^ t F(z t ))(Q)f 

mod O^ 1 + AT 2 |t|) 

with the second term in the right-hand expression 0(efT 1 t 2 ). 

From now on we use rescaling x h-> fix, t h-> Jit, 9 h-> /i9. Our analysis 
implies immediately 
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Proposition 2.1.3. Let conditions (2.1)-(2.8) be fulfilled 6 ^. Then (after 
reseating x — > fix, t ^ fit) as \ t\ < Cqe^ 1 

(2.21) \D^ t \ < Cp V/3. 

2.2 Semiclassical approximation to U(x, y, t) 

Proposition 2.2.1. Let conditions (2.1)-(2.8) be fulfilled. Then 

(i) Uniformly with respect to \t\ < ce^ 1 e lh ltA is an h-Fourier integral 
operator corresponding to Hamiltonian flow V|/ t ; 

(ii) As | sin(26>)| > e with 6 defined by (2.20) 

(2.22) U(x, y, t) = (47rn)- 1 /(sin(^)- 1 e'^ 1 ^' t ) £ b m (x, y, t)h m 

m 

with (j> defined by (2.25) -(2.29) below and satisfying (with all derivatives) 

(2.23) ( p = 4>(9) + 0( f i~ 2 e- 1 ), 

(2.24) b m = & 0m + 0(/i- V 1 ) 

with 4> defined by (1.10). 

Proof. Both assertions of proposition are standard as |r| < T = const but 
we need to extend them for larger t. Further, is an action 7 ^ 

(2.25) <f> = - L(x(t'),x(t'))dt' 

Jo 

with the Lagrangian 

(2.26) L(x, x) = Q2 x k C k - a(x, f )) ^=(2.29) = 

k 

\ Yl Sjk{x)xjk k + ^J2xjVj{x) - V{x). 

j,k j 

6 ) Before rescaling. 

7 ) Sign "— " is due to considering of propagator e' h tA rather than e~' h tA . 
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which is the Legendre transformation of the Hamiltonian 

(2.27) a(x,0 = XyOOte - Vj(x)) fa - V k (x)) + V(x) 

j.k 

(2.28) x i = 2^^(x)(6-\/,(x)), 
(2-29) &=JEE^ +V *W' 



f := c/f / c/t' and we are talking about trajectories from y as t' = to x as 
t' = t. 

To extend assertions (i), (ii) to t: \t\ < T = ce^ 1 note that 
(2.30) e 7rtn-M = _ e ,>-is- 1 e 

with ft-pseudo-differential operator 6; 6 = if A is a pilot-model operator 
with a = 0. Therefore as usual 

(2 31) e irth- 1 A _ e ikh- 1 A e i(t-wk)hr 1 A _ ^_^k e ikfj,- 1 hr 1 B ^(t-TT^h^A 

with k = |_ t/ 7rJ and the right-hand expression is h- Fourier integral operator 
as \k\ < ce^ 1 . Assertion (i) is proven. 

Note that the corresponding canonical manifold is well- projected to (x, y) 
space (1R 4 ) as | sin(2#)| > e and compare it with the canonical manifold for 
the pilot-model operator; this proves (ii) as either F = 1 or \t\ < e\e with 
sufficiently small e 1 = €i(e). 

But then one can represent e lh tA as e' h tlA e lh t2A ■ ■ ■ e lh tnA with n < 
Ci and \tj\ < e ly u, | s i n ( 2 1,- ) | > |e which implies (ii) in the general case. □ 

So far we exclude both vicinities of sin(#) = which matches x = y 
(poles) and cos(9) = which matches x and y being antipodal points of 
the trajectory (equator) but we need to approach both of them. Actually 
exclusion of the latter was no more than a precaution but poles require a 
modification: 

Proposition 2.2.2. Let conditions (2.1)-(2.8) be fulfilled. Then 
(i) Decomposition (2.22) remains valid as |cos(#)| < e; 
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(ii) Decomposition (2.22) remains valid as 

(2.32) Cmax(ft,// -1 £|t|) < |sin(0)| < e 
albeit with an error not exceeding 

(2.33) Cr 1 |sin(^)|- 1 (V|sin(^)|) / 

and with b m such that 

(2.34) \DP((f)-4>)\ < C^e|sin(0)|-^l V/3, 

(2.35) |D> m -6 m0 )| < C m ^|sin(^)|- m -^ V/3, m. 

Proof. As 

U(x, y, t) = fi' 2 J U(x, z, t') U(z, y, -t") dz 

for t = t' — t" (where factor u~ d is due to rescaling U as a function) we 
need to consider this oscillatory integral. If we consider oscillatory integral 
with the propagator for the pilot-model operator we note that the standard 
stationary phase method applies with an effective semiclassical parameter 
h as (cot(t') — cot(t")) = — sin(t) sec(t') sec(t") disjoint from i.e. also as 
| cos(t)| < e. 

On the other hand, an effective semiclassical parameter is h/\ sin(t) | near 
poles. 

In virtue of proposition 2.2.1 both these claims remain true for a general 
operator as well albeit with t replaced by 9. □ 

Remark 2.2.3. So far we need only \W / F\ < jie rather than \W / F\ x 

HE 6 ). 



2.3 Semiclassical approximation to e(x, x, r) 

Therefore in zone {| sin(^)| > Ch} all arguments of the pilot-model theory 
work (with obvious modifications) under non- degeneracy assumptions 

(2.36) V - r < -e , 

(2.37) |V(r - V)F _1 | x e 

we conclude that 
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(2.38) Contribution of zone {|sin(0)| > Ch} to F t _^-i T (l — xi(t))r x (7 does 
not exceed Cfih^ 1 + Cpae^br^ . 

Consider now zone {|sin(6 l )| < e|}. Both operators A and B (denned 
by (2.30) are £-microhyperbolic and unless factors they are coming with 
(namely (8 — irk) (with the closest irk) and et respectively are of the same 
magnitude we can use £-microhyperbolicity to prove 

(2.39) Contribution of zone {\t\ x T > e^ 1 max(fih, C\ sin(0)|)} to F t ^ h -i T V x U 
does not exceed 

Cf/h- 2 xTxeTx (h/sT) l+1 = C^h^T^/eT) 1 

(with arbitrarily large /) while its contribution to expression (1.59) does not 
exceed the same expression albeit without factor T i.e. 

Ci/h-^h/eT) 1 , 

where factor u 2 is due to rescaling, T x eT is the measure of the zone 
: |^| x 7, | sin < eT} and other factors are standard; recall that B 

comes with the factor e and therefore effective semiclassical parameter is 

h/e. We leave easy details to the reader. 

After summation with respect to T we conclude that 

(2.40) As e > Ch contributions of zone {\t\ > 1, |sin(0)| < ch} to both 
F t ^ h -i T V X U and expression (1.59) do not exceed Cfih^lh/e) 1 

and 

(2.41) As e < C h < 1 contribution of zone {|t| > C h/e, |sin(0)| < ch} 
to F t _^-i r r X U does not exceed C/z 2 ^ 1 while its contribution to expression 
(1.59) does not exceed Cfih^ 1 . 

So, as e > Ch we are done but as e < Ch we need to consider zone 
{1 < \ t\ < ch/e, | sin(0)| < ch} and the same arguments imply 

(2.42) As e < C h < 1 contribution of zone {1 < \ t\ < ch/e, | sin(0)| < ch} 
to F t ^ h -i T V X U does not exceed 

C^hT 2 x e~ x h x h = Cf^e^ 1 

where factor u 2 is due to rescaling, e~ 1 h x h is the measure of the zone 
{9 : \9\ < e^h, | si n (6*) | < h} and factor hr 2 is the standard one. 
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Combining with the results for zone {| sin(0)| > ch} we arrive to 

Proposition 2.3.1. For magnetic Schro dinger operator in domain X , B(0,£) C 
X C IR 2 under standard smoothness assumptions and non- degeneracy as- 
sumptions (0.6) and (2.36) -(2.37), (2.3) after rescaling x i— >■ fix, t i— >■ //t 
fas /x < estimates (1.49)-(1.50) of proposition 1.2.4 an d estimates 

(1.55)-(1.57) of proposition 1.2.7 hold as T < ce -1 . 
Further, estimates (1.78) -(1.81) /io/d as we//. 

Now we apply the standard Tauberian arguments. We are interested 
mainly in the case of a x 1; then automatically £ x 1 and T* x ^u 2 inside 
of domain. However in more general case t x a x fie automatically and 
7"* x i/efifi 2 x fi 2 again and we arrive to 

Proposition 2.3.2. For magnetic Schro dinger operator in domain X , B(0,£) C 
X C IR 2 under standard smoothness assumptions and non- degeneracy as- 
sumptions (0.6) and (2.36)-(2.37), (2.3) Tauberian estimates (1.51) of 
proposition 1.2.6 and estimates (1.57)— (1.58) of corollary 1.2.8 hold with 
the Tauberian expression (1.52) and t x a. 

2.4 Stationary phase calculations 

Let apply the stationary phase method to the Tauberian expression. To 
do so we need to remember that assumptions \t\ e^ 1 , tp' = yield 
ip" x | si n (6*) | 1 and all other estimates hold; therefore the stationary phase 
construction is available there. Also this construction works as long as 
\t\ < ce^ 1 and \tp"\ > e i.e. \t — t*\ > ee^ 1 where T* = e^ 1 . 

Moreover, \<p"\ > e as \ t — t*\ < ee~ l and </?' = and therefore all our 
estimates work here, in the near equator zone as well where 

Definition 2.4.1. Equator just moves to the point where \ip"\ is minimal 
for ip' = 0. Let us in the near equator zone redefine 9 in such way that 
equator is cos(#) = 0. 

Let us introduce 

Definition 2.4.2. Consider r-term representation (2.22) (i.e. with summa- 
tion over m < r) and plug it into Tauberian expression, excepting \ t\ < e and 
calculate it by the stationary phase method with r terms again. Let us call 
the result r-term semiclassical approximation and denote it by h~ 2 Af X:Corr ( r ). 
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Then its main term is delivered by modified (1.60)-(1.62) 
( 2 - 43 ) E T^=^" l (t*sin(^))-V"(tfc)r | 6b(tfc)x 



exp^sign<p"(tfc) + /7i V(tfc)) 



with 



(2.44) ^"(tfc) ~ 2rcot(0 fc ) as | cos(0*)| > Ce* 
and 

(2.45) 3:={/c^0: | cos(0 fc )| > C max(eK /p)} 

where tk is a time of /c-th return to x (along /c-th loop), 9k = 0(tk) and <£>(t/() 
is the corresponding action. 

Then repeating arguments of the previous section we arrive to 

Proposition 2.4.3. For magnetic Schro dinger operator in domain X , 6(0,1) C 
X C 1R 2 under standard smoothness assumptions and non- degeneracy as- 
sumptions (0.6) and (2.36)-(2.37) 7 (2.3) 

(i) As e > h estimate (1.64) of proposition 1.3.1 holds; 

(ii) Estimates (1.66)-(1.67) of proposition 1.3.2, estimates (1.68) —(1.71) of 
proposition 1.3.3 and estimates (1.82) -(1.85) of proposition 1.5.1 hold. 

2.5 Approximation by a pilot-model 
operator 

2.5.1 Weak magnetic field case 

We cannot do better in this framework; however using the pilot-model 
operator as an approximation we could improve these results. Without any 
loss of the generality we can assume that (2.5)-(2.7) are fulfilled. 
Assume temporarily that F = 1. Then one can prove easily that 

(2.46) tk = 9{t k ) = t k , </>{t k ) = 4>{t k ) mod/i- 2 |/c| 
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where bar denotes objects related to the pilot-model operator. Then 
— 0(tfc)) = 0(u.~ 2 h~ 1 \k\) and it is less than 1 as 

(2.47) \k\ < k := fi 3 h. 

Therefore it makes sense to apply this approach only as k > 1 i.e. 

(2.48) ^>h~^. 

Under this assumption let us compare r-th terms for our operator and the 
pilot model as \k\ < k. 

One can see easily that their difference does not exceed 

(2.49) C^iih-\hle\k\y^\k\- x x /i" 3 /)- 1 ^! 

k 

which as r > 1 does not exceed C fi~ 2 h~ 2 (h/e\k\y + ^\k\ calculated for 
the smallest possible \k\ which is 1 provided H < e. Therefore we get 
Cfi~ 2 h~ 2 (h/e) r+ 2. On the other hand in (2.49) one should replace the last 
factor by 1 as \k\ > k. 

Therefore, if we start from /-term approximation and apply estimate 

2 2 

(1.66), (1.67) as e > h*, h < e < fa respectively and pass from /-term to 
r-term approximation we get estimate 

(2-50) R^ )7 : = 

|e T (x,x, r) - /) _2 A/" XiCorr(r ) - e T (x,x, r) + h' 2 M XtCon{r) \ < 

Cfi-H- 1 + Cfi- 2 h- 2 {h/ey +1 i + Cu,h-\h/e) l+1 * + 

Cjj,eh~ 1 [h/e^) r+2 + C^eh~^ as e > /P 
C/J/Ts as h< e < 

where terms —h~ 2 Af^ and h~ 2 N^ obviously cancel one another. 
2.5.2 Successive approximations 

Now we need to get rid off the third term in the right-hand expression of 

(2.50) without condition e 3> h which allows to eliminate it. To do this we 
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need more sophisticated arguments. Namely, let us apply the successive 
approximations. 

Assuming that the original operator is perturbed by 0(/x~ 2 ) we arrive 
to the heuristic conclusion that 

(2.51) Contribution of \ t\ x T to R w " does not exceed the previous estimate 
of the contribution to R W/ multiplied by u~ 2 Th~ y \ 

this statement needs justification. We can achieve a better estimate by 
taking more terms in the successive approximation but it leads to a rather 
overcomplicated formula; so we take just one term. 
Then, in frames of (2.51) we have three cases: 

(a) e > h; then all what successive approximation does is to eliminate the 
third term in the right-hand expression of (2.50). 

(b) fi^h^ 1 < e < h; then 1 > k < e~ l and we need to adjust our arguments 
treating separately \k\ < k and \k\ > k and recalculate the second 
term in the right-hand expression of (2.50); however equator zone is 
not affected. 

(c) e < vn\n(h, fi~ 3 h~ 1 ). In this case in addition we need to recalculate the 
last term in the right-hand expression of (2.50). 

Let us reconsider the second term in the right-hand expression of 
(2.50). First, expression (2.49) with summation over k : \k\ > k re- 
turns Cfi~ 2 h~ 2 (h/e\k\) r+ 2 \k\ calculated as k = k; so we get Cfi~ 2 h~ 2 k = 
Cn~ x h- X e- X . 

Further, as \k\ < k we need to replace {h/e\k\) r+ ^ by 1; so we sum 
C fih~ 1 \k\~ 1 x /x 3 /7 1 1 A" | which returns C[i~ 2 h~ 2 k again. Therefore 8 ^ we 
conclude that 

(2.52) As h > e, r > 1 one should replace the second term in (2.50) by 

Furthermore, as e < mm(h, h^ 1 ) one should multiply the last term in 
(2.50) i.e. C/i5/)~5 by fi~ 3 e~ 1 h~ 1 resulting in C/i - ?^ 1 /)^. Therefore 8 - 1 we 
conclude that 



8 ) Pending (2.51). 
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(2.53) As e < min^,/!- 3 /)- 1 ), r > 1 one should replace the second term in 
(2.50) by Qrie^/Ts. 

Therefore we (heuristically, the proof follows) arrive to 

Proposition 2.5.1. Consider two magnetic Schrddinger operators A and A 
satisfying the standard smoothness assumptions, non- degeneracy assumptions 
(0.6) and (2.36)-(2.37) ; (2.3) and such that 

(2.54) g* = gi" + 0(fi- 2 ), V=V+ 0(fi~ 2 ), Vj = Vj + 0{^ 3 ). 
Then for r > 1 

x ^w,, ^ i.i ^ I u~ 2 h~ 2 (H/s) r+ i as e > h 

2.55 R^ r) < Cfi-'h- 1 + C 1 \ 1 ' ~ + 

qJ/Ts a s mm(h, uT 3 h^ 1 ) < s < & 



C/j, 36 1 h 3 as e < m\n(h, fr 3 h 1 ) 



and 

5 W// 



(2-56) R x 7 )7 : = 

7~ 2 | J (e T (x, X, T) - /7 _2 jV x ,corr(r) ~ e T (x, X, t) + /l~W x ,corr(r)) ^7 dx| < 

{u.- 2 h- 2 (h/8) r+1 i(h/e~iy as ~ie>h 

u.~ 2 h- 2 {h/ey +1 2 as e>h + 

/x _1 /7 _1 £ _1 as e < h 

{[i,eh~ l (h/s^) r+2 + Cu,^eh~^ as e > TP 
^i/i-i as min^,^ 3 /?- 1 ) < e < 

tu/wZe /or r = 

(2.57) R^ := |e T (x,x,r) - e T (x,x,r)| < C//- 1 ^ 



C 



H'^h'^e' 1 as fi 3 h > e^ 1 
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(2.58) R™' )7 := 7 - 2 | J(e T (x,x,T)-e T (x,x,T)^,dx\<C f i- 1 h- 1 + 
/? _1 £~ 1 7~5 as u.~ z e~ l < 7 < min(/7£ _1 , /x -1 ), 

/i/7- 1 min(l,(///7/£)^)(V £ 7)' as 7 > 

2.5.3 Justification 

To justify (2.51) and thus to prove proposition 2.5.1 we need 

Proposition 2.5.2. Consider two magnetic Schrddinger operators A and 
A satisfying (2.54). Then as \t\ < efi 3 h and |sin(r)| > max (ee|t|, Cti) 

(2.59) U{x,y,t)-U{x,y,t) = 

(47rn)- 1 /(sin(^(t)))" 1 e / ^^ (x ' y ' t) ^ b' m {x, y, t)h m 

m 

mod O^-V^llsin^t))!" 1 x (h/\ sin(t)|)') 

(2.60) |D^J < C m ^- 3 |t|r 1 |sin(^(0)r m - |/31 . 
Proof. Obviously in our assumptions 

U(x, y, t) = (47rn)- 1 /(sin(^(t)))~ 1 e / ^^ (x ' y ' t) ^ b'^(x, y, t)h m 

m 

mod O^lsin^t))!" 1 x sin(0(t))|)') 

with 6^, satisfying 

|D a b / ;|<C ma |sin(^(0)r m - H 
and plugging it into the right-hand expression of 

(2.61) e ih ~ ltA - e'*- 1 * = ih' 1 f e ihrH ' A {A - ~A)e ih ~^-^ A dt' 

Jo 
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we arrive to (2.59). Here to cover zones where either |sin(t')| < e or 
| sin(t — t')\ < e we just pass to the standard representation through non- 
singular oscillatory integrals (with one or two extra variable) and then apply 
the stationary phase method. We leave easy but tedious details to the 
reader. □ 

Corollary 2.5.3. In frames of proposition 2.5.2 equality (2.59) holds with 
U(x,y, t) and U(x,y, t) replaced respectively by 

(2.62) {4nh)- 1 i{sm{9{t)))- 1 e ih ~ 1(t>{x ^ t) ^ b m {x, y, t)h m 

m 

and 

(2.63) (47r^)- 1 /(sin(^(0))- 1 e' 7, "^' t) £ b m (x, y, t)h m . 

m 

2.5.4 Reduction to the pilot-model operator 

From now we are interested only in the estimates without micro-averaging. 

The problem however is that A at this moment is a more general operator 
than the pilot-model: it contains linear with respect to x terms in the 
original (before rescaling) g' k and V, and also it contains quadratic terms 
in V 2 = xi + \(3uxl + /3i 2 xix 2 + \ + ■■■ all of which generate 0(fi~ 1 ) 
perturbation. 

However one select a conformal coordinate system in such way that 
Vw(0) = 9) . Then if we know that F = F(0) + 0(|x| 2 ) we conclude that 
/3n = (3i2 = and by a gauge transformation we can make [3 2 2 = 0. 

Further, we can assume that r = and then make V(0) = — 1 by division 
V(0) !->■ —1, h i->- h(— V(0)) _ 2 5 11 1->- ^(_\/(o)) - 2 before rescaling. Rescaling 
/i (and changing coordinate orientation if needed we can make F(0) = 1. 

Furthermore, we can assume that F = 1. Really, as we can consider 
energy level r = we starting from equation —hD t u = Au, go to 

(2.64) F- l Au = -hF^DtU = -hD t + hF~ 1 {F - l)D t u 

9 ) Really, changing conformal system to another conformal system is done with yi = f, 
Vy2 = { — f X2 , f Xl ) which can be satisfied iff Af = 0; then u) is replaced by wlVfl and if 
f — axi + bxix 2 + jc(xj — xf ) then IV/^ — a + bx 2 + cx x + 0(|x| 2 ) and selecting a, b, c wc 
can make uj(0) — 1, Vw(0) = 0. 
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and on an interval \r\ < hh^^T^ 1 (which is completely microlocally admissi- 
ble as \ t\ x T the last term is 0({i~ 1 hh~ s T~ 1 ) and multiplying by Thr x we 
get 0(fi^ 1 h^ s ) which allows us to use the successive approximation method. 
We leave to the reader to use our standard methods and solve a rather 

easy 

Problem 2.5.4- Prove that effectively in estimates one can make 5 = 0. 

Note that or 1 < /i -3 /? -1 as // < h~i and therefore this new error (F is 
replaced by F(0)) does not exceed the right-hand expressions of (2.55) and 
(2.57) as r > 1 and r = respectively. 

Meanwhile our new error does not exceed C/z/? _1 x ^r 1 = Ch~ l anyway 
and as /i > /7~2 ( =>. e <h) it is less than u,3h~3 leave alone e"/? -1 which 
are parts of right-hand expressions of (2.55) and (2.57) respectively. Thus 
we arrive to 

(2.65) Estimates (2.55) and (2.57) hold with A which is our original pilot- 
model operator. 

2.6 Superstrong magnetic field 

Finally, let /i > h -1 . Then as only matter is what is inside of the circle of 
the radius C^iT^h^ (we can get rid off h~ 5 factor in estimates) we get in 
successive approximations C/i/? -1 x ji~ l h x ji'^h^e' 1 x Z?" 1 where \i~ x h is 
the size of the perturbation, iT^h^e^ 1 is T (without rescaling); so we get 

We leave to the reader a rather easy 

Problem 2.6.1. Prove that effectively in estimates one can make 5 = as 
well in /i"2 +,5 /j2 radius and iT^h^e^ 1 time. 

So we arrive to 

Proposition 2.6.2. Under assumptions (2.36), (2.37), F = 1 and (2.54) 
estimate 

(2.66) |e T (0,0,r) - e T (0,0,r)| < C^-^/)-^- 1 as fi > h' 1 . 

We leave to the reader yet another easy 

Problem 2.6.3. Write correction to r and may be a as F ^ const and 
gJ k ^ const. 
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2.7 Main theorem 

Now we can write our main theorem of these two sections: 

Theorem 2.7.1. For magnetic Schro dinger operator in domain X G M 2 
such that 8(0, 1) C X and satisfying conditions (0.5) -(0.7) the following 
estimates hold for x e 6(0, |) 

(i) As 1 < n < h-\ 

(2.67) |e(x, x, 0) - /T 2 A/-f (0)| < Cfi^h' 1 + C/J/H + Cf/h^ 

and 

(2.68) R™, := |e(x,x,0) - h~ 2 (J^(0) + A^ corr(r) (0)) | < 

C/i- 1 /?- 1 + C(ih-s + Cfih~ 1 ( f i 2 h) r+1 2 + 

j (h- 1 (hu. 5 2) r+1 > + as u.<h-l, 

| 5 1 2 

^ fi3h 3 as (j, > h 5. 

fnj ^4s < n < \r\ 

(2.69) R^:= 

| (e(x, X, r) - /? _2 A/" x , corr ( r) - e x (x, X, r) + /rW x , cor r(r)) | < 

- 2 h' 2 {fi 2 h) r+1 2 as fi< h~s 

i "I" 
as fi < h 2 _ 

/T 1 (/J/?) r+ * + /i^/j"! as /x < 

/- J 5 , _1 2 1 

<- ^ //3/7 3 as h $ < fi < h 2, 

_ 1 _ 4 1 

H in 3 as n> h 2 

w/w/e /or r = 

(2.70) R^ := |e(x,x,r)-e x (x,x,r)| < + 

as n< h 2 , 
as as fi < h~z ; 
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where here and in (Hi) e x is constructed for a pilot-model in x; 
(Hi) For magnetic Schrddinger-Pauli operator with p > h^ 1 

(2.71) |e(x,x,r) - e x (x,x,r)| < Cuhh'^. 

2.8 Problems 

Finally, let us formulate a series of the problems with probably simple parts 
(i) and really difficult parts (ii): 

Problem 2.8.1. As p < h^ 1 drop condition (0.5) in the pointwise asymp- 
totics: 

(i) Use the simple rescaling technique p h- >■ pp\p, h h-> /?p~i with the scaling 
function p(x) = max(eo| \/(x)|, /^/?, /z -1 ) reducing operator to the similar one 
either satisfying condition (0.5) or with ph x 1 or with p = 1 and considered 
in subsection 5.2.1; in two latter cases condition (0.5) is not needed; 

(ii) Using propagation, improve what follows from the simple rescaling tech- 
nique. 

Problem 2.8.2. Derive asymptotics without condition (0.6) replaced by the 
non- degeneracy assumption (14-1.1) i.e. \F\ + \ VF\ x 1: 

(i) Use the simple rescaling technique fi t->- fip 2 , h \-t hp^ 1 with the scaling 
function p(x) = max(e |F|, p~?); 

(ii) Improve what follows from the simple rescaling technique. 

Problem 2.8.3. Derive asymptotics under non- degeneracy assumption (0.8) 
rather than (0.7). The results of this section provide them at point x where 
|W/F| > 7 := max(/i-\//~5/ji). 

(i) Can one improve those results using that under condition (0.8) + and 
(0.8) - (when (0.8) holds but (0.8) + does not) the drift dynamic is elliptic 
elliptic and hyperbolic respectively? 

(ii) Can one improve a trivial estimate 0(/i/) _1 ) at point x where |VV/F| < 
j? Will be results affected by the difference between cases of (0.8) + and 
(0.8)"? 
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3 Dirac energy: 2D-estimates 

In this and the next sections we consider asymptotics of expression I defined 
by (0.2). 

3.1 Tauberian formula 

Let us consider first contribution of zone {|x — y\ > C7}. 

Proposition 3.1.1. Under conditions (0.5)-(0.7) the contribution of zone 
{\ x — y\ > C7}, to the remainder is 0(/i~ 1 /)~ 1 7~ K ) while the main part is 
given by the Tauberian approximation \ T , i.e. by the same expression (0.1) 
with e(x, y, 0) replaced by its standard implicit Tauberian approximation 
with T x en (0.11). 

Proof. Recall that e lh 1{A is the propagator of A and U(x, y, t) is its Schwartz' 
kernel. 

Consider expression (0.1) with u(x,y) replaced by w 7 (x, y) which is a 
cut-off of cj(x, y) in the zone {|x — y| x 7} and with the original functions 
■01, vb 2 replaced by 1. Let us replace one copy of e(x, y, r) by e(x, y, r, r') = 
(e(x, y , t) — e(x, y , r')) with r' < r and the second copy by e(x, y, t") and 
denote the resulting expression by l 7 (r, r', r"). 

Now let us use decomposition 

(3. 1) w 7 (x, y) = 7' d - K J ^i, 7 (x, z)vb 2n (y, z) dz 

with some ipi,ip2 £ ^o°°i as before subscript 7 means rescaling. 
Then l 7 (r, r', r") does not exceed 

(3.2) ^C 7 -1|^(r,r> ; || 1 

j 

where E(t,t') = E(t) — E(t'), ipj are real-valued 7-admissible functions 
supported in C 7- vicinities of zj and balls B(zj, 2Q7) cover domain X with 
the multiplicity not exceeding Co. Here we used that ||£(t")|| = 1. Since 
E(t, t') is a positive operator and tpj = ipj, one can replace the trace norm 
by the trace itself and get 



(3.3) 



Cr K Tr Vj E(t, r>j = C 7 " K Tr E(r, r')vJ, 

j 
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with if> = T.jVj- 

Further, we know from the standard theory of Chapter f 3 that under 
conditions (0.5)-(0.7) 

(3.4) || E{t,t')$ ||i < Ch- 2 (\T-r'\ + CT^h) 

Vt, t' G [— e, e], T = e/i 

and therefore 

(3.5) |l 7 (r, r', t")\ < C 7 - K /)- 2 (|r - r'\ + CT^h) 

in the same framework and therefore due to the standard Tauberian ar- 
guments we conclude that the contribution of zone { |x — y | x 7} to the 
Tauberian remainder estimate does not exceed C / u _1 /? _1 7 _K which implies 
the statement immediately. □ 

However we need to consider also zone {|x — y| < C7}, complementary 
to one above. Assume that 

(3.6) Q K (z) = DzPk-U + *Vi.o 

j 

with the first subscript at Q showing the degree of the singularity. Then 

(3.7) u K (x, x - y)ip 7 {x - y) = ^ D Xj (w K _ij^ 7 ) + u K i(/ y + uj^ip" 

j 

where if)^ = ip{{ x — y)l 1 ) with if) supported in 6(0, 1) and equal 1 in 6(0, |) 
while if)' is defined similarly with if)' supported in 6(0, 1) \ 6(0, |) and the 
last term gains 1 in the regularity. 

After integration by parts expression l Ki7 , defined by (0.10) with Q. 
replaced by Q.ij)^, becomes 




(x, y){hD Xj ) (e(x, y, r) • e(y, x, r)) dxdy 



plus two other terms: the term defined by (0.1) with the kernel Q' K j of the 
same singularity n, albeit without factor /j -1 and supported in the zone 
{\ x — y\ > an d the term defined by (0.1) with the kernel ^' K _ij, also 
without factor h^ 1 and of singularity (k — 1). 
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The former term could be considered as before yielding to the same 
remainder estimate 0(n~ 1 h l ~ dr y~ K ). To the latter term we can apply the 
same trick again and again raising power (and these terms are treated in 
the same manner (but simpler) as we deal below with (3.8). 

So, one needs to consider (3.8) and thus, denoting the second copy of 
e(y, x, r) by f(y, x, r) and without using that they are equal 

(3.9) (hD Xj )(e(x,y,r)-f(y,x,r)) = 
{hD Xj e(x, y ,r))f(y,x,r) - e(x, y, r) (f (y, x, r) \hD Xj )) = 
(Pj, x e(x,y,r))f(y,x,r) - e(x, y, r) (f (y, x, r) f P,- x ) . 

Recall that Pj = hDj — u.Vj(x) and f Pj = —hDj — fJ>Vj(x) is the dual operator. 
Also recall that if e(x,y,r) and f(y,x,r) are Schwartz kernels of E(t) and 
F(r), then Pj x e(x,y, r) and f(y,x,r) t Pj x are those of PjE(r) and F(t)Pj. 
Therefore we are interested in the expressions of the type 

(3.10) / T 1 J J oj K ^ 1 (x,y)e(x,y,r)f(x,y,T)^dxdy. 

If k < 1 then replacing e(x, y, r) and f(y, x, r) by their standard Taube- 
rian expressions one gets an error not exceeding Ch' 1 x ^ 1 h~ 2r y 1 ~ K because 
|| PjE(t) II < C , ||P//"(t)|| < Co where F(r) is an operator with the Schwartz 
kernel f(x,y, r) and also because 



(3.11) 52\\<p J PjE{T,T')<pj\\ 1 <Y,hjPjE(T,r')\\ 2 . ||E(r,T> y || 2 < 
j j 

Y^TripjPjE{ry)Pfipj + ^Tr^£(r,r> i 
j j 

< Ch- 2 (\r-r'\ + Cfi^h) Vr,rG[-e,e] 

which also easily follows from Chapter 13. 

So, in this case one gets remainder estimate 0( / u~ 1 /7 _1 7 _K + /i _1 /7 _2 7 1_K ) 
which is optimized to 0(/i _1 as 7 x h. 

On the other hand, as 1 < k < 2 one can apply the same trick again since 
we did not use the fact that e(., ., .) and f(., ., .) coincide; then we arrive to 
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the same estimates with Pj replaced by PjPk or even by P J := Pj x P} 2 ■ • • Pj,: 
(3.12) JrP J E{T,T'){P J )* < Ch- 2 {\T-r'\+ fi^h) W, r G [— e, e] . 
Finally, 

Remark 3.1.2. note that as k ^ 1, decomposition (3.6) is always possible. 
Further, as k = 1 this decomposition is possible as well provided one adds 
term x(x)|x — y| _1 with an appropriate coefficient. On the other hand, if 
k = 1 and u(x,y) = x(x)|x — y\~ l then this decomposition is also possible 
but with uoj(x,y) = x(x)(x,- — yj)\x — y| _1 log |x — y\. 

So we arrive to 

Proposition 3.1.3. Let conditions (0.5)-(0.7) be fulfilled. Then 

(i) As < k < 2 and either k ^ 1 or k = 1 anc? w(x,y) is replaced 
by u(x, y) — x(|(x + y))|x — y| _1 with an appropriate smooth coefficient 
x(x), with the error 0(fi^ 1 h^ 1 ^ K ) one can replace e(x,y, r) by its standard 
Tauberian expression (0.11) in the formula (0.11) for I. 

(ii) As k = 1 and us = x(|(x+y))|x— y| _1 , w^/i the error 0(/x 1 /? 1 K | log /?|) 
one can replace e(x,y,r) &y its standard Tauberian expression (0.11) in the 
formula (0.11) for I. 

Remark 3.1.4- (i) The arguments above show that in an appropriate sense 
one can consider arbitrary k e M and even in C. 

(ii) One needs only (3.12) rather than (0.5)-(0.7), and (3.12) holds as (0.7) 
is replaced by a weaker non-degeneracy condition (0.8) + ; 

(hi) Furthermore (3.12) with an extra factor (1 + u,h\ log h\) in the right-hand 
expression holds under condition (0.8). 

Problem 3.1.5. Can one prove the similar result for \ m defined by (1.0.6) 
V. Ivrii [Ivr3] with m > 3 ? 

Thus we arrive to 

Proposition 3.1.6. Let conditions (0.5) and (0.6) be fulfilled. 

(i) Further, let either condition (0.8) + be fulfilled or condition (0.8) be fulfilled 
and fi < log Then (0.11) and statements (i), (ii) of proposition 

3.1.3 hold. 
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(ii) On the other hand, let condition (0.8) be fulfilled and /7 _1 | I og | 1 < 
\i < h^ 1 . Then (0.11) and statements (i), (ii) of proposition 3.1.3 hold with 
an extra factor (1 + fih\ log h\) in the right-hand expressions. 

Remark 3.1.7. Under certain assumptions (see [IR04]) this result could be 
generalized for d > 4. However in calculations we use that d = 2. 

3.2 Superstrong magnetic field 

Consider Schrodinger-Pauli operator as fih > 6q. Then we arrive immediately 
Proposition 3.2.1. Let fxh > €q. Then 

(i) Contribution of {\x — y\ > 7} to I does not exceed Cfih~ 1/ ~f~ K ; 

(ii) Further, 

(3.13) |i| < cy+H-^k 

(Hi) Furthermore, I = 0(/i _o °) provided 

(3.14) V + /jhF>e Vx 6 6(0,1); 
in particular it is the case as 3 < 1 and fih > Q. 

Proof. Statement (i) trivially follows from the fact that S£ 2 norm of e(.,,,r) 
does not exceed C[i*h~i. 

Meanwhile contribution of {(x, y) : x G 6(0, 1), y G 6(0, 1), |x — y| < 7} 
does not exceed C/i 2 /? _2 7 2_K . Adding to (i) we conclude that the sum 
reaches its minimal Cfi 1+ 2 K h~ 1 ~2 K value as 7 = /i _ 5/?2. 

Finally, statement (hi) is due to the fact that e(x,y,0) = 0(/i~°°) as 
condition (3.14) is fulfilled in one of points x,y; see subsection 13.5.1. □ 

Proposition 3.2.2. Let conditions (0.1)-(0.5) be fulfilled. Let fih > e and 

one of the nondegeneracy conditions 

(3.15) \j + (2m + l)fj,h\ + > e 

or 

. , V V V 

(3.17) + |— + (2m + l)fih\ + |V — | < e =^ detHess— >e 
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(3.16) |l -l T | < Cuh K h~ l 2 



holds as k 7^ 1; as k = 1 this estimate holds with an extra factor | log /x| in 
its right-hand expression; 

(Hi) If instead of (3.15) or (3.17) + we assume that 



then (i), (ii) hold with an extra factor | logyu| in the right-hand expressions; 

(iv) I — l T = 0(/i~°°) under ellipticity condition \ V + (2m + l)fihF\ > e. In 
particular, it happens as 1 < 3 2Z + 1 and jih > C . 

Proof. To prove statement (i) recall that the drift speed does not exceed 
C/x -1 and therefore Hilbert- Schmidt norm of vbE(r)vb' does not exceed C as 
vb,vb' are ^ °°-functions with dist(supp vb, supp vb') x 1. 

Really, it is true for a Hilbert-Schmidt norm of (E(t) — E(r'))^ with 
|r — t'\ < iT x h and then by Tauberian theorem it is true for a Hilbert- 
Schmidt norm of (E(t) — E T (r))^' with E T operator with the Schwartz 
kernel e T with time T x [i. However vbE T vb' is negligible as T < efi due to 
propagation results. 

Then rescaling x i->- x/7, y i->- y/7, h i->- /1 4 117 we conclude that 
for 7-admissible vb,vb' with dist(supp -0, supp vb') x 1 Hilbert-Schmidt norm 
of rbE(r)rb' does not exceed C and therefore the contribution of 

K(z) = {(x, y) :\x- z\< 27, \y - z\ < 2j, \x - y \ > 7} 

does not exceed 7~ K ; therefore taking partition we conclude that the con- 
tribution of zone {(x,y) G 6(0, 1) x 6(0, 1), |x — y| > 7} does not exceed 



Statement (iv) follows from the fact e — e T = 0(/i 00 ) under ellipticity 
condition. 



(3.17) 




C 7 ' 



-2-k 
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To prove statement (ii) we apply the same arguments as in proposition 3.1.3. 
While (3.1)-(3.3), (3.6)-(3.10) remain true, (3.4), (3.5) should be modified: 
factor h~ 2 must be replaced by \i\r x . 
We claim that 

(3.18) As < /t < 1 (3.11) holds with factor br 2 replaced by /^i Ai — i . 

Really, the same arguments are applied as before albeit now PjE(t) is 
bounded by C (/i/?)2 rather than by C . 

Then we arrive to the remainder estimate 0(^i~ K + //2/7~2^ 1_K ) and 
optimizing by 7 = 7 := C/j,~^h^ we arrive to the remainder estimate 

As k = 1 we get the same estimate albeit with a logarithmic factor. 
To tackle the case 1 < k < 2 one should replace (3.12) by 

(3.12)' Tr P J E{r, r'){P J )* < C^-\\t - r'| + a,' 1 ti)(u,ti)\ J \ 

Vr.r'e [-e,e] 

as fih > 1; one can prove it easily as ||P j £(t)|| < C(l + /i/?)^- 7 !. 

Then as 1 < k < 2 we arrive to the remainder estimate 0(7 _K + / u/7 _1 7 2_K ) 
and optimizing by 7 = 7 we again arrive to the remainder estimate 0(7~ K ). 

Finally, the above arguments imply (iii). □ 

4 Dirac energy: 2D-calculations 

4.1 Pilot Model 
4.1.1 Transformations 

Consider pilot-model operator (1.1). Let us rescale as before 10 -*. Then as 
U(x,y, t) is defined by (1.9)— (1. 10), and e(x,y, r) is given by (1.13), and t 

lf) ) We need to add factor ^~ 4+K = /i~ 2 x pT 2 x pT K where two factors fi~ 2 are coming 
from dx and dy and p K from lj. 
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means a complex conjugation we get 

(4.1) I := (27r)-V 4+K // J J ' oo^x, ^y){t 't'T'x 

l -{2it)- A ?T 2 ti K JJ J J cu{x,y)csc{t')csc{t")(t't")- 1 x 

exp^/TT^fx.y, t') - 0(x,y, t") - (t' - t») c/t / c/t // c/xc/y 

with 

(4.2) 0(x,y,t')-0(x,y,t")-(t / -t> = 

-l(cot(t')-cot(t"))(xi-yi) 2 
- 1 cot(t / )(x 2 - y 2 + 2t'e) 2 + J cot(t")(x 2 - y 2 + 2t" £ ) 2 

+ (X! + yi + 2e)(t' - t")e - (f - t")(r + e 2 ). 

Recall that in the original coordinates u>(x, y) = fi(|(x + y); x — y) where 
.) is uniformly smooth with respect to the first variable and positively 
homogeneous of degree — n with respect to the second one as |x — y| < 1. 
Therefore without any loss of the generality one can assume that 

(4.3) u{x, y) = Q(^(x + y); x - y)^(x - y) 

where now Q is positively homogeneous of degree — k and ip G ^ °°(B(0, 1)) 
equal 1 in 8(0, |); the difference would be of the same nature but with k 
replaced by 0. 

We can replace variables x,y with new variables x := |(x + y) and 
z := (x — y) and rescale. Note that (4.2) depends on x in a very specific 
way: it does not depend on x 2 at all and it is linear with respect to xi; thus 
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(4.1) becomes 11 ^ 

^(27r)- 3 rV +2 J J J I ft^/T^r' - t"),z) csc(r / )csc(r")(r / r")- 1 x 

exp^TT^-^cotfO - cot( t"))z 2 - ^cot{t'){z 2 + 2t'a f i- 1 ) 2 + 
l - cot(t")(z 2 + 2r" £ ) 2 - (t' - t")(r - e 2 ))) dt'dt'^^z) d Zl dz 2 

where ft(A; z) = F Xl ^ A ft is a partial Fourier transform. Rewriting the third 
line in (4.2) as 

- - (cot( t) - cot( t')) (x 2 - y 2 + ( t + t')ef 

- l -{cot{t)-cot{t')){t-t'fe 2 - 

^ (cot(t) + cot(t')) (x 2 - y 2 + (t + t>)(t - t> 
and replacing t' = t + s, t" = t — s we arrive to 
(4.4) 1 = 

i(27r)- 3 /rV +2 // // ^(4rVs,z)csc(t + s)csc(t-s)(t + s)- 1 (t-s)- 1 x 



exp(/7* 1 -^(cot(t + s) - cot(t - s)) (z 2 + (z 2 + 2te) 2 + 4s 2 e 2 ) 
- (cot(t + s) + cot(t - s))se - 2s(r - e 2 ) ) dtds i)(uT 1 z)dz l dz 2 . 

Remark J^.l.l. (i) For a sake of simplicity we assume (without any loss of 
the generalization) that a > 0. 

(ii) If a <C 1 in the above analysis applied to the general operators we 
need to take a-admissible with respect to x function ft Q instead of ft; thus 
up to the shift we replace ft(x, •) by ft(xa _1 , •) and thus Q(2h~ 1 es, •) by 
a 2 Q(2h~ 1 eas, •). In this case we will use notation \ a instead of I. 

In the same time we need to consider separately |z| < fia and |z| > fia. 
For a while we will not assume that e x fi~ l a. 



n ) In these calculations we skip X2 as an argument of Q and integration by dfi X X2 
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4.1.2 Case a 2 ^> fih 

Assume first that 

(4.5) a 2 > fih 1 ' 5 

where for a sake of simplicity we assume that a > (and therefore e > 0). 

Remark J^.1.2. (i) In the virtue of the factor Q.(2h~ 1 a 2 s, •) under assumption 

(4.5) we need to consider only |s| < h 6 and therefore we can consider 
separately \t'\ < e , \ t"\ < e and \ t'\ > e , \ t"\ > eo; 

(ii) Note that due to section 6.3 contribution of zone {\t'\ < eo, \ t"\ < e } de- 
fined by integral expressions (4.1) or (4.4) with an extra factor Xe (O 
or Xe {t") or Xeo(^) differs from the same expression for non-magnetic 
Schrodinger operator by 0(/i/) _1_K x h K ) = 0(/i K+1 /? _1 ) as k / 1 and 
by 0(/i 2 /? _1 | log/i|) as k = 1. 

(iii) Furthermore, if we remove from this expression for a non-magnetic 
Schrodinger operator cut-off {\t'\ > eo} then the error would not exceed the 
same expression as well. 

Let us consider contribution of zone {\t'\ > eo, \ t"\ > e } defined by an 
integral expressions (4.1) or (4.4) with an extra factor (l — x eo (t')). Due to 
remark 4.1.2(i) we need to consider only t' , t" belonging to the same tick. 

Let us consider first zone 

(4.6) {\s\ > h/a 2 , |sin(t)| > C\s\}. 

Then cot(t') — cot(f") x — sin _2 (t)s and integration by parts with respect 
to z delivers one of the factors 

(4.7) kr^N 2 + l z 2l ■ l z 2 + ^1)^1 sin(t)| 2 , 

(4.8) |s|- 1 % 1 | 2 + |z 2 + et| 2 ) _1 |sin(t)| 2 . 

Thus integrating by parts many times in the zone where both of these factors 
are less than 1 we acquire factors 

(4.9) (l + M/r^Izi! 2 + |z 2 | • |z 2 + et\)\ sin(t)| 2 ) 

(4.10) (l+ |s|r 1 (|z 1 | 2 + |z 2 + £t| 2 )|sin(t)r 2 ) 
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respectively. Multiplying by |z| K and integrating we get after multiplication 
by |sin(t)|- 2 

(4.11) x |sin(t)|- K |s|- 1+ HH« 

and integrating by t over one tick intersected with {t : | sin(t) | > |s|} we 

get 

k < 1, 

(4.12) Cti-^l \s\^ K K >1, 




s\ 2 (1 + I log |s||) K=l. 



This expression (4.12) must be either integrated by with respect to s: 
|s| < h/a 2 or multiplied by (h/a 2 ) 1 ^^ 1 due to factor ft and integrated over 
|s| > h/a 2 , resulting in both cases in the same answer which is the value of 
(4.12) x |s| calculated as s = h/a 2 i.e. 

{{h/a 2 )i K «<1, 
(Va 2 ) 1 -^ «>1, 
{h/a 2 )Hl+ \\og{h/a 2 )\) k=1. 

In addition to zone (4.6) we need to consider zone 

(4.14) {|sin(t')|x|s|, |sin(t")|<|s|}; 

its tween {| sin(t')| x |s|, | sin(t")| < |s|} is considered in the same way. 

In zone (4.14) | cot(t') - cot(t")| x | sin(t // ) | 1 and in this case (4.7), (4.8) 
are replaced by 

(4.7) ' n(|z 1 | 2 + |z 2 |-|z 2 + ^|)- 1 |sin(0|, 

(4.8) ' n(|z 1 | 2 + |z 2 + 5t| 2 )" 1 |sin(t // )| 
and (4.9), (4.10) by 

(4.9) ' (l + ^(N 2 + |z 2 | • |z 2 + et\)\ sin(Or 1 ) _/ , 

(4.10) ' (l + r 1 ^! 2 + |z 2 + et| 2 )| sinC^)! X ) 
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Then, multiplying by \z\ K and integrating we get after multiplication by 
| sinCt'') I 1 1 

(4.11) ' n 1 -2 K |sin(t")r5 K |sr 1 ; 

then integrating by \t"\ over one tick but intersected with {| sin(t")| < |s|} 
we get 

(4.12) ' ftH K |s|-k 

Finally, either integrating over |s| < h/a 2 or multiplying by |s|~'(/i/ a 2 ) 1 
and integrating over |s| > h/a 2 we get in both cases the same answer 
h 1 ~2 K, (ha 2 ) 1 ~2 K not exceeding (4.13). 

Therefore the total contribution of zones (4.6) and (4.14) is given by 
expression (4.13). Then multiplying by \k\~ 2 u. K h~ 2 a 2 we get after summation 
with respect to k : \k\ > 1 the value as k = 1 i.e. (4.13) x fi K h~ 2 a 2 . 

Therefore we arrive to 

Proposition 4.1.3. For the pilot-model operator with 1 < u. < hr x and 

a > toji^ 1 under additional assumption (4.5) 

(4.15) \\l-\l'\<CR w (a): = 

{/i 1+re /rV- re 0<k<1, 

u. 2 h- K a K 1<k<2, 

/i 2 /?- 1 a(l + |log(/i/?/« 2 )|) « = 1 

where ^ is a Tauberian expression for I albeit with T = tQ^T 1 . 

Now we need to explore the difference between I J' and defined by 
(0.9). 

Proposition 4.1.4. For the pilot-model operator with 1 < /i < hr x and 

a > eQjjb 



(4.16) I'a-^CI < C^h~ K a 2 . 

Proof. Repeating arguments of subsection 6.3.4 one can prove easily that 

(4.17) - £ ^)j,J\ < C(fih) L h- 2 - K a 2 + C^h-'a 2 

/:0</<Z.-lJ=0,l 
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where Z^-^i is defined by (0.9) albeit with ft(|(x + y). x — y) multiplied 
by a homogeneous polynomial of degree (/ +j) with respect to (x — y) and 
for I = j = this polynomial is 1 and with integrals taken only over zone 
{M:\x-y\<^}. 

Since without any loss of the generality we can assume that ft is even 
with respect to the second argument, all terms with odd (/ +j) vanish and 
since replacing u, by — u, we should arrive to the same result, all terms with 
odd / vanish and therefore 

(4.18) In (4.17) all terms vanish except those with j = and even /. 

Picking L = 2 we arrive then to (4.16). □ 

Combining propositions 4.1.3 and 4.1.4 and noting that the right-hand 
expression of (4.16) does not exceed CR w (a) we arrive to 

Corollary 4.1.5. For the pilot-model operator with 1 < u, < h^ 1 and 

a > eo/i -1 under additional assumption (4.5) 

(4.19) \\l-l™\<CR w (a). 



4.1.3 Improvement 

In a certain case (under assumption (4.22) below) we can improve the results 
of the previous subsubsection 4.1.2. To do this we note that (an easy proof 
is left to the reader) 

Remark 4-1-6. Let £ < e\t\ with 

(4.20) C:=(VM)*|sin(t)|. 

Then the sum of (4.9) and (4.10), multiplied by \z\~ K dz and integrated, 
does not exceed 

U 2 (e\t\y« 0<k<1, 

(4.21) C I ( 4 - 2K (e\t\) K - 2 1<k<2, 

{( 2 (e\t\y\l + \\og((/e\t\)) K = l. 
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We need to multiply by | sin( t) | 2 dt and integrate. Then as < k < 1 
main contribution came from sin(t) x 1 and improvement is possible iff then 
( < e (as \k\ = 1 provided the large part of contribution). But as s x h/a 2 , 
sin(t) x 1 (which was the main contributor in the previous subsubsection) 
( x a while e = fi~ 1 a. 

So, improvement is possible only for 1 < k < 2. In this case main 
contribution came from sin(t) x s and improvement is possible if then 
( x (ft|s|)5 < fi~ 1 a. Setting \s\ = h/a 2 we have ( x h/a and therefore 
improvement is possible if 

(4.22) a 2 > fi 2 h. 

Consider zone (4.6) first. Then multiplication by | sin ( t) | — 2 dt and inte- 
gration results instead of (4.12) in 



(4.23) 



h 2 \s\ 2 £ 



c < 



h 



1—K I _| 1 — K,_K — 2 



/p|s|-^(l + |log(ft|s|/£ 2 )) 



3 
2' 

< k < 2, 

3 

" 2 



1 < K < 

3 
2 

K 



where we set \k\ = 1. Finally, either integrating over \s\ < h/a 2 or mul- 
tiplying by \s\~ ! {h/ a 2 ) 1 and integrating over \s\ > h/a 2 we get instead of 
(4.13) 



(4.24) 



/j, K ha K 



C I 



3 
2' 

< k < 2, 
3 

~" 2' 



1 < K < 

3 
2 

K 



H~2ha~2(l + | log( / u 2 /?/a 2 )|) 
Consider zone (4.14); we apply remark 4.1.6 albeit with 
(4.20)' C:=(ft|sin(t")l)*- 

Again multiplying by | sin(t / )| 1 1 sin(t // )| 1 and integrating by \t"\ over one 
tick but intersected with {|sin(t // )| < |s|} and setting |sin(t')| = |s|, and 
finally either integrating over |s| < h/a 2 or multiplying by |s| _/ (^/a 2 )' and 
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integrating over |s| > h/a 2 , we arrive to the same answer, not exceeding 
(4.24). 

Finally, multiplying by \k\~ 2 fi K h~ 2 a 2 we get after summation with respect 
to k : \k\ > 1 the following improvement to proposition 4.1.3 

Proposition 4.1.7. For the pilot-model operator with 1 < // < h^ 1 and 
a > max( / u~ 1 , fihi , /jihi -5 ) 



(4.25) |l^-iy|<C/?r(«) 



2kl-1 2-k 



1 < K < 



2' 



^6-2^2-2^-4 



/ u 3 /T 1 oi(l + llog^ 2 /?/" 2 )!) 



K 



< K < 2, 

3 
2' 



Remark 4-1-8. Right-hand expression of (4.25) in comparison with (4.15) 
has factor 

3 
2' 

< k < 2, 

3 
2' 



(4.26) 



tfh/a 2 )*- 1 
(V 2 h/a 2 ) 2 - 



^ 2 h/a 2 )Hl + \\og(^h)\) 



1 < K < 

3 
2 

ft 



While we cannot directly apply estimate (4.16) now as fi 2 h K may be 
greater than R™(a), we can apply (4.17)-(4.18) with L = 4 resulting in 

Proposition 4.1.9. For the pilot-model operator with 1 < /i < h' 1 and 
a > max(/i _1 , /i/?5 , ^2h^ & ) 



(4.27) 



If -Z 

I a a 



W 



2^1 < CR?{a) 



W/ 



mt/i X^ rr := Z^ Q but taken only over zone {\x — y\ < u, x }. 



4.1.4 Case a 2 ^> fih 

Assume now that a > ^T 1 but (4.5) fails. Then in contrast to the previous 
we will need to compute contributions of pair of ticks (k', k") with k' ^ k". 
Then if t', t" belong to k'-th and k"-th ticks respectively we denote r = k' — k" 
and s=t'-t"- 2vrr. 
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Contribution of k' = k" 7^ 0. Consider case k' — k" = r = first. As 
a 2 > h applying the same arguments as before we get the same answer 
(4.13) as before. Meanwhile, as a 2 < H we need to integrate (4.12) over 
{s : |s| < 1} and we arrive to /t 1_ 2 K for all k : < k < 2. 

Multiplication by fj, K h~ 2 \k\~ 2 a 2 and summation with respect to k results 
in the right hand expression of (4.15) as a 2 > h and 

(4.28) C/i K /r 2 ftH« a 2 = Q> +1 /rH« a 2 

as a 2 < h. 



Contribution of 7^ k' 7^ k" 7^ 0. Consider now case k' — k" = r 7^ 0. As 

a 2 > h expression (4.13) acquires factor (\k' — k"\a 2 /h)~ l and multiplying by 
fi K h~ 2 a 2 \k'\~ 1 \k"\~ 1 we get after summation with respect to k' 7^ 0, k" 7^ 
the same right-hand expression of (4.15) 12 \ 

Meanwhile, as a 2 < h instead of (4.13) we get 

c/H^i + i*' - k"\a 2 /hy'. 

Multiplication by ^h^a^k'^ 1 ^"^ 1 and summation with with respect to 
k', k" returns 

(4.29) Cii 1+1 2 K h~ l ~^a 2 {l + (log(ft/a 2 ))+) 2 . 

Contribution of = k' 7^ k". Consider now k' = 0, k" 7^ 0; its tween 
case k' 7^ 0, k" = is addressed in the same way. To do this we need to 
modify our expression fo I: namely, we do not divide by t' and then take 
Fourier transform F t /_^-i r ; instead we take Fourier transform F t /^-v, then 
integrate by r' to r and divide by h; modifying this way (4.1) we arrive 
instead of (4.4) to 

(4.30) \ a = 

^(27r)- 3 rV +2 « 2 J jj JJn{Ah~ 1 a 2 s,z)x 
csc(t + s) csc(t - s)(t + s) _1 (t - s) _1 x 

exp^/TT 1 -^(cot(t + s) - cot(t - s)) (z 2 + (z 2 + 2tef + 4s V) 
-(cot(t+s)+cot(t-s))se-2s(T-s 2 )+(t+s)(T-T') ) dtdsip(fi~ 1 z)dz 1 dz 2 dr. 



12 ) Actually, with an extra factor (H/a 2 ) 1 but we do not need it here. 
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Then repeating our above arguments we conclude that 

(4.31) As a 2 > h the resulting term (instead of (4.13)) does not exceed 
(4.13) x {h/a 2 \k"\) 1 x fr 1 and as a 2 < h the resulting term does not exceed 
Ch 1 - 1 ^ x (l + a^/c"^- 1 ) - ' x h' 1 

and we can rewrite the result in both cases as 

oH^i + o^i/c"^ 1 )-'. 

Multiplying by ji K h^ 2 \k"\^ 1 a 2 we get after summation with respect to 
k" ^ 

(4.32) QiH-MV^/a 2 )' as a 2 > h, 

(4.33) C^H- 2 -5 K a 2 (l + (log(ft/a 2 ))+) as a 2 < h. 

So we have proven 

Proposition 4.1.10. For the pilot-model operator with 1 < /j, < h~ x , a > 

(4.34) I fa) I < (4.15) + (4.32) as a 2 > h, 

(4.35) \\J a) -Xf a) \ < (4.28) + (4.29) + (4.33) as a 2 < h. 

4.2 General operators 

Consider now general operators satisfying (0.5)-(0.6) and either (0.7) or 
(0.8) or (0.8) + , r = 0. 

First, using proposition 3.1.1, remark 3.1. 4(h), (ii) and propagation results 
we conclude that 

Proposition 4.2.1. Let conditions (0.5)-(0.6) be fulfilled. Then 

(i) Contribution of zone {|x — y| > Cq^ 1 } to l T does not exceed Cii^^hr 1 
and 

(ii) Under condition (0.7) contribution of time {\t\ > C } (before rescaling) 
to l T does not exceed Cji^h" 1 . 



Chapter 4. Dirac energy: 2D -calculations 



63 



Note that 0(/i K+1 /? _1 ) is smaller than the remainder estimate /? w (l) or 
R^(l) and as long as we are interested in Weyl approximation X w (even 
with correction term) we should be completely happy with it. 

Under condition (0.7) we are done, immediately arriving to statement 
(i) and related part of (iii) of theorem 4.2.3 below. 

Under condition (0.8) let us introduce 



and consider covering of 8(0, 1) by a-admissible elements. In virtue of 
proposition 4.2.1(i) we need to consider only pairs (x,y) belonging to the 
same element. Due to rescaling, contribution of one such element to |/ T — X w | 
does not exceed Cfih~ 1 ~ K a 2 and therefore 

(4.37) Under condition (0.8) contribution of zone 

{(x,y) : \x-y\ < C Q n~\ \VVF~ l \ < a} 

to |/ T — X w | does not exceed Cfih~ 1 " K a 2 . 

Therefore we need to consider only zone where |WF _1 (x)| > a and 
\ x — y\ < Cofi^ 1 . We are going to prove that 

Proposition 4.2.2. Let conditions (0.5) -(0.6) be fulfilled. Further, let 
fi < h -1 . Then 

(i) For each z with a(z) > a contribution of zone 



to |I T — X w | does not exceed R w (a) which is defined as the right-hand 
expression of (4.15) or (4.25). 

(ii) Further, as a 2 > fi 2 h, 1 < k < 2 this contribution does not exceed /?J /V (a) 
which is defined as the right-hand expression of (4.25). 

Note that both R w (a) and R^ J (a) contain a in the positive powers and 
thus after integration with respect a -1 da we get their values as a = 1. 

Therefore we arrive to statement (ii) and related part of (iii) of theo- 
rem 4.2.3: 



(4.36) 




a=Comax(/x 1 ,(/ih)^ 5 ) 



{(x,y)eB(z,a(z)), |x - y| < C // _1 |} 
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Theorem 4.2.3. Let conditions (0.5) -(0.6) be fulfilled. Then as /i < h 1 
(i) Under condition (0.7) 

(4.38) |l -X w | < C^/T^H- R w 

with 

{//^/T 1 as k <1, 

u 2 h~ K as k > 1, 

H 2 h~ 1 \ log(yu/?)| as /c = 1. 

Under condition (0.8) 13 ^ 

(4.40) |l - X w | < C f i- 1 h- 1 - K + /? w + Cfi 2 h- K {fih)- 5 ; 

(Hi) As fi 2 h < 1 anc? 1 < k < 2 m t/iese estimates one can replace X w fey 
X w +X™ r and /? w 6y 

fi"- 2K h 2 - 2K \<k<2, 



(4.41) /?* 



w ._ 



/^(l + iiogOu 2 /))!) 



2 
3 

K= 2- 



Remark 4-2.4- (i) Actually as fi 2 h < 1, k > | this correction term does not 
exceed Ca~ 1 h~ l ~ K . 

(ii) Consider other terms in the estimates. Then (4.28), (4.29), (4.33) contain 
a in the positive powers and (4.32) contains a in the negative power and thus 
after integration with respect a -1 da we get their values as a = (/i/?)? i.e. 
(4.28), (4.29) become Cah K+2 h~^ K which is less than (4.39) and (4.41) while 
(4.32), (4.33) become C/z2 K+1 /)~ 1 ~2 K which is larger than (4.39). Therefore 
selecting a > (fih)2~ 5 saves us from all these problems for the price of 5 ^ 0. 

Even as h~ 1+s < /i < h^ 1 it works as C/i2 K+1 /7 _1_ 2 K acquires factor (fih) 1 
which makes it subordinate to (4.39). 

Problem 4.2.5. Prove above results with 5 = 0. 



13 ) If (0.8) + fails the first term should acquire factor (1 + fih\ log h\) but it is important 
only as fi is close to h~ x but then the other terms dominate. 
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Proof of proposition 4-2.2. First of all applying proposition 2.2.2 and the 
same arguments as in the analysis of section 2 we arrive to 

(4.42) Contribution of zone {| sin(0(t'))| > C h, |sin(0(t"))| > C h} to the 
error |I T — X w | does not exceed expression (4.15) or (4.41) (in its frames). 

Therefore we need to explore contributions of three remaining zones 

(4.43) {| sin(0(t'))| < 2C h, | sin(0(t"))| < 2C h} 
and 

(4.44) {| sin(0(t'))| > 2C Q h, \ sin(0(t"))| < C Q h} 

(the same analysis would cover its tween zone { | sin (0(t')) \ >2C h, |sin(0(t"))| < 

C h}). 

Note first that the standard propagation arguments imply that we get 
a factor C{\z-y\ + |z 2 |) _, (| s\n(9(t'))\ + e\t'\)' and the similar factor with t' 
replaced by t". 

Zone (4.43) as e > Qh. Then in the classical propagation shift e\k'\ is 
always greater than H. 

Obviously contribution of zone (4.43), intersected with (k' , k") windings 
x{(x, y) : |x — y\ > p} to the Tauberian expression does not exceed 

(4.45) Cfi 2 h 2 \k'\- 1 \k"\- 1 p- K sup JJ \U{x,y,t)\ 2 dxdy < Cfi 2 h 2 - d p- K a d 

where factor h 2 is due to integration with respect to t' , t" and U means 
cut-off of U to the energy level \r\ < c. 
Plugging (as d = 2) 

(4.46) p = e^ 1 max(|/c'|,|/c // |), 
e = p~ x a we get 

(4.47) C// 2K+2 |/f / |- 1 |^ / r 1 « 2 " K (max(|/c'|,|/c // |))" K 



Summation by k', k" returns Cp 2K+2 a 2 K which is less than the right-hand 
expression of (4.15) or (as a 2 > p 2 h, 1 < k < 2) the right-hand expression 
of (4.25). 
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Meanwhile, estimating U(x,y, t) by h d and integrating by |x — y] < p 
and then integrating by t', t" over \9(t') — irk'\ < h, \9(t") — nk"\ < H we get 

(4.48) Ch 2 - 2 W~ K I k'\~ 1 1 /c" | - 1 (/i/p) « d 

where the last factor is due to propagation results; plugging d = 2, (4.46) 
we get 

(4.49) C/j 2 h- K \k'\- 1 \k ,, \- 1 a 2 (h/emax(\k / \, \k"\))'; 

after summation we get Cn 2 h~ K (h/e) 'a 2 which is again less than the right- 
hand expression of (4.15) or (as a 2 > fi 2 h, 1 < k < 2) the right-hand 
expression of (4.25). 

Zone (4.43) as e < Qh. Taking sum of (4.47) over k', k" : max(|/c'|, \k"\) > 

h/e we get 

Cfi 2K+2 a 2 - K (l + | \og fi 2 h/a\) x (fi 2 h/a)' K = Cfi 2 h- K a 2 {l + | lo g/ u 2 /?/«|) 

while taking sum of (4.49) we get the right-hand expression as well. Obviously, 
it is less than the right-hand expression of (4.15). 

Thus we need to consider only k', k" : max(|/c'|, \k"\) < h/e. 

Using the same arguments as before, we conclude that the contribution 
of zone (4.43), intersected with (k',k") windings x{(x,y) : |x — y| > h} 
to the Tauberian expression does not exceed Cfi 2 h~ K \k'\~ 1 \k"\~ 1 ; similarly 
contribution of zone (4.43), intersected with (k', k") windings x{(x,y) : |x — 
y| < h} to the Tauberian expression does not exceed Cn 2 h~ K \k'\~ 1 \k"\~ 1 a 2 
as well. 

After summation we get Cfx 2 h~ K (l + | log p 2 h/a\) 2 which again is less 
than the right-hand expression of (4.15). 

Zone (AAA) as e > Qh. Now shift with respect to (x— y) is pr x max(£|t|, \et 
sin(t')l); then repeating analysis in zone (4.43) we arrive to (4.47) for 
1 < k < 2 and to (4.49) for < k < 2; as < k < 1 one should replace 
(4.47) by 

(4.50) Cu. 2K+2 h\k'\- 1 \k"\- 1 

with a logarithmic factor as k = 1 and summation with respect to k', k" 
results in C/x 3k+2 /?(1 + | log/?|) 2 which again is less than R w (a). 

Zone (AAA) as e < C-Ji. With the above adjustment arguments of zone 
(4.43) work here as well. □ 
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4.3 Perturbations 
4.3.1 General scheme 

Contribution of k' = k" ^ 0. Consider first case k' = k" = k. Assume 
first that a 2 > h. 

Recall that /c-th tick' contribution to l T was not exceeding expression 
(4.13) multiplied by fx K h~ 2 a 2 \k\~ 2 i.e. 

{{h/a 2 ) 1 ^ k < 1, 

{h/a 2 ) 1 - 1 ^ k > 1, 

(V« 2 )^(l + |l0g(n/« 2 )|) K=l 

(as | sin(r')| < h or | sin(r")| < ft one should use proposition 4.2.2). Then we 
can suspect that the successive approximations bringing factor aT/h as U 
is replaced by U° simultaneously brings the same factor when we replace 
U by U in I; so contribution of /c-th tick to the error is estimated by the 



previous expression multiplied by Ca\k\/h i.e. 

{{h/a 2 )^ K «<1, 

{h/a 2 ) 1 - 1 ^ k > 1, 

(ft/a 2 )^(l + |log(/i/a 2 )|) k = 1. 
and summation over /c : \k\ < k := /l/cr results in 

fih/a 2 )^ k < 1, 

(ft/a 2 ) 1 -^ k > 1, 



[ (n/a 2 )5(l + | log(V« 2 )l) « = 1. 

Remark 4-3.1. (i) Here a is the magnitude of perturbation and as we approx- 
imate by pilot-model it is 0(fi~ 2 ) but if we are trying to use the magnetic 
Weyl approximation it is 0(e). 

(ii) Expression (4.53) provides us with the answer as a 2 > h; as a 2 < K we 
need to recall that if we are interested only in zone {|x — y\ < C yU -1 } before 
rescaling then |z| < Co and we need therefore there reset h/a 2 to 1 and also 
\s\ to 1 resulting in 

(4.54) Cair K+2 a 2 \\ogh/a\h- 2 -^ K 
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(iii) In these settings we need to assume that T < C s 1 i.e. k < C s 1 i.e. 

(4.55) fxh < cr/e. 

This inequality always holds as a = e but it may fail as a = ji~ 2 ; in this 
case we either reset k to Ce^ 1 or reevaluate a and analyze separately in 
both cases (see in the next subsubsection). 

(iv) Note, as fih < a/e summation (4.51) for k : \k\ > k brings the same 
expression albeit without logarithmic factor; otherwise if we reset k = Ce^ 1 
this sum returns Cu. K h~ l . 

(iv) However for the pilot-model approximation we reexamine zone {\x — y\ > 
Cofi^ 1 } and improve the remainder estimate of (iv) 

Contribution of ^ k' ^ k" ^ 0. Next we need to consider ^ k' ^ 
k" 7^ 0. Recall that as a 2 > h the contribution of such pair to l T does not 
exceed 

(4.56) CpL K+2 a 2 \k'\- 1 \k"\-\h/a 2 \k' - /^'D'/T^x 

{h/a 2 )^ K k < 1, 

{h/a 2 ) 1 - 1 ^ k > 1, 

(h/a 2 )Hl + \\og(h/a 2 )\) « = 1; 

multiplying by Ca\k'\/h (as perturbation goes to the first "factor") we get 

(4.57) Ca^ +2 a 2 \k"\- 1 (H/a 2 \k' - k"\)'fr 2 -* K x 

{h/a 2 ) 1 ^ k < 1, 

{h/a 2 ) 1 - 1 ^ k > 1, 

(V« 2 )kl + |log(V« 2 )|) « = 1. 

Then summation with respect to k' of [h/a 2 \k' — k"\)' returns (h/a 2 ) 1 as 
h < a 2 and summation with respect to k" then returns (4.53) with an extra 
factor (h/a 2 ) 1 and therefore is less than (4.53). 

On the other hand, as h > a summation with respect to k' of 1 returns 
h/a 2 and as we recall that big left brace expression should be replaced by 1 
we arrive after summation with respect to k" to 

(4.58) Ca^ +2 h- l -^ K {l + (log h/a)) 
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which is (4.54) with an extra factor (h/a 2 ) and therefore is greater than 
(4.53). 

Contribution of = k" ^ k! . If k! ^ 0, k" = we recall that l/c"^ 1 
should be replaced by fT 1 in (4.56)-(4.57) and perturbation brings factor 
a\k'\/h so we get instead of (4.57) 

{h/a 2 )^ K «<1, 
C^ K+2 a 2 (ft/a 2 |/c'|)'r 3 ^ <J (n/a 2 ) 1 -^ « > 1, 

(ft/a 2 )5(l + |log(ft/a 2 )|) « = 1. 

Then summation with respect to k' returns 

(4.59) C<T f j, K+ W{h/a?)'tr 3 -2 K 
as h < a 2 and 

(4.60) Cafi K+2 h- 2 -^ K 
as h > a 2 ; (4.60) is even greater than (4.58). 



Contribution of = k! . With perturbation factors are no more of the 
same rights and we need to consider this case as well. 

(i) If k' = 0, k" ^ we recall that {k'l^ 1 should be replaced by fT 1 in 
(4.56)-(4.57) but perturbation brings factor a so together we get factor ahr 1 
which falls into into case k' = 1, k" ^ 0. 

(ii) Similarly, if k' = k" = we fall into case k' = 1, k" = 0. 

In total we get (4.53)+(4.59) as a 2 > h and (4.60) as a 2 < h. So, we 
arrive to 

Proposition 4.3.2. Let conditions (0.5) -(0.6) be fulfilled in 6(0, 1) and let 
a := e |W/F(x)| > Q/i" 1 where x e 6(0, \). 

Let e z (x, y, r) be a Schwartz kernel of the spectral projector of either the 
magnetic Weyl approximation (then a = air 1 ) or the pilot-model approxi- 
mation (then a = iT 2 ) at point z. Let T > min(eo^/(J, Co/ia" 1 ) > 1. 

Let^M e^°(6(0,l)). 
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(i) Then as a 2 > K, 7 = C /i 1 

(4.61) R 2 = R 2 (a, 7 ): = 

I J Q a , 7 (^(x + y),x-y)(|e T (x,y,0)| 2 -|e7 (x+y) (x,y,0)| 2 )c/xc/y|< 

(4.53) + (4.59) 

where 

Q Q , 7 (x, z) := vb(a~ 1 (x -x))^i(7~ 1 ^); 

(ii) As a 2 < h, 7 = Co/i -1 we get 

(4.62) R 2 < (4.54) + (4.60). 

fw,) Further, in this inequality one can replace e T by e. 

Proof. The easy albeit tedious proof following our standard technique we 
leave to the reader. □ 

4.3.2 Improvement 

We want to apply arguments of the previous subsubsection "Improvement" 
to improve results of the previous subsubsection. One can see easily that 
improvement makes sense only if we can improve k' = k" = 1 and thus as 
a 2 > ji 2 h. This excludes magnetic Weyl approximation. 

In this case R^(a) defined by (4.25) acquires factor /i -3 /? -1 : 



(4.63) C < 



3 

,,2k-3l-2 2-k 1 ^ ^ ^ 

^3-2^1-2^-4 ^<K<2, 

h- 2 aHl + \\og(fi 2 h/a 2 )\) k=\. 



This replaces (4.53) and (4.59) becomes negligible. 
So, we arrive to 

Proposition 4.3.3. Let conditions of proposition 4-3.2 be fulfilled, we con- 
sider the pilot-model approximation and let a 2 > fi 2 h. Then 

(i) R2(«,7) does not exceed (4.63). 

(ii) Further, in this inequality one can replace e T by e. 
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Proof. The easy albeit tedious proof following our standard technique we 
leave to the reader. □ 

4.3.3 Results. I 

Let condition (0.7) be fulfilled. 

Magnetic Weyl approximation. If we are looking for the magnetic 
Weyl approximation, we should assume that h~2 < /i < br 1 and pick up 
<7 = e = a = 1 and look what happens with expressions (4.53) and 
(4.59). The former returns the last term in the estimate (4.64) below and 
the latter returns the lesser expression C^i K ~ 1 h~ 2 h l . So we arrive to estimate 
(4.64) of theorem 4.3.4 below albeit with an extra factor (l + (logii 2 /7) + ) in 
front of the big left brace. 

To get rid off it let us consider two-term successive approximation. Then 
to estimate a remainder one needs to sum the same expressions as before 
but with an extra factor which is either fi~ 2 h~ 1 \k'\ or /x 2 /? 1 1 /c" | . 

This really does not matter for k' = k" ^ and we sum u~ 2 hr x instead 
of \k\~ x from \k\ = 1 to \k\ = fi 2 h resulting in 1 rather than logarithm. As 
^ k' ^ k" 7^ we get the same answer albeit with a factor (fih) 1 as a x 1. 
If k! = and/or k" = we replace \k'\~ l and/or \k"\~ 1 by FT 1 but this is 
more than compensated by (/i/))' but effectively T' = h and/or T" = h so 
in fact we reset these cases to k' = 1 and/or k" = 1. 

On the other hand, the second term in the successive approximation 
refers to perturbation being strictly a(xi — yi) but this leads to as ft is 
even with respect to (x — y). Thus we arrive to estimate (4.64) of theorem 

4.3.4 below. 

Theorem 4.3.4. Let conditions (0.5)-(0.7) be fulfilled. Then as br \ < u. < 
h- 1 

{^h- 2 < k < 1, 

h- 1 -* 1< k < 2, 

/)- 2 (l+|l0g(/i/))|) K = l. 

Pilot-model approximation. If we are looking for the pilot-model ap- 
proximation, we should assume that h~* < fi < h^ 1 and pick up a = ji~ 2 , 
e = a = 1. We need to consider variants: 
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(i) h~3 < fi< h~*. As < k < 1 (4.53) becomes C/i K - 2 /r 2 (l + | \og p 3 h\) 
(with another logarithmic factor as k = 1 and it is larger than Cp K h~ 1 
(which is the contribution of {|x — y| > CqP^ 1 }). 

As 1 < k < 2 we can replace (4.53) by (4.63) which is then 0(u.~ 1 h~ 1 ~ K ) 
and p K h~ 1 is also 0(p' 1 h' 1 ' K ). 

Noting that (4.59) is negligible we arrive to statement (i) of theorem 
4.3.6 below. 

(ii) < p < h~ x . In this case we note that (4.59) is still less than (4.53) 
but p K h^ 1 may become leading term. To improve this term let us consider 
contribution of zone |x — y| x 7 > Cq^T 1 into an error: we need to take 
a = 7 2 and T x ^7 so we get the estimate of perturbation 

(4.65) Cfih- 1 7"- V K x a T/h x C7T 2 7 2 ~ K 
recall that unperturbed term does not exceed 

(4.66) Ciih- l T- 1 1 - K x C//- 1 /7- 1 7 - K_1 

because the drift propagation does not expand B(x, p) with p > C yU _1 
but merely shifts it. We sum (4.65) as 7 2 7~ < h i.e. 7 s < pT l h i.e. 

7 < 7 = /i~3/)3 = yU _1 (/i 2 /))3. 

4 1 12 

Then we get Ch~3~3 K p3 K ~3 and we arrive to statement (ii) of theorem 
4.3.6 below. 

Remark 4-3.5. Obviously arguments of the magnetic Weyl approximation 
allowing us to get rid off logarithmic factor still work albeit the second term 
does not necessarily vanish as the perturbation is quadratic. The logarithmic 
factor in estimates (4.67) and (4.69) does not seem be worth of trouble to 
write this second term. 

Theorem 4.3.6. Let conditions (0.5)-(0.7) be fulfilled. Then 
(1) As h~l < p < /T5 

(4.67) |l -T| < Cp-H- 1 -^ 

( p K ' 2 h' 2 
C\\ogp 3 h\ I 

[ A^- 1 /7- 2 (l + j l0g/7|) 



< K < 1, 

1 < K < 2, 
K = 1 
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where 

(4-68) 7:= |Q 1 , 1 (^(x + y),x-y)|e. (x+y) (x,y,0)| 2 c/xc/y 
(ii) As h~* < \i < h' 1 

(4.69) |l -T| < C/i^/T 1 ^ + Ch- 4 3- 1 3 K / I 1 3 K - 2 3 + 

{fi K ~ 2 h- 2 
/.-^(l + iiogC^)!) 

4.3.4 Results. II 

Assume now that condition (0.8) is fulfilled instead of (0.7) and estimate an 
approximation error again. 

Magnetic Weyl approximation Let us consider the magnetic Weyl 
approximation first. Then we should assume that h~2 < /i < Z?" 1 and pick 
up a = e = a/i -1 . Recall that the contribution of {(x,y) : |x — y| > C yU -1 } 
does not exceed Cfi^h^ 1 and therefore we need to consider zone 

(4.70) {(x, y):\x-y\< C ^\ a{x) > C^ 1 } 

where also (1 - e) < a{y)/a{x) < (1 + e) M \ 

Therefore we can use a-admissible partition in zone (4.69) and we need 
just look what happens with expressions (4.53) and (4.59) after integration 
with respect to a~ 1 da as a > (/i/?) 5 and (4.60) /i -1 < a < {[rh)*. 

Obviously integral of (4.53) resets to its value as a = 1 as under assump- 
tion (0.7); integrals of (4.59) and (4.60) reset to their values as a = (/i/?) 5 
which are equal and less than what we got already. 

Therefore estimate (4.64) still holds 15 ) and we arrive to statement (i) 
of theorem 4.3.7 below after we apply the same arguments with two-term 
approximation as in the previous subsubsection to get rid off the logarithmic 
factor. 

14 ) And add contribution of the zone {(x,y) : a(x) < a :— Ci/U _1 , a(y) < a} which 
would be C/.ih~ 1 ~ K x a 2 = C[i~ 1 h~ 1 ~ K in virtue (0.8). 

15 ) Albeit the first term Cfi^ 1 h^ 1 ^ K acquires factor | log h\ if assumption (0.8) + fails: 
see Tauberian estimate. 



< K < 1, 

1 < K < 2, 
K = 1. 
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Pilot-model approximation Let us consider the pilot-model approxi- 
mation; then again o = ji~ 2 . Then we need to consider two cases. 

(a) < /x < In this case we are completely happy with the zone 
{(x,y) : \x-y\ > C /x -1 }. 

As < k < 1 we need to integrate with respect to a' 1 da (4.53) and 
(4.59) as a > (fih)2 and (4.60) as a < (/i/j)?. The first integral resets 
to (4.53) as a = 1; integrals of (4.59) and (4.60) reset to their values as 
a = (fih)2 which are equal (up to a logarithmic factor) and less than what 
we got already. Thus we arrive to statement (ii) of theorem 4.3.7 below. 

As 1 < k < 2 we need to replace integral of (4.53) by its integral over 
zone (/i/j) 5 < ol < fiti? plus integral of (4.63) over zone fiti? < a < 1. 
These integrals are reset to their integrands values as a = (/i/)) 5 and a = 1 
respectively; the latter is larger than the former albeit not necessarily larger 
than integral of (4.60). 

Thus we arrive to statement (iii) of theorem 4.3.7 below. 

(b) < /j, < h -1 . In this case the contribution of {(x,y) : |x — y| < Cq/x -1 } 
does not exceed the right-hand expression of (4.69) but we need to reexamine 
zone {(x,y) : |x — y| > C yU -1 } which we split into two: {|x — y| < e o;(x)} 
and {|x — y| > eo«(x)}. 

Repeating corresponding arguments of the analysis under assumption 
(0.7) we conclude that the contribution of the former zone does not exceed 
what we got then. In the latter zone T x u 2 and its contribution does not 
exceed Cfi~ 1 h~ 1 ~ K . 

Thus we arrive to statement (iv) of theorem 4.3.7 below. 

Theorem 4.3.7. Let conditions (0.5)-(0.6) and (0.8) be fulfilled. Then 15 *) 

(i) As < ji < /r 1 estimate (4.64) holds; 

(ii) As h~i < /i < h~i , < k < 1 estimate (4.67) holds; 

(iii) As h~\ < fi < 1 < k < 2 

(4.71) |l - 7| < Cfi- 1 h- 1 - K + Cah K - 2 h-^- 2 {l + | logii 3 /?|); 

(iv) As h~i < ji < h^ 1 estimate (4.69) holds. 
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4.4 Superstrong magnetic field 

Consider under assumption (0.8) the Schrodinger-Pauli operator with u, > 
br 1 . Recall that according to propositions 3.2.1(h), 3.2.2(h) respectively 

ill < cy+H-H* 

and 

|l - l T | < Cuh K h-^ 

(as k ^ 1; as k = 1 the latter estimate holds with an extra factor | log /x| in 
its right-hand expression). 

Let us assume that F = 1 to avoid some unpleasant correction terms as 
we reduce operator to one with F = 1. Recall that now h := a~^h^ is the 
magnitude of the cyclotron radius. 

Magnetic Weyl approximation. Consider hrst the magnetic Weyl ap- 
proximation. Then a = ah = a/i~2/?2 and T = a _1 h = a~ 1 fi2ti2- 1 then the 
the shift (unrescaled) due to magnetic drift does not exceed Can~ 2 T = Ch, 
which justifies the choice of a. 

Therefore as a > a = C\u~2\)2 the contribution of a-ball intersected 
with {(x,y) : |x — y| < 7 := Qj/i^s/js} to |I T — X MW | does not exceed 
Cfih~ 1 a x fi2 K h~2 K x Ca / U2 + 2 K /7~2~2 K . Then summation with respect to all 
such balls does not exceed C / U2 + 2 K /?~2~2 K ) as in estimate (4.72) below. 

Meanwhile the contributions to both I and T of zone {(x, y) : |x — y\ > 7} 
does not exceed C^/~ 1 ~ K which is also of the same magnitude. 

Finally, the contributions to both I and I of zone {(x,y) : a(x) < 

a(y) < a} do not exceed C/i 1+ 2 K /? -1- 2 K x a 2 x C/i2 K /j~2 K which is less 
than the right-hand expressions in both (4.72) and (4.73). 

Thus we arrive to estimate (4.72) below. 

Pilot-model approximation. Meanwhile for the pilot-model approxima- 
tion a = a~ x h as long as T < 7"* := a~ 1 ^h2- 1 so we conclude that as a > a, 
the contribution to (I — I) of a-ball intersected with {(x,y) : |x — y| < a} 
does not exceed Ca[j,2 K h~2 K . After summation with respect to a~ 1 da we 
get Cu.2 K h~2 K which is less than the right-hand expression in (4.73). 

Now we need to consider larger T and 7; then a = 7 2 and T x a^u 2 ^. 
Then we should dehne 7 so that aT < h iff 7 < 7. Therefore we pick up 

111 — 412 _ 222 

7 = yU~3/73Q;3 and matching T = u.^h^oT^ and 7 = /i~3/j3Q;3. 
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Using the same arguments as above we conclude that the contribution to 
(I — I) of the zone {(x,y) : a < \x — y \ < 7} and contributions to both I and 
I of the zone {(x, y) : a > \x — y\ < 7} also do not exceed the right-hand 
expression in (4.73). 

Thus we arrive to estimate (4.73) below. 

Theorem 4.4.1. For Schrddinger-Pauli operator under conditions (0.5)- 
(0.6), (0.8), n > /r 1 and F = 1 

(4.72) ||-J MW | < C^+H~H« 
and 

(4.73) |l -T| < Cuh +1 i K h-h-h. 

4.5 Problems and remarks 

Remark 4-5.1. The main difference between cases 1 < k < 2 and < k < 1 
as fih -C 1 is that in the former case the main contribution to the remainder 
is delivered by (x, y) close to one another (|x — y | fi^ 1 ) while in the latter 
case by (x, y) with |x — y | x 

We can calculate X MW and I plugging corresponding expressions e^ w (x, y, 0) 
and e z (x, y, 0) with z = |(x + y) into I. 

Problem 4.5.2. Find nice expressions for X MW and\. 

Problem 4.5.3. As fih < 1 get rid off condition (0.5); according to Chap- 
ter 13 we do not need it for estimate |l — l T | but we want to get rid of it in 
estimates for \ I — X w | , 1 1 — X MW | and \ I — 1 1 . To do this 

(i) Assume first that condition 

(4.74) |V| + |W|xl 

is fulfilled and consider the scaling function a(x) = max(e| V(x)|, \ih, C/^ 1 ) . 

(ii) Assume then that condition 

(4.75) \V\ + |W| + |detHess\/| x 1 
is fulfilled and consider the scaling function 

a{x) = max(e{\V(x)\ + \VV{x)\ 2 )K (iih)K CfT 1 ). 
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Problem 4.5.4. Get rid off condition (0.6) assuming instead that \F\ + 
|VF| > e. 

(i) Repeating arguments of section 3 and using results of Chapter 14 one can 
prove easily that \ \ — \ T \ < Cfi~ 1 / 2 h~ 1 ~ K as ji < hr 1 . Derive similar estimates 
for Schrddinger and Schrddinger-Pauli operators as h^ 1 < fi < br 2 and for 
Schrddinger-Pauli operator as fi > h~ 2 ; 

(ii) Introducing scaling function £(x) = e|F(x)| in zone {x : £(x) > /I-*} one 
can rescale fi h> fi£ 2 , h >->■ Mr 1 , fih" 1 ^ K 1— > fi~ y br 1 ~ ti £~ 1 + K , after integration 
over £~ 2 d£ we get Cfi~ 1 h~ 1 ~ K £~ 2+K calculated as £ = fi~2 i.e. fi~2 K h^ 1 ^ K . 
This is an estimate for the contribution of the regular zone. 

Prove the same estimate for the contribution of the singular zone as 
fi < h . Derive similar estimates as h^ 1 < fi < br 2 and as fi > hr 2 ; 

(Hi) Is it possible to upgrade the above estimate to Cfi^^h^ 1 ^^ as k < 1? 
Explore also cases h' 1 < fi < br 2 and fi > br 2 . 



5 Pointwise asymptotics: 3D-pilot-model 
5.1 Pilot-model in R 3 : propagator 

Consider the pilot-model operator 

(5.1) A = A:= h 2 D 2 + {hD 2 - fix^ 2 + h 2 D 2 + 2a Xl + 2/3x 3 

which is the sum of 2D-pilot-model (1-1) A(2) and ID Schrddinger operator 

(5.2) B := h 2 D 2 + /3x 3 
and therefore propagator of A is 

(5.3) U{x,y,t) = U (2) {x',y',t)U w {x 3 ,y 3 ,t) 

where U(2)( x ' , y' , t) is the Schwartz kernel of the propagator e lh A P) and 
therefore after reseating it is defined by (1.1) while U(i)(x 3 , y 3 , t) is the 
Schwartz kernel of the propagator e' h lfB and then one can prove easily that 
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before rescaling 

(5.4) U {1) (x 3 ,y 3 ,t) = 

(27T/))- 1 J exp(// ) - 1 ((c + fit)* - (C - pt)ys - K 2 t - \p 2 ?)) d( = 

\{2vht)-h expjV 1 ^ +y 3 ) + \t-\x,-y,f - ^ 2 t 3 )); 
in particular 

(5.5) U {1) (x 3 ,x 3 , t) = l -{2itht)--> exp(/^ 1 (2/3tx 3 - \^))- 
Therefore 

(5.6) After standard rescaling x \— > fix, t i— >■ fit U(x,x, t) in comparison 
with U(2){x' ,x' , t) acquires factor 

(5.7) i/i(27r^)^exp(/r 1 (2/3tx 3 -^- 2 /3 2 t 3 )) 

and therefore (1.23) and (1.24) are replaced respectively by 

2 

(5.8) (p(t) := -t 2 fT 2 a 2 cot(t) + tfi~ 2 a 2 + -fi~ 2 /3 2 t 3 - tr 
and 

(5.9) fi- 2 a 2 (t 2 - t sin(2t)) + 2^ 2 /3 2 t 2 sin 2 (t) = (r - fT 2 a 2 ) sin 2 (t); 
the latter equation could be rewritten as (1.24) with r replaced by t' 

(5.10) fi- 2 a 2 {t 2 - tsin(2t)) = (r - fi~ 2 a 2 ) sin 2 (t) 
and 

(5.11) 2 / x- 2 /3 2 t 2 = (r-r'); 

the former equation shows the return times of 2D-movement on the energy 
level r' and the former the looping time of ID-movement associated with B 
(after rescaling) on the energy level (r — r'). 
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Remark 5.1.1. We should set T = eofi (in contrast to 2D-case). Really, 
unless we want to make assumptions outside of the ball 6(0, £) with a small 
constant £ we need to take T < e (before rescaling) and then T < 
(after rescaling). This absolves us from the analysis of the equatorial zone. 

Let us compare solutions t* k of (5.9) and solutions t k to (1.24). One can 
see easily that 

(5.12) As \k\ < e|/3|"V 

(5.13) t* k = t k (l + 0(fi- 2 p 2 k 2 )) 
and 

(5.14) <p"(tl) = ^ 2) (t k ){l + 0(f,~ 2 (3 2 k 2 )) 

with (/?( 2 ) defined by (1.23). Furthermore other properties of <f( 2 ) are fulfilled 
for ip as well. 

5.2 Tauberian estimates 

Then as \ k\ < e/i _1 in virtue of (5.6) and calculations of section 1 contribution 
of /c-th tick to F t _^ h -i T V x U does not exceed 



(5.15) C 



rih~ y x (^/h\k\y as \<\k\< k, 

Hh-\n 2 h/a\k\)* x (n/h\k\)* as k < \k\ < e/x _1 , 



where as before k = e^ 2 ha~ x . Recall that without an extra factor (////? |/c|) 2 
we would have estimate for F t _+ h -i T Y X 'U(2)- 

However, this extra factor is a game changer. As before there are three 
cases: 

(a) a > fj^h; then we set k = 1 and the summation with respect to 
k : 1 < \k\ < e\il returns 

Ciih~\ii 2 h/a)^ ■ (n/h) 1 * |og(^) x Cfih^a- 1 * log(^); 
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(b) fih£ 1 < a < fi 2 h; then 1 < k < efi£ and summation with respect to 
k : k < \k\ < e\x£ returns 

C/J^oH \og{fi£)/k x Cfih^a-^l + (\og{a£/ fih)) + ) 

while summation with respect to k : 1 < k\ < k returns the same 
expression albeit without logarithmic factor; 

(c) a < nh£~ x \ then k > efi£ and we reset set k = e\i£ and summation with 
respect to k returns 



Qi/T 1 x (fi/h^k 1 * x Qi 2 /rM. 



So we have proved 



Proposition 5.2.1. (cf. proposition 1.2.4)- (i) After rescaling as T = e\x£, 
£ > Cofi- 1 

(5.16) \F t _ h -, T xr{t)VM < C/2h- 2 + 

{i^-bT x a^ (l + (\og(a£/ fih)) + ) as a£>fih, 

fi 2 h~h^ as a£ < fih; 

(ii) Therefore as operator A coincides with the pilot-model operator (5.1) in 
B(0,£) 



(5.17) |e(0, 0, r) - e T (0, 0, r)| < Ch~ 2 £~ 1 + 

3 1 

/i/?~2£-2 as a£ < \ih; 



C 



3 5/i _1 a-5(l + (log(a£//i/?))+)r 1 as a£>fih, 



(Hi) In particular, as a x £ x 1 the right-hand expression in (5.17) does not 
exceed Ch~ 2 + C/iS/r^l + | \ogfih\) which is 0(h~ 2 ) as /i < (h\ log 



5.3 Weyl estimates 

Even simpler are Weyl estimates: in comparison with 2D-case they acquire 
factor /!2/?~2 and (1.64) becomes 

(5.18) := |e T (x,x,0)-/)- 3 AC /v | < CfTh~^ + C^a'h' 1 

and more generally we arrive to statement (i) below 
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Proposition 5.3.1. (cf. proposition 1.3.2). For a pilot-model operator 
(5.1) with r x 1, 

(i) As fi 2 h < a < 1 estimate (5.18) holds; moreover if we replace Af^ by 
■A/^ + A/^corrfr) with the correction term h~ 3 Af X:CO „( r ) delivered by r-term 
stationary phase approximation then an error does not exceed 

(5.19) Cfi-h-l + Ciih- 3 2(ii 2 h/a) r+l , 

(ii) As a < fi 2 h estimate holds: 

(5.20) R™ := |e T (x,x,0) - h~ 3 Af™\ < Qj/ri; 

(Hi) In particular, as /i < we have R w = 0(h~ 2 ) without any assump- 
tions to a, (3: \a\ < 1, \j3\ < 1. 

Recall that here we do not need to analyze the equatorial zone and 

2 2 

therefore the difference between cases a > fih~3 and a < fih~3 disappears 
together with corresponding terms in (1.66) or (1.67). 

5.4 Micro- averaging 

Micro-averaging adds yet another layer of complexity. Let us consider an 
isotropic micro- averaging with function ?/> 7 (x) = ifj(x/ / ~f) and an anisotropic 
micro- averaging with function ip-y(x) = ip{x' X3/73) where 7 = (7,73) is 
a scale with respect to (x',X3). 

Then one needs to integrate by parts with respect to X3 taking in account 
factor exp(2//?~ 1 /3tx 3 ) in (5.5) and rescaling with respect to t which brings 
factor (/i/)/|/3|7 3 |/c|) / as {jxhj |/3|7 3 |/c|) < 1 in addition to factor (/ih/\a\^\k\y 
as (pih/\a\-f\k\) < 1. Thus 

(5.21) Micro-averaging brings factor (/x/j/z/(7)|/c|)' as (fih/u(-y)\k\) < 1 with 
z/( 7 ) : = | a |7+ |/3|7 3 . 

Then we need to modify (1.53) replacing \a\^ by ^(7): 

(5.22) \k\>k(j):=fih/u(j). 

To understand how it affects a Tauberian estimate we need to consider 
the following cases: 
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(a) Micro-averaging has no effect as k(-y) > efil i.e. as 1/(7) < hi 1 ; 

(b) Micro-averaging has a minimal effect (only decreasing the logarithmic 
factor) as efil > k(-y) > k. Obviously in this case k = e max(/z 2 /j/a:, l). 

So this case holds iff hi -1 < 1/(7) < min(/i/), a/i -1 ). Then the logarith- 
mic factor becomes log(/c(7)//c). 

(c) Micro-averaging kills logarithmic factor and further leads to k^ replaced 
by £(7)2 as 1 < k(j) < k. Obviously in this case k = e min (fi 2 h/a, fit) . 

So this case holds iff max(/?£~ 1 , ctfT 1 ) < 1/(7) < fih. Then the right- 
hand expression in (5.24) becomes C/z 2 /? _1 z/~2^~2; 

(d) Micro-averaging resets k to 1 and brings factor (fih/u^)) 1 as £(7) < 1 
i.e. as 1/(7) > /i/j. Then the right-hand expression in (5.24) becomes 
C^h-l^h/v^)) 1 . 

Thus we arrive to 

Proposition 5.4.1. For the pilot-model operator (5.1) with r x 1, < 
Oi < 1, |/3 1 < 1, /i < h' 1 in M 3 estimates 

(5.23) l-V\Ft^nM(xAt) ~ Xt)nU^))\ < Qx/? T ( 7 ) 
and 

(5.24) r^lFt^MxAtnU^l < Cfih- 2 + C>/? T ( 7 ) 
hold as T = efil where 

(5.25) /? T ( 7 ) = 

' /J/7 _1 aH(l + (| g(a£//i/)))+) as al>fih, 1/(7) < W -1 , 

fih~^l^ as a£ < fih, 2/(7) < 

/J/7 _1 aT^ (1 + (log(o;//iz/(7))) + ) as M" 1 < 1/(7) < min(/i/?, a/i -1 ), 

/i/?~ 1 z/(7) _ 5 as max(/?£~ 1 , a/i" 1 ) < 2/(7) < fih, 

fi^h~^{fih/v{^f)) ! as 2/(7) > fih. 

In the right-hand expression of (5.25) two first lines correspond to the 
case (a). 
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Corollary 5.4.2. In frames of proposition 5.4-1 the following estimate holds 

(5.26) 7 " 2 7 3 ~ 1 |r(e(., ., r) - e T (., .,r))^\ < C(/T 2 + R^j))^ 

Consider now how micro-averaging affects Weyl estimate (with the 
right-hand expression (5.19)). There is no effect as ^(7) < \ih and factor 
(/i/?/z/(7))' comes out as 2/(7) > jih. Thus we arrive to 

Proposition 5.4.3. For the pilot-model operator (5.1) with r x 1, < 
a < 1, \i < h' 1 in R 3 



i) As fi 2 h < a < 1 estimate holds: 

(r)( 



(5-27) ( 7 ) 



l~ 2 li l \ f (e T (x,x,0)-/)- 3 (AC /v +Ar x , corr(r) ))^ 7 c/x| < 

Ch + C^h 2(/x h a) < v 

II as z/(7j < /xn 

fnj As a < fi 2 h estimate holds: 

(5.28) R?( 7 ) := 7 - 2 7 3 - 1 | / (e T (x,x,0) - h~^)^dx\ < 

Ch - 2 + c »h-l{^ h/u[l))l as 

1 1 as 2/(7) < /x/?. 

5.5 Strong and superstrong magnetic field 
5.5.1 Tauberian estimate 

Consider strong /ih x 1 and superstrong /ih > 1 magnetic field for the 
Schrodinger-Pauli pilot-model operator 



(5.29) A = A := h 2 D{ + (hD 2 - /ixi) 2 + h 2 D 2 + 2ax x + 2(3x 3 - /i/? 

which alternatively means that we assume that |r — | < C in the 
framework of the standard pilot-model (5.1), 3 = 1. 
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Proposition 5.5.1. For a pilot-model operator (5.1) with fih > 1 in IR 3 as 

\ x \ < Qh |y| < Cq, r < cfih 

(i) e(x, y, r) = mod 0(fih°°) as r < fih — eo (lower spectral gap); 

(ii) As a > i/ie following estimates hold 

(5.30) |e(x,x,r + /?) - e(x,y,r)| < 

C^/T 1 max(max(| \/(x) + 2y>/? - r|, /?))~^ 

and 

(5.31) \F t ^ h -i T U{x,y,t)\ < C /i 2 h- 1 max(max{\V{x) + 2j [ih - r\, h))' 1 * 

and these estimate are sharp as x = y and r is close to Landau level fih. 
Recall that in (5.31) t is rescaled. 

Proof. Statement (i) is obvious. To prove (ii) note that 

(5.32) e(x,y,r) = y e e (x 3 , y 3 , r - r') d T ie [2) {x', y' , r') = 

J e (2 )(x', y', r - t') d T/ e B (x 3 , y 3 , r) 

and therefore 

(5.33) e(x, y, r + h) — e(x, y, r) = 

y (ee(^3, 73, r - r' + h) - e e (x 3 , y 3 , r - r')) 9 r /e (2) (x / , y', r) c/r'. 

Recall that due to our analysis in section 1 d T 'e^){x', y' , t') is essentially 
supported in Cfi^h^ a- vicinity of Landau level and fast decays out of it. 
Finally recall that due to the analysis of subsubsections 5.2.1-5.2.1 

(5.34) |e s (x 3 ,y 3 ,r + h) - e B (x 3 ,y 3 ,r)| < C(r')^ as T > > Qh. 

Combining it with estimate (1.73) for |9 r /e(2)(x', y', t')\ we conclude that 
the left-hand expression of (5.30) does not exceed 

C/Ja: -1 /)^) x ayT^h^ x max(max(| V(x) + 2j/ih — r|, /?)) 2 . 
Estimate (5.31) is proven; estimate (5.31) follows from it. □ 
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Further, as before T* = tofii. Then we immediately arrive to 

Corollary 5.5.2. For a self-adjoint operator in domain X, B(0,£) C X C 
M. 2 , £ > Coh, coinciding in B(0,£) with the pilot-model operator (5.29) with 
t x I, < a < 1, \0\ < I, fi > h' 1 

(i) Statements (i)-(ii) of proposition 5.5.1 remain true; 

(ii) Formula (5.32) holds modulo 0(/i/?~i£ _1 ). 

5.5.2 Micro-averaging 

Let us consider micro-averaging. First let us estimate Fourier transform 
where as before we rescaled t h-> fit; formulae (5.16)-(5.17) imply immedi- 
ately 

Proposition 5.5.3. For a pilot-model operator (5.29) in M 3 with \xh > 1, 
73 > h 

(5.35) rVlFt^rXAtnU^l < C^h-^p 
Corollary 5.5.4. In frames of corollary 5.5.2 as i > 73 > h 

(5.36) 7" 2 7 3 - 1 |r(e(., ., r) - e T (., ., r))^| < C^/T^r 1 . 

5.6 Magnetic Weyl approximation 

We can try to use a host of the different approximations but restrict ourselves 
now to the magnetic Weyl approximation. Recall (5.32) and note that 

(5.37) e MW (x,y,r) = J e^(x 3 ,y 3 , r - t') d T ,e™? {x> ,y> , r') = 

J e^(x',y',r-r')d Tl e^(x 3 ,y 3 ,r'). 
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5.6.1 Pointwise asymptotics 

Proposition 5.6.1. For a pilot-model operator (5.1) with r < c, /ih < 1 
and for a pilot-model operator (5.29) with r < c, /ih > 1 

/,o^ i mw/ x mw/ /M^r/ //i/7_1 as fih<l, 

(5.38) |e (x,y,r)-e (x,y,r)|<C<^ _ 3 

[ /i/? 2 as fih > 1 

as \t — t'\ < h. 

Proof. Obviously the left-hand expression does not exceed 

with J = C max((/i/?) _1 , 1) which implies (5.38). □ 

Therefore one can hardly expect that the magnetic Weyl approximation 
provides a better error than the right hand expression of (5.38). 

Proposition 5.6.2. Let \r\ < e. Then for a pilot-model operator (5.1) as 
\xh < 1 and for a pilot-model operator (5.1) as fih > 1 

(5.39) R MW := \e T {x,x,0)- h- 3 Af MW \< Ch- 2 + Cfih- 1 + Cfih'l\f3\l 

Proof. Without any loss of the generality we can assume that x = 0, (5 > 0. 
As contribution of k-th tick to the Tauberian expression is 0(/i5 h~2 1 Ac | 5 ) , 
its contribution to the error when we replace (3 by does not exceed 
Qz2/)~2 1 /c I 2 x /i~ 3 /) -1 /3 2 |/c| 3 and summation with respect to 

(5.40) k : |/f| < /fi :=ji/H/i3 
returns 

(5.41) CAt-i/7-i^ 2 |/f|i = Cnh~*fi 

3 3 3 ~ 3 3 ~ — - 

On the other hand, summation of C/!2/7~5 |/c|~5 as |/c| > /ci returns C[i^h~^k x 2 
which is the same expression because /<i was defined from /i~ 3 /? _1 /3 2 K 3 = 1. 

Meanwhile, contribution of /c-th tick to the error when we replace a by 
does not exceed C^h~2 1 A" | 2 x fi~ 2 h~ 1 a\k\ and summation with respect to 

(5.42) k : \k\ < k 2 ■= f/ha- 1 
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returns 

(5.43) Cfi-h-la~kl = C^h^aS 

which is less than Cbr 2 . On the other hand, summation of C/J/7~l ]/c|~l as 

~ 3 3 ~ 1 ~ 

\k\ > k 2 returns Cfi2h~2k 2 2 which is the same expression because k 2 was 
defined from equation ji~ 2 h~ 1 ak 2 = 1. □ 

Remark 5. 6. 3. (i) We summed with respect to all k while it would be enough 
only with respect to k : \k\ < efi. However it provides us by no improvement. 

(ii) So far we have not used factor C(fx 2 h/a\k\)^; using it we acquire in (5.41) 
factor (fx 2 h/aki)2 = ^2/73 (3ict~2 resulting in 

(5.44) R MW < C/T 2 + Crib' 1 + 
provided 

(5.45) a>/i/?i/?5. 

(iii) In contrast to 2D the magnetic Weyl approximation now is better than 
Weyl approximation as fi > h~^ (when it matters) but Weyl approximation 
with the correction terms may be better still. 

5.6.2 Micro-averaging 

Consider now micro-averaging. We need to redo only the first step of the 
proof of proposition 5.6.2 and only in the case ii > h~3f3~ 2 and it is useful 
only if ki > k(-y) = fih/is(-y) or equivalently 

(5.46) = oq + /?7 3 > /3§ /?i . 

In this case in the left-hand expression of (5.41) should be reset to k\ replaced 
by k(~7) as k(-y) > 1 or by 1 with an extra factor (iih/v(~f))' as k(-y) < 1 
resulting in C/? 1 MW ( 7 ) with 

(5.47) /?! MW (7) := 

/3 2 z/(7)~2 (fih/v^)) 1 as 2/(7) > fih, a<fi 2 h, 

/3 2 z/(7) _ 5 as z/( 7 ) < ^h, ^(7) > a, 

_5 1 , 

/3 2 z/(7) 2 (j^h/a) 2 (/i/?/z/(7)) as > fih, a>fi 2 h, 

/3 2 z/(7)~5( ;(i //(^)/ a )5 as //(^) < ^h, /iz/(7) < a 
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and therefore we arrive to 

Proposition 5.6.4. Let |r| < e and (5.46) be fulfilled. Then for a pilot- 
model operator (5.1) as fih < 1 and for a pilot-model operator (5.1) as 
t*h > 1 a 

(5.48) R MW ( 7 ) := 7 - 2 7 3 - 1 | J (e T (x, x, 0) - /T 3 ^ A/" MW )^ 7 dx\ < 

ch- 2 + c^/r 1 + C/? 1 MW ( 7 ) 

uk*7i /?! MW (7) denned by (5.47). 

Remark 5.6.5. If we replace only e( 2 )(., .,.) by e MW (.,.,.) we can skip the 
last term. 



5.7 Geometric interpretation 

Note that there are two classical dynamics: in x' and in Xi. The former has 
return times t' k := t' k (r') where r' is a corresponding part of energy and t' k 
are defined in section 1 as tk(r) while the latter is 

(5.49) x 3 (t) = x 3 (0) + 2^(0) t - (3t 2 , &{t) = 6(0) - (it. 

and has return times (to 0) t k (r — r') = 2/3 _1 (r — r') and therefore the 
total system has return times 16 ) t k and return energy partitions (r' k ,r — r' k ) 
defined from the pair of equations 

(5.50) t k = t' k (r' k ) = t' k \r - t>) 

which is equivalent to (5.9)-(5.10). 

Therefore relatively sets of return directions is thinner in 3D than in 
2D which explains errors acquiring lesser factor (in comparison with the 
principal part) in 3D than in 2D. 



16 ) To point 0. 
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In particular, as f3 = ID-dynamics does not return at all unless it stays 
at 0, therefore in this case (5.49) mean exactly that r' k = r and t' k (r) = t k (r). 

6 Pointwise asymptotics: general 
3D-operators 

6.1 Set-up 

We assume that 

(6-1) F:=Q2gj k F j F k )^ 

jk 

where F = (F 1 , F 2 , F 3 ) and F are vector and scalar intensities of the magnetic 
field respectively. Temporarily we assume that F = 1 (as we can reduce the 
general case to this one by dividing operator by F and using our standard 
arguments). 

Further, without any loss of the generality we will assume locally that 

(6.2) F 1 = F 2 = 0. 

Remark 6.1.1. (i) Due to Frobenius theorem we can make locally g 3J = 
for j = 1, 2 and F 1 = F 2 = simultaneously if and only if A df* = where 
f * = Y.j,k F J gjkdx k is 1-form. 

(ii) As 

(6.3) F =(\Y1 S iJ S klF ikFji) 1 , F Jk = d k Aj - djA k , 

ij.k.l 

under assumption (6.2) 

(6.4) F=(g n g 22 -g' 2 g 21 ^\F 12 \. 

However 

(6.5) We can assume that F 1 = F 2 = locally and simultaneously that 
gi k = 8j k along a single magnetic line {x' = 0}. 
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Really, we can satisfy condition /* A df* = changing g J but preserving 
them at the chosen magnetic line. Note that F then acquires some scalar 
factor which does not affect this condition. So, we can assume g 3 -* = 
for j = 1,2 and F 1 = F 2 = at the chosen magnetic line which (after 
shift) becomes then {x' = 0}. Changing x 3 i->- 0(x 3 ) makes g 33 = 1 along 
this magnetic line. Then F 3 = ±F = 1 (for an appropriate orientation). 
Changing x' = £>(x 3 )x' with an appropriate matrix B(x') makes gi k = bjk 
along this line. 

Further 

(6.6) In frames of assumption (6.5) we can assume that g 31 = 0. 

Really, one can achieve it by x' i— >■ x' , x 3 i— >■ 0(x) with 0(x) = x 3 + 0(|x'| 2 ). 
Further, one can make 

(6.7) g jk = 5 jk + 0(\x'\ 2 ) asx 3 = j,k = 1,2 

by x' i-> x 1 + Q{x') with Q quadratic respect to x'. 

By a gauge transform one can make A 3 = 0. Then F 1 = F 2 = imply 
that Ai, A 2 do not depend on x 3 and therefore we can assume that Ai = 
as well: 

(6.8) A 1 = A 3 = 0, A 2 = A 2 {x')=x 2 + 0{\x'\ 3 ) 
where the last equality is due to F = 1 and assumption (6.5). 



6.2 Classical dynamics 

Let us consider a classical dynamics. We know that for time \ t\ < T = e 
Hamiltonian trajectory starting from {|x| < Cofi" 1 } remains confined in 
C 0/ u~ 1 -tube {|x'| < C /i -1 }. Then in virtue of our assumptions 

(6.9) ^ = 26 + oor 1 ) ^ = - V X3 (0, 0, x 3 ) + o^- 1 ) 

and therefore 



(6.10) 



(*3,6)(0 = (4 e 3 °)(0 + 0(^-^1) 
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where here (x° , £ 3 )(t) denotes dynamics for a lD-Schrodinger £| + V(0, 0, X3) 
with the same initial data (as t = 0). 
Then one can prove easily that 

(6.11) < C(\0\ ■ \t\ + 16(0)1), |*| < C{\(5\t 2 + 16(0)1 • \t\) 

if x 3 (0) = and \t\ < e 

with 

(6.12) ^ := _v x3 (0). 
Then 

(6.13) 6 = 6(0) (1 + 0(t 2 )) + 2/3t(l + 0(t 2 )), 

(6.14) x 3 = 6(0)t(l + 0(t 2 )) + /3t 2 (l + 0(t 2 )). 

Let us analyze evolution of (x', £')(t) but we start from (pi, P2)(0- Recall 
that 

(6.15) p-^Pi, P2} = 1, g n g 22 ~ g 12 g 21 = 1 mod 0(p~ 2 ) 

in the tube in question. 

Let p[ := pi + 7p 3 ; then p' 2 := {g 12 p x + g 22 p 2 + g 23 p 3 ) = -\^{a, p[} 
mod 0(// _1 ). Then 

\lT\a, p' 2 } = (g^pi + ^ 12 P2 + g 13 Ps)g 22 ~ (g 12 Pi + g 22 P2 + g 23 Pa)g 12 = 

pi + UV 2 -^ 23 ^. 

So we pick up 7 = g' 13 g' 22 — g 23 g 12 and then 

(6.16) p[ := (g 22 rk Pl + TPs), 7 = sV 2 - g"V 2 . 

(6.17) P^:=U 22 )^U 12 Pi + g 22 p 2 +^ 23 P3), 

satisfy 

(6.18) ^{a, pi} = -p 2 , ^{a, P2} = Pi. A* _1 {pi- P 2 } = 1 
modulo 0(p _1 ). We set p' 3 = p 3 . 
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We do not assume that g 13 = identically; instead its choice will be 
different and it would be only 0(/i~ 1 ) in the tube in question. Note that 

(6.19) l = g 13 g 22 mod 0(fi- 2 ) = g 13 mod 0(/i- 2 + /U - 1 |x 3 |), 

(6.20) l Xj =g 13 mod O^ 1 ). 

Really g* = b jk + 0(|x'| 2 + |x'| • |x 3 |) 0' = L 2), g 3J = 0(\x>\). 

Let us consider again {a, p[}, this time more precisely. Note that 

{a, p[} = -2/ip 2 + 2 ^ 7xj .p J p 3 - \/ Xl + 0(/i" 1 ) as x 3 = 

i 

as |x'| = 0(/i _1 ) due to our assumptions. 

Therefore we arrive to the first equation below and the second is proven 
in the same way 

(6.21) {a, pi} ee -2M + 2^4 3 PjP3- V*. 

(6.22) {a,p 2 } = 2M + 2 5> 2 x >,-p 3 - ^ mod O^" 1 ). 

j 

We want 

(6.23) g 23 = gl 3 = 0, -g 23 = 4 3 = «(*) as x' = 

which is possible to arrange as the only value at x' = which is fixed is 

(6.24) k{x 3 ) := X - (gl 3 - gl 3 ) | x , =0 = - 1 - (g 13 , X2 - g 23 , Xl ) | x/=0 ; 

one can prove easily this by transformation x 3 \-> x 3 + (?(x 3 ; x') with Q(x 3 ; x') 
quadratic with respect to x'. Note that 

(6.25) /* A </£ = 2/t(x 3 )c/xi A c/x 2 A c/x 3 as x = 0. 
Then (6.21)-(6.22) become 

(6.26) {a, p'J ee -2(/i + /tp 3 )p 2 - V X1 , 

(6.27) {a, p 2 } ee 2(/i + «p 3 )pi - V X2 mod O^" 1 ). 

Obviously 
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(6.28) We can replace in equations (6.26)-(6.27) V Xj by V Xj (x 3 ) : = V Xj .| x /= 
and p 3 by p%(t). 

If k = and V XJ did not depend on x 3 we would get the pilot-model 
operator. Unfortunately, it is not the case and (6.26), (6.27) coincide with 
those for the pilot-model only modulo 0( / u~ 1 + |p 3 | + |x 3 |). However note 
that «(x 3 )p 3 = ^K(x 3 ) where 

(6.29) K = ^JkcJx 3 . 
Therefore 

(6.30) p'(t) = fi(t)n _1 (o)p'(o) + f o.{t)Qr\e)V dt 

Jo 

mod 0(// _1 |t|) 

with 

P = (Pi- P2) f , V = (V'l/)', ft(t) fundamental matrix of system (6.26)-(6.27). 
Expressing Q _1 (t') via its derivative and integrating by parts we conclude 
that 

(6.32) p'(t) = ft(t)p'(0) +^fi(t)l/'(x 3 (0)) -pTV(x 3 (t)). 
Meanwhile, as 7 = 1,2 

(6.33) {a, (-l)Wi + Pi) = - X) &>^' " ^ mod 

k,i 

and we can replace in the right-hand expression x' by 0, p K by p' k and as 
/c = 1,2 we can replace p' k by the corresponding component of ft(t)p'(0). 
Then integrating by parts we arrive to 

(6.34) ((-l)J m _ j + Pj )\l=l = 

- / t GE(^ 11 +^ 2 )x,(Pl+P2 2 )+^(x3(0)) dt EE - f V Xj (x 3 (t)) dt 
JO j JO 

modulo 0( / u~ 1 ) as (g -11 + g 22 ) x = 0(p _1 ) in virtue of our assumptions 
including (6.15). Further, considering out-of integral terms we conclude that 
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(6.35) Equality (6.34) holds modulo (/x 1 ( 1 1 1 + | sin(2//0) + \p 3 \) (where 
p 3 (t) = p 3 (0) mod |t|). 

6.3 Semiclassical approximation 

So, after rescaling, the classical dynamics of the general operator is close to 
one of the generalized pilot-model 

(6.36) A := h 2 Dl + {hD 2 - /ixj 2 + h 2 D\ + V°{x), 

V°(x) := V(x 3 ) + ^j(x3)xj 

with 

(6.37) V(x 3 )= V(0,0,x 3 ), « i (x 3 ) = \/ x .(0,0,x 3 ), 7 = 1,2 

and the latter one is close albeit with a larger error to the dynamics of the 
pilot-model (5.1) with V(x 3 ), <Xj(x 3 ) replaced by V(0) + (3x 3 (f3 = V X3 (0)), 
cxj = otj(0) respectively. 

Proposition 6.3.1. (cf. proposition 2.2.1). Let conditions (2.1)-(2.3), 
(6.7) -(6.8) be fulfilled . Then after rescaling 

(i) Uniformly with respect to \ t\ < eii the propagator e lh 1{A is an h-Fourier 
integral operator corresponding to the Hamiltonian flow W t ; 

(ii) As |sin(26»i)| > e, \t\ > e 

(6.38) U{x,y,t) = 

M27rn)-ir5,Xsin(^ 1 ))- 1 e / ^ 1 ^ t )(^b m (x,y,0)^ m 

m 

with 6\ defined by 

(6.39) 0!(t) = 9(t) + K(x 3 (t)) - K(x 3 (0)) 

with 9(t), K(x 3 ) defined by (2.20), (6.29) respectively and with <fi defined by 
(2.25) -(2.29) and satisfying (with all derivatives) 

(6.40) (f) = 4>{6) + 0{fi- 1 t), 

(6.41) b m = 8 0m + O{fi- 1 t) 
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with 

(6.42) = Pn-Hte + y 3 ) + ^(x 3 - y^ft- 1 - p 2 fj,- 2 t 3 + 2 
with 4>2 defined by (1.10). 

Proof. Easy but tedious details we leave to the reader: comparing with 
the pilot-model or generalized pilot model we see that the structure of the 
canonical graph is close to one of the pilot-model. □ 

Proposition 6.3.2. (cf. proposition 2.2.2). Let conditions (2.1)-(2.3), 
(6.7) -(6.8) be fulfilled^. Then after resettling 

(i) Decomposition (6.39) remains valid as \t\ > e and |cos(#i)| < e; 

(ii) Further, after rescaling decomposition (6.39) remains valid as \ t\ > e 
and 

(6.43) C max(h, fi- 2 \t\) < |sin(0i)| < e 
albeit with an error not exceeding 

(6.44) Crt|t|-^|sin(^ 1 )|- 1 (V|sin(^ 1 )|) / 
and with <p, b m such that 

(6.45) \D p {<f> - 0)| < C /3/ u- 1 |t|«- 1 |sin(^ 1 )|-l /3 l V/3, 

(6.46) |D> m -6 m0 )| < C m ^- 1 |t||sin(^ 1 )|- m -l /3 l V/3, m. 

Proof. Again a tedious but an easy proof we leave to the reader. □ 

6.4 Tauberian estimates 

Remark 6.4-1. (i) Let us set x = y. Then we get expression (6.38) for 
U(y,y, t) albeit now due to definition Q\ = 9. 

(ii) Now we can continue in the same manner as before but to get a bet- 
ter estimate one needs to observe that the shift in x' (after rescaling) is 
now or 1 JV X /V • F _1 dt and its difference with ar x at does not exceed 
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C/i _1 |t| • | osc(V' V)\ where osc(l/l/) means an oscillation of 1/1/ along trajec- 
tory (until time t). 

However this trajectory must return to x 3 = y 3 and therefore one can see 
easily that along it p 3 = 0(/3/i _1 t) and x 3 — y 3 = 0((3ur 2 t 2 ) and we need to 
assume that this is less than ea i.e. 

(6.47) |t| < T, = e^min(l,ai/H), a := \W{y)\, (3 := \d y3 V{y)\. 

Therefore only as \ t\ < T* one can use nondegeneracy condition but it 
makes sense to consider t : 7"* < \t\ < T* = e/i. Then we have the following 
counterpart of proposition 5.2.1 where we need to sum (5.15) to k = and 
compare k = a 2 ha' 1 with T: namely, k < iff 

(6.48) a > max(fih,0^(fih)^,fi~ 1 ). 

Setting here /3 = 1 (as without spatial averaging or micro- averaging d\\fV/F 
is our foe and we need only estimate it from above) and thus 7~* = e^iah we 
arrive to estimate (6.50) below. 

On the other hand, if we do not use non-degeneracy condition we need 
to sum 

(6.49) Qj/ri|fr|"5 

3 3 1 

which results in Ca^bT^ 7~2 and we arrive to estimate (6.51) below. 

Proposition 6.4.2. (cf. proposition 5.2.1). Let conditions (2.1)-(2.3), 
(6.7)-(6.8) ; (0.5)-(0.6) be fulfilled. Then 

(i) Assume that conditions (6.47) -(6.48) be fulfilled as well. Then after 
rescaling as T is given by (6.47) with (3 = 1 

(6.50) |F t _ s -i T xr(t)r x U| < Cu.h~ 2 + 

j /J/7 _1 aH(l + (log(a/(^/i)i))+) as a > , 

I 3 3 1 2 

[/i2/? _ 2 72 as a<(u,h)3; 

(ii) In the general case as \ t\ < T* 

(6.51) \F„h-1tXtWxV\ < C ^ 2 + C ^ h-1 Tl 
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Then picking-up T = 7"* and dividing the left-hand expression of (6.50) 
by 7"* or picking-up T = T* and dividing the left-hand expression of (6.51) 
by 7"* we arrive to estimates (6.52), (6.53) below. 

Corollary 6.4.3. (i) In frames of proposition 6.4-2(i) with a > (/z/j) i 
(which implies fi < hr 1 ) 

(6.52) R T := |e(0,0,r) - e T (0,0,r)| < 

C/T 2 a-5 + Cfih^a-^l + (log(a/(///i)i)) + ) ; 

(ii) In the general case 

(6.53) R T < C/T 2 + C>/H; 

(Hi) In particular, according to (6.52), as ji < hr 1 , a x 1 

(6.54) R T < OT 2 + C/J/T 1 (l + \\ogfih\) 
which is 0(h~ 2 ) as /i < (h\ \ogh\)~i; 

(iv) In particular, according to (6.53), R T = 0(h~ 2 ) as fi < ; 

(v) Estimate (6.52) is better iff br^ < fi < h^ 1 and 

(6.55) a> fi^h- 1 + {fih) 1 2\\og{fih)\. 



6.5 Weyl estimates 

The proof of the following statement is rather obvious: 

Proposition 6.5.1. (cf. proposition 5.3.1). Let conditions (2.1)-(2.3), 
(6.7)-(6.8), (6.47) -(6.48), (0.5)-(0.6) be fulfilled. Then 

(i) As fi 2 h < a < 1 estimate (5.18) holds and correction term h~ 3 N XtC on 
is delivered by r-term stationary phase approximation /? _3 A/" X|Corr ( r ) with an 
error not exceeding (5.19). 

(ii) As a < ji 2 h estimate (5.19) holds; 

(Hi) In particular without any nondegeneracy condition R^ = 0(fj,^h~^) and 
it is 0(h~ 2 ) as < h~3 ; 

(iv) On the other hand, = 0(/i~5/j~§ 4- ^l/? -1 ) as a x 1 and fi < h~^ 
and therefore R^ = 0(h~ 2 ) as a x 1 and ji < h~s. 
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Combining with corollary 6.4.3 we conclude that 

Theorem 6.5.2. Let conditions (2.1)-(2.3), (6.7)-(6.8), (6.48), (0.5)-(0.6) 
be fulfilled. Then 

(i) In the general case 

(6.56) |e(x,x,0)-/7- 3 A/^(x,x,0)| < C/T 2 + oJ/ri; 

(ii) Under non-degeneracy condition 

(6.57) \V ±F V/F\^1 
as jj, < br^- estimates 

(6.58) | e(x, x, 0) - h~ 3 N^{x, x, 0) | < Ch~ 2 + C^h' 1 

and 

(6.59) |e(x,x,0) - r 3 AC(x,x,0) - /r 3 A/" x , corr ( r) | < 

c/r 2 + Cu,h- 3 2(n 2 h) r+1 2 

hold. 

6.6 Successive approximations 

Let us try successive approximation method. As an approximation we 
successively try the generalized pilot-model operator A , the pilot-model 
operator A and the magnetic Weyl approximation. As for fi < Weyl 
approximation delivers 0(h~ 2 ) error we assume that 

(6.60) H>h-s. 

6.6.1 Generalized pilot-model approximation, fih < 1 

We claim that 

(6.61) Effectively if x = y(= 0) we can estimate \\A — A°\\ by 

(6.62) £(k) := C(fi- 2 \k\ + (u,h) 1 2u,- 1 \k\- 1 2). 
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Really, due to our assumptions g-' 1 = bji + 0(|x'| 2 ) + 0(\x'\ ■ |x 3 |) with 
|x'| = O^ 1 ), x 3 = Oifi-^k]) = 1,2), g 33 = 1 + 0(/i~ 2 ) (before 
rescaling) the only term which does not allow such estimate is gj 3 DjD 3 with 
J 1.2. 

However for a classical trajectory to return back to time t x k we need 
to have p 3 = 0( / u~ 1 /c); in our analysis we will need also to satisfy uncertainty 
principle: p 2 \k\ > p/?; so we estimate p 3 by Cfi~ 1 \k\ + C(u,h)i \k\~* and as 
g j3 = 0(|x'|) this leads to (6.62). 

Then the error in k-th winding does not exceed 

(6.63) C/J/7-i|/r|-i x (fT 2 \k\ + (nh)^- 1 ]^) x p^/T^/clx 

f (^) r+l as «-V 2 ^<^l<(^ 3 ^)^ 

/ q;|k| 

[l as \k\ < m\n(a^ 1 u 2 h,(a 3 h)^). 

Really, the successive approximation makes sense only as £(/c)|/c| < fih which 
in view of (6.60) is equivalent to \k\ < (/z 3 /?)2. 

Here as \k\ > oT^^h we apply r-term stationary phase approximation 
as well. Also recall that this interval originally was a _1 /i 2 /) < \k\ < efia^ 17 ) 

2 

but it is non-empty iff a > (/u/?)?. 

Then (/i 3 /j)2 < /io^ and further (p, 3 /j)2 > or 1 a 2 h iff a > {fih)~2. Only 
under this assumption we need to consider the first case in (6.63) but it is 
the most important it happens for a x 1. 

Let us break (6.63) into two expressions: one with the second factor 
u,~ 2 \k\ and another with the second factor (/x/?) 2 /x" 1 1 /c | 2 . 

In the general case (when we do not know if a > (/i/?) 2 ") summation of the 
first expression returns C/? -2 (/! 3 /7)4 and of the second expression returns a 
lesser Ch~ 2 (l + log(/i 3 /))). Contribution of k : \k\ > (/i 3 /)) 2 " to the Tauberian 

3 3 3 3 3 1 

expression is a sum of Cp^/j^ 1 A: | 2 , which gives us Cp2 h~2~(fi 3 h)~4 i.e. 
exactly the same answer. Therefore we arrive to 

(6.64) The Tauberian expressions for the original operator and its general- 
ized pilot-model approximation differ by no more than Ch~ 2 (fi 3 h)^ . 

17 ) Here the upper limit meant that |x3| < a along trajectory before k-th winging. 
Since we look at the trajectory returning at k-th winging, |x3| < (yu _1 |/(|) 2 before it and 
we need to satisfy (^ x \k\) 2 < a. 
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Assume now that a > (f^h)^. Again break (6.63) into two expressions 
the same way as before. Obviously summation with respect to k of the 
second expression returns Ch~ 2 {l + (log(/i 2 /?/o;)) + ) no matter what r is. 

As r = summation of the first expression with respect to k returns 
Cfih~2 a ~2. One can see easily that the contribution of k : \k\ > (/U 3 /?)5 to 

3 1 

the Tauberian expression does not exceed Cfih~2a~2. 
Therefore we arrive to 

(6.65) As a > (/i/j)? the Tauberian expressions for the original operator 
and its generalized pilot-model approximation differ by no more than C/T 2 + 
Cnh~2a~2 . 

As r > 2 summation of the first expression with respect to k returns its 
value ( x | A" | ) as k = max(l, u. 2 h/a): 

la rr\ r 3 - u-1 i (~ Y * aS 1 ^ t^h/a, 

(6.66) C[i2h a 2 x < K a J 

y 1 as 1 < ffh/a 

and as r = 1 we get C/ii/?~ 1 a~i(l + log(a:(/i/?) _ i). One can see easily that 
the contribution of k : \k\ > (ft 3 h)2 to the Tauberian expression does not 
exceed these expressions as r > 1. 
Therefore we arrive to 

(6.67) As a > (^h)^ the Tauberian expressions with the subtracted h~ 3 Af x corr ( r 
i.e. 

(6.68) e T (x, x, r) - h~ W x , corr(r) 

for the original operator and its generalized pilot-model approximation differ 
by no more than 

(6.69) C/T 2 + Ch- 2 {\og{fi 2 h/a)) + + 

(l + log(a:(/!/))~5)) as r = l,a>n 2 h, 

„ 3 , _3 (l + \og( a 3 h)) as r = l,a<u 2 h, 

Cfj,2h 1 a 2 J v o Vf ^ " 

(/rh/a) 1 " 1 as r > 2, a > fTh, 

1 as r > 2, a < fi 2 h. 
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Combining with corollary 6.4.3 we estimate expressions (6.70)-(6.72) 
below for different a; in particular we arrive to 

Theorem 6.6.1. Let conditions (2.1)-(2.3), (6.7)-(6.8), (6.48), (0.5)-(0.6) 
be fulfilled 6 ^. Let e°(., ., .) be a Schwartz kernel of the spectral projector for 
operator A defined by (6.36) -(6.37). Then as h~3 < /i < h^ 1 

(i) In the general case 

(6.70) |e(x,x,0) - e°(x,x,0)| < qJ/T?; 

(ii) Under non- degeneracy condition (6.57) estimates 

(6.71) |e(x,x,0) - e°(x,x,0)| < Ch~ 2 + C>/H 
and 

(6.72) |e(x,x,0) - e°(x,x,0) - h~ 3 (j\f XiCon{r) - A/" x ° corr(r) ) | < 

Ch- 2 (1 + (log ///>)+) + Cu\h-\l + | log/i/)|) + 

f (1 + \og(^h)) as r = 1, 
C/J/T 1 \ {^hy- 1 as r>2,fi<hK 

[l as r > 2,/i > /T5 

hold. 

Remark 6.6.2. Replacing a(x 3 )xj by a(0)xj we estimate the norm of the 
perturbation by Ci(^) = A* 1 ' (a^ 2 !^! 2 + A* -1 ) which is less than expression 
(6.62) and thus leads to the lesser error. On the other hand, this new 
approximation is already a direct sum A = A( 2 ) + B with the pilot model 
operator A( 2 ) and 

(6.73) B = h 2 D 2 + V°{x 3 ). 

Therefore all above estimate remain valid with A defined by (6.36)-(6.37) 
albeit with aj = V Xj (0, 0, 0). 
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6.6.2 Pilot-model approximation, fih < 1 

Consider now the pilot-model approximation A which differ from the ap- 
proximation A of remark 6.6.2 by 0(x|). Then 

(6.74) Effectively if x = y we can estimate \\A° — A\\ by 

(6.75) £(k) := C(u,-*k* + u,- 1 h\k\). 

Really, we need to estimate 0(xf ) and according to above arguments we 
estimate |x 3 | by C[^r 2 k 2 + (u.h)i 1 A" | 2 x /i _1 |/c|) which leads to (6.74). 
Then we have estimate of the error in /c-th winding not exceeding 

(6.76) C/J/7-i|/r|-i x (u~ A k* + u- x h\k\) x fi^h^l^x 

{(l j2h V+\ 1 , 1 — -12, 

( — 7—7 ) 2 as \k\ > a urn, 
y a\k\' 

1 as \k\ < a -1 // 2 /?. 

The perturbation factor is 0(1) iff |/c| < /x/?s and we need to compare 
it with a: -1 // 2 /?: we have cases a > u.hs when a^/j^h is less and a < u.hs 
when /i/75 is less. 

4 

Case ct > /x/75. Then we need to sum 

(a) Expression (6.76) with the last factor 1 from \k\ = 1 to \k\ = a~ 1 u 2 h 
which returns Cu-ih~z \k\* x (ii~ 4 /c 4 + // _1 /)|/c|) calculated as |/c| = 
a _1 u 2 /?: 

C/i 2 /7 Z Q; 2 -I- C/i2a 2 ; 

as a x 1 the first term does not exceed Cu%h~ x iff /1 < /)~* while for 
the second term it is always true; 

(b) Expression (6.76) with the last factor (u 2 h / a\k\) r+ 2 from \k\ = a^ 1 u 2 h 

1 1 5 1 1 1 

to k = fihs which returns what we got earlier plus value of C/x5 /7 — 2 1 /r 1 2 x 

(/x~ 4 /c 4 + /i _1 /j|/c|) x (/i 2 /j/a|/c|)'' + 2 calculated at k = u.hs (may be with 

the logarithmic term for exceptional values of r = 4, 1) i.e. 

as r = 1 we get Cu\h~^a~i] as a x 1 this is less than Ch~ 2 as \x < h~^; 
note that the second term always is 0(/i2/7~ 1 + h~ 2 ); 
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(c) Qj/j2|/c| 2(fi 2 h/a) r+ 2 as \k\ > Cjih^ which returns the expression 
above. 

In total, we get modulo 0(h~ 2 + /ii/? -1 ) 

(6.77) C^tfa-i + QJcH + C/i/rt^ /a) r+1 2. 

4 

Case a < /ihs. Then we need to sum 

(a) Expression (6.76) with the last factor 1 from \k\ = 1 to \k\ = phi which 
returns Cu.2 h~2\k\2 x (fi~ 4 k A + calculated at k = u,hs i.e. 

Cu.h~* + Cph~i; 

here the first term is larger than CjAhi^ + Ch~ 2 as ji > h~i but the 
second term is always smaller than this; 

(b) Ca%\n as \k\ > Cp\v= which returns an above expression. 
So, we get in this case 

(6.78) Qi/ri. 
Thus we arrive to 

Theorem 6.6.3. Let conditions (2.1)-(2.3), (6.7)-(6.8), (6.48), (0.5)-(0.6) 
be fulfilled 6 ^. Let e(., ., .) be a Schwartz kernel of the spectral projector for 
the pilot-model operator A. Then as h~3 < ji < h^ 1 

(i) In the general case 

(6.79) |e(x,x,0) - e(x,x,0)| < C/J/T* + Cph'l; 

4 

(ii) Under non- degeneracy condition (6.57) as jj, < h~5 estimates 

(6.80) |e(x,x,0) - e(x,x,0)| < Ch' 2 + C '///H + Cfi^ h 2 + C pS + C /J h~l 
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and 

(6.81) |e(x,x,0) - e(x,x,0) - h~ 3 (Af x , corr{r) - A/J corr(r) ) | < 

Ch- 2 (1 + (log/i 2 /))+) + Calh-\l + | \ognh\) + 
( (1 + log^ 3 /))) as r = 1, 
C/J/T 1 i {/Shy- 1 as r > 2, fi < + 

[ 1 as r > 2, u. > h'^ 

C^h 2 + C / u5C/i/)~i(/i/?^) r+ 3. 

Let us improve the above results. Note that we need only consider 
components of the estimate which are (in unrescaled x 3 , p 3 ) due to p 3 = 
0(fi^ 1 k), x 3 = 0(fi~ 2 k 2 ) rather than those which are due to the uncertainty 
principle p 3 = 0((u,h)2\k\~2), x 3 = 0(fi~ 1 (fih)2 1 /c 1 2 ) as the latter brought a 
proper estimate in the above analysis. 

Consider only case a x 1, f3 < 1. Then we need to consider only p > /7~i. 
Note that actually instead of p 3 = O(t), x 3 = 0(t 2 ) we can use estimates 
p 3 = 0((3t), x 3 = 0(/3t 2 ) and then our estimates acquire factor (3 2 (recall 
that we are not discussing terms which are due to the uncertainty principle). 

Case /r 16 / 25 < u, < /r 4 / 5 . Consider (6.77) with r = 1; then the last term 
is the largest and the estimate with above improvement is 0((3 2 /j, 5 ^ 2 h~ 2 ^ 5 ) 
and it is OOu 3 / 2 /)- 1 ) as (3 < {^)—2. So we need to consider case 

(6.82) (3 > (M/7§)-i (> h 1 ' 10 ). 

As r = we get pJ/7 _ i which is less than in (6.71); 

Case h^ 5 < fi< ft" 1 . Then estimate is (modulo O^/ 2 ^ 1 )) O^V/r 8 / 5 ) 
and it is OOu 3 / 2 /)- 1 ) as /3 < (///»§)*. So we need to consider case 

(6.83) (3 > {u,hl)* (> ft 1 / 10 ). 
Therefore we conclude that 
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Theorem 6.6.4. As p := |V|| F V/F| < h' 1 /™ estimates (6.70) -(6.72) of 
theorem 6.6.1 remain true for A replaced by A. 

So we need to consider the case when lD-Schrodinger operator is almost 
non-degenerate and the corresponding dynamics is confined to zone x 3 < (3; 
rescaling x 3 h-> x 3new = x 3 //3 we arrive to operator 

(6.84) /3 2 (/> 2 D 3 2 +\/ 3 (x 3 )), /7 3 = V/? 2 «1- V 3 (x 3 ) = P' 2 V(x 3 ) 

satisfies d 3 V 3 (x 3 ) x 1. Then we can construct its solution by WKB method. 
Note that before any rescaling 

(6.85) U w {x 3 ,y 3 ,t) = 

(27T/))- 1 J exp (V 1 (S(x 3 , t, C) - y 3 C)) B(x 3 , t, () d( 

where 5 solves 

(6.86) S t = S 2 + V(x), S| t=0 = x 3 C 

and B ~ B^h k in the standard way. 
Then one can prove easily that 

(6.87) S = S + o(x 2 |t| + |x 3 |(|C| + \t\)\t\ + ( 2 \t\ 3 + |t| 5 ) 
where 

(6.88) S = x 3 ( + t(C 2 - /3x 3 ) + C/?t 2 + ^/3 2 t 3 

solves the same problem for pilot-model V = — /3x 3 . Rescaling x 3 — > x 3 /(3, 
h i — y h/(3 2 shows one can calculate 5 and 5 for (3 = 0, then plug x 3 = x 3 /(3, 
C = (/P and multiply the result by /3 2 ; (6.87) obviously survives but (6.88) 
improves to 

(6.87)* S = S + o(x 3 2 |t| + |x 3 |(|C| + P\t\)\t\ + C 2 |t| 3 + |t| 5 ) . 

Obviously phase function 

(6.89) 0(x 3 ,y 3 ,t,C) :=5(x 3 ,t,C)-y 3 C 



Chapter 6. Pointwise asymptotics: general 3D-operators 



106 



is equivalent to phase function 

(6.90) := (C + (3t)x 3 - (C - (3t)y 3 - 2( 2 t - ^ 2 t 3 

we used earlier. 

Furthermore, obviously 

(6.91) Phase function 

(6.92) 0(x 3 ,y 3 ,t,C) :=S(x 3 ,t,C)-y 3 C 

is ID-action and therefore the total exponent is also 3D-action associated 
with the generalized pilot-model. 

Note that for a phase function (6.89) C = iff x 3 - y 3 + 2t( + (3t 2 = 
and then = — At and therefore for \t\ >> h we can apply a stationary 
phase method for both S and S and while the principal term would be of 
magnitude (/?|t|)~2 the error in /-term approximation will be of magnitude 
(/?|t|) _ 2(/)/|t|) / or after all the substitution to U and rescalings it will be 

(6.93) C/J/7-i|/c|-i x {^h/\k\)' 

in the Tauberian expression. While almost useless for fih close to 1 and 
k = 1 the last factor is very important for larger k. 

Then the contribution of zone {k : \k\ > m} to such error does not 
exceed (as / = 1) 

(6.94) Cfih-^m-l 

and it is 0( / u _ i/7~ 1 ) as m > jAh^. Obviously m <^ fih^ and m <^ ji 2 h which 
means that in zone k : \k\ < m we could perfectly deal with successive 
approximations. 

On the other hand, replacing B by 1 and <fi^ by —At leads to an error 

(6.95) C/J/7-i|/c|-i x fi~ 2 k 2 x min((/i 2 /)/|/<|)^ l) 

which after summation becomes /is/? -1 x ji which is not as good as the 
Tauberian estimate. This spares us from more complicated formula. We 
leave to the reader 

Problem 6.6.5. (i) Write the approximation due to described combination 
of successive approximations and WKB method. 

(ii) Consider a < 1. 
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6.6.3 Magnetic Weyl approximation 

In virtue of the above results the result of subsubsection 5.6.1 hold for 
general operators. 

6.7 Strong and super-strong magnetic field 

The standard method of successive approximations shows that replacing 
operator by the generalized pilot-model brings an error with the contribution 
of {t : \ t\ x T} (before rescaling) not exceeding 

(6.96) C>/T 2 x (n^h + Cfi-h^T 1 ! + h 6 )) 

and summation with respect to T running from h to 1 returns 

C>/T 2 x (u.~ 1 h+ Cu.-h-*(l + h s \\ogh\)) x qJ/t§ + log//| 
which is less than the Tauberian estimate. Therefore 

Theorem 6.7.1. Let conditions (2.1)-(2.3), (6.7)-(6.8), (6.48), (0.5)-(0.6) 
be fulfilled 6 ^ . Then as pih > 1 

(6.97) |e(x,x,0) - e°(x,x,0)| < C>/H. 

Similarly, transition to the pilot model brings an error not exceeding 

(6.98) Cfih- 2 x (fi 2 T 4 + hi) 

and summation with respect from h to T* returns C/i/?~i + Cjih~ 2 f3 2 T^. 
On the other hand, Tauberian approach shows that the contribution of 
|t| > T* should not exceed C/i/?~3 7"" 1 and minimizing it sum we arrive to 

2 1 

7~* = /3 - 5/)io and it is 

(6.99) C/x/ri(l + /3§/T*>). 
Therefore 

Theorem 6.7.2. Let conditions (2.1)-(2.3), (6.7)-(6.8), (6.48), (0.5)-(0.6) 
he fulfilled . Then as fih > 1 

(6.100) |e(x,x,0) - e(x,x,0)| < C///ri. 
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Further, if we pick up T* = min(/3 2 hi, 1) we conclude that the contri- 

3 

bution of zone {t : T* < \t\ < e} to an error does not exceed Cfih~?. So 
we need to consider zone {t : \t\ > h«}. However in this zone we can apply 
stationary phase expression for L/( 31 ) which gives an error 

c |V/,- 2 x(^^xC^- 2 (A)i 

3 2 

which is 0(/i/j~2) as T > h^. 

Problem 6.7.3. Construct this stationary phase approximation. 

Again for the magnetic Weyl approximation error we refer to subsubsec- 
tion 5.2. 

6.8 Micro- averaging 

We are interested only in Tauberian, Weyl and magnetic Weyl estimates as 
we consider pilot-model and generalized pilot-model approximations a bit 
too complicated. We consider here only isotropic micro- averaging (i.e. with 
73 = 7)- 

6.8.1 Tauberian estimates 

Basically everything remains as in section 5.4 but we need to answer a 
question: "what is £?" Or rather "what T* we need to take?" Clearly 
T* = t min(o;, f3) fits but can we do better than this? 

(i) Let (3 < a. Then as we take 73 = 7 we have 2/(7) x 07 so averaging 
with respect to x 3 does not matter and therefore we do not need to assume 
that 73 < e(3 (which in this case would mean 7 < e(3). In this case our 
assumption is 7 < ea. 

Then as V±fV/F has a less oscillation then ea with a = \V±fV/F\\ 
for time 7"* = e^a we can take 7"* = e^cx. 

(ii) Let (5 > a. Then as we take 73 = 7 we have u(-y) x ^7 so averaging 
with respect to x' does not matter as well as the shift and therefore we do 
not need to assume that 7 < ea (which in this case would mean 73 < e(3). 
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Then as V\\fV/F has a less oscillation then e(3 with (3 = |V||f V/F\ \ t=Q 
for time 7"* = ei/3 we can take T* = e\f3. However in this case we should not 
include factor (fj,h/a\k\)2 which is a loss as \ihja < h//3j i.e. as /xj < a(3. 

So, in both cases we can take 73 = 7 and 7"* x |V V/F\ (but we obviously 
need 7 < |W/F|). 

6.8.2 Weyl and magnetic Weyl estimates 

As we do not care about 7"* anymore we just take 73 = 7 < |W/F|. 

Problem 6.8.1. As (5 is rather an obstacle in the pilot-model approximation, 
it would be interesting investigate this approximation with micro- averaging. 



7 Dirac energy: 3D-estimates 

In this and the next sections we consider asymptotics of I defined by (0.2). 

7.1 Tauberian formula 

Let us consider first contribution of zone {(x, y) : |x — y| > C7}. 

Proposition 7.1.1. (cf. proposition 3.1.1). Let fi < hr 1 . Then under 
conditions (0.5)-(0.7) the contribution of zone {(x,y) : |x — y| > C^}, 
to the remainder is 0(/) _2 7 _K ) while the main part is given by the same 
expression (0.1) with e(x,y,0) replaced by its standard implicit Tauberian 
approximation with T x e (0.11). 

Proof. Proof repeats one of proposition 3.1.1 albeit according to Chapter 13 
estimate (3.4) is replaced by 

(7.1) || E{t,t')$ ||i < Oi- 3 (\t-t'\ + CT^h) 

Vr, t' e [-e, e], T = efi 



with T = e. 



□ 
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Proposition 7.1.2. (cf. proposition 3.1.3). Let conditions (0.5)-(0.7) be 
fulfilled. Then 

(i) As < k < 3 and either k ^ 1,2 or k = 1,2 and u(x, y) is replaced by 
u(x,y) — x(|(x + y))|x — y|~ K with an appropriate smooth coefficient x(x), 
with the error 0(h~ 2 ~ K ) one can replace e(x,y,r) by its standard Tauberian 
expression (0.11) in the formula (0.11) for I. 

(ii) As k = 1,2 and u = x(|(x + y))|x — y|~ K , with the error 0(h~ 2 ~ K \ log h\) 
one can replace e(x,y,r) by its standard Tauberian expression (0.11) in the 
formula (0.11) for I. 

Proposition 7.1.3. (cf. proposition 3.1.6). Let conditions (0.5) and (0.6) 
be fulfilled. 

(i) Further, let either condition (0.8) 18 - ) be fulfilled or < h 5 ~ 1 \ log | 1 
with an arbitrarily small exponent. Then (0.11) and statements (i), (ii) of 
proposition 1.1.2 hold. 

(ii) Furthermore, let h^ 1 ^ 6 ^ < br 1 . Then (0.11) and statements (i), (ii) 
of proposition 7.1.2 hold with an extra factor (1 + /x/? 1-5 ) in the right-hand 
expressions. 

7.2 Superstrong magnetic field 

Consider Schrodinger-Pauli operator as fih > eo which is a bit more tricky: 

Proposition 7.2.1. (cf. proposition 3.2.1). Let fih > eo- Then 

(i) Contribution of zone {(x,y) : |x — y| > 7} to I does not exceed Cfih~ 2/ -f~ K ; 

(ii) Further, 

{/T K as < K < 1, 

/)- 1 (l + (log/i/?) + ) as k = 1, 

^{K-l) h -\{n+l) as 1<K<3 

(Hi) Furthermore, I = 0(fih 00 ) under condition (3.14); in particular it is the 
case as 3 < 1 and \xh> Cq. 

18 ) Which actually could be replaced by much weaker non-degeneracy condition of 
Chapter 13 for d — 3. 
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Proof. Statement (i) trivially follows from the fact that Jzf 2 -norm of e(.,.,r) 
does not exceed C^h" 1 . 

Meanwhile estimate of the contribution of the zone 

{(x,y) : x G 8(0, 1), ye 8(0,1), |x - y\ < 7} 

is more subtle: if < /« < 1 u = d X3 u' with 00' = 0(7 1_K ) as |x — y| < 7. 
Then 

/ 7 = J \e(x,y,r)\ 2 uj(x,y)^(x - y) dxdy = 

— 2h~ 1 Re / J hD X3 e(x, y, r) • e(y, x, r)cu'(x, y)ip!y(x — y) dxdy— 

J \ e(x,y,r)\ 2 uj(x,y,T) ■ d^vb'^x - y) dxdy 

with the first and second terms not exceeding C , -f 1 ~ K h~ 1 x fih~ 2 and C7~ K x 
/i/? -2 respectively where we used the fact that Jzf 2 -norm of P/e(., ., r) does 
not exceed Cfr^h' 1 . Setting 7 = /? finishes the proof of (ii) in this case. 

As 1 < k < 3 we consider a partition of {(x,y) : |x — y| < h} to subzones 
{(x, y) : |x' — y'| > a} and {(x, y) : |x' — y'| < a} and in the former zone 
we repeat the same arguments as before albeit with u' = O(log \x' — y'\/h), 
u)' = 0(\x' — y / | 1_K ) for k = 1 and 1 < k < 3 respectively; then contribution 
of this zone does not exceed C/7 _1 | \oga/h\ x u,h~ 2 and Ch~ 1 a 1 ~ K x /i/7~ 2 
respectively; as a = or 2^ we arrive to estimate (7.2) for contribution of 
this zone as well. 

Finally, contribution of the zone {(x,y) : |x' — y'\ < a, |x 3 — y 3 < 7} does 
not exceed Cfi 2 h~ A x J |z|~ K dz with the integral taken over cylinder {\z'\ < 
°~, l z 3 < 7} and not exceeding C(l + | loga/7|)cr 2 and Ca 2 ~ K respectively 
as |e(x,y, r) = 0(fih~ 2 ). Again for selected 7, a we get estimate (7.2) for 
contribution of this zone as well. 

Finally, statement (iii) is due to the fact that e(x,y,0) = 0({ih°°) as 
condition (3.14) is fulfilled in one of points x,y; see subsection 13.5.4. □ 

Proposition 7.2.2. (cf. proposition 3.2.2). Let conditions (0.1)-(0.5) be 
fulfilled. Let fih > e and one of the nondegeneracy conditions (3.15) or 
(3.17J be fulfilled. 
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(i) As a > Co/j, 2/72 contribution of the zone 

{(x,y): \x -y'\ > C fi~h^, |x — y| > 7} 

to I does not exceed Cii\r x ^ K 7 

(%%) As 7 > C /? contribution of the zone 

{(x,y) : |x' -y'| < C /x~M, |x 3 - y 3 | > 7} 

to I does not exceed C/i/) _1 7~ 1 ~ K ; 
(Hi) Estimate 

(7.3) |I-I T |< 



C///T 1 < 



' h~ K os < k < 1, 

+ (log /!/))+) as K = 1, 

k uhh~^ (l + (log/i/?) + ) as k = 2 



holds. 



Proof. To prove statement (i) recall that drift speed does not exceed C//" 1 
and therefore Hilbert- Schmidt norm of vbE(r)vb' does not exceed C/i5/7 _ 2 as 
■0, '0' are -^o 00 "^ 110 ^ 0118 w ith dist(supp yb, supp?//) > C /i~2/)2. 

Really, it is true for a Hilbert-Schmidt norm of (E(t) — E(r / ))^ with 
|r — t'\ < h and then by Tauberian theorem it is true for a Hilbert-Schmidt 
norm of (E(t) — E T (r))^' with E T operator with the Schwartz kernel e T 
with time T x 1: 

(7.4) \\(E(r)-E T (r))vb'\\ HS < C^h'l. 

However ipE J ip' is negligible as T < e due to propagation results. This 
implies (i) obviously. Statement (ii) is proven in the same way but now 

rx 7 . 

To prove statement (iii) we need to estimate contribution to I — l T of the 
zones {(x, y) : \x' — y'\ < h, h < |x 3 — y 3 | < e} and {(x, y) : |x — y| < /?}. 
However estimate of the contribution of the former is trivial due to estimate 
(7.4). So, consider zone {(x,y) : |x — y| < h}. 
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(a) As < k < 1 we can replace as in the proof of proposition 7.2.1 u(x, y) by 
cu'(x, y) = 0(/j 1 ~ k ) in this zone and one copy of (e(x, y, r) — e T (x, y, r)) to 
D 3 (e(x,y, r) — e T (x,y,r)); however subsection 13.5.4 implies the following 
generalization of (7.4): 

(7.5) \\P*P' j (E{t) - E T (r))^|| HS < C^h^^hy 

with P' = (Pi, P2) and therefore contribution of this zone to (I — l T ) does 
not exceed Cfih~ 1 ~ K . In this case estimate (7.3) is proven. 

(b) Let 1 < n < 3. Then the above arguments imply 

(7.6) As 1 < k < 3 contribution of the zone {(x,y) : \x' — y'\ > a} with 
a > Co/i~2/?2 to I does not exceed Cfih^a 1 ^^; as k = 1 it does not exceed 
Cfih- 2 {1 + |log/i/?|). 

Note that with a = Cofj,~ih$ we get exactly the right-hand expression 
in (7.3). Therefore we need to estimate contribution of the zone {(x,y) : 
|x' — y' < a, |x 3 — y 3 | < h} with o = C fi~2hi. 

Consider case 1 < k < 2 first. Note that u = (3 — k)~ 1 J2j d Zj ZjUJ 19 -* 
and therefore we can replace u either by (xj — yj)u) = 0(\x' — y'\ ■ \x — y|~ K ) 
with j = 1,2 or by (x 3 — y 3 )w = 0(|x — y| 1_K ), simultaneously applying 
h~ x Pj to one of the factors (e(x, y, r) — e T (x, y, r)) . Then we can replace 
(xj — yj)uj by (xj — yj)u' simultaneously applying h~ x Pj to one of the factors 
(e(x, y, r) — e T (x, y, r)) (the same as before or another one). 

So as J = 1,2 we estimate the term by Cfih^ 1 x h~ 2 (u,h)^ x a x 
with factors (/i/?) 5 and cr coming from Pj and (x,- — y); as k = 1 the last 
factor is replaced by | log cr/'y|- As 7 = 3 these factors are replaced by 1 
and h respectively with the same product. Note that when we drag the 
last d X3 through x 3 they may cancel one another but the power becomes 
(1 — k) which we treated already (and while the gain in this power is exactly 
compensated by an extra factor h^ 1 .) 

Plugging a we get a proper estimate (7.3). 

(c) The same arguments work for 2 < k < 3 as well albeit now u = 
(3 — K,)^ 2 J2jk^Zjd Zk ZjZkU. Note that when we drag the last d X3 through 
x 3 they may cancel one another but the power becomes (1 — re) which we 
treated already (and while the gain in this power is exactly compensated by 

19 - 1 Recall that we replace (x, y) by (x, z) — (|(x + y), x — y). 
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an extra factor h 1 .) However as k = 2 we get power —1 and then logarithm 
comes in. □ 

8 Dirac energy: 3D-calculations 

8.1 Pilot Model 

8.1.1 Transformations 

Consider first the pilot-model operator (5.1). 

Let us rescale as before 20 ^. Then U(x, y, t) in comparison with U( 2 ){x' , y', t) 
has factor l/(i)(x 3 ,y 3 , t) defined by (5.4). Let us consider effect of all these 
changes ignoring other variables; again without any loss of the generality 
we assume that u(x,y) is given by (4.3); 

We can replace variables x,y with new variables x := |(x + y) and 
z := (x — y) and rescale. Note that the phase of U(i) is linear with respect 
to x 3 and we can get rid off [i~ 1 dx 3 integration in the same manner we 
got rid off dxi. we replace Q(... ,x 3 , .) by its partial Fourier transform 
x 3 — > — /i~ 1 /)~ 1 (t / — t") with an extra factor 2ir; thus in comparison with 2D 
case we have 

- Q(2 / u- 1 /)- 1 a(t / - t"), -fi-^ptf - t"), z) 

where Q(., ., z) is a partial Fourier transform with respect to (x 1: x 3 ) (instead 
of Q(2 /U - 1 /)- 1 a(t'-t // ),z)); 

- an extra factor which after easy reductions becomes 

-/)- 1 (t / t // )-^expf/r 1 (t / -t // )(-J(t / t // )~ 1 z 3 2 -^/3V 2 (t /2 +t / t // +t //2 ) N ), 
4 V 8 3 / 

- integration with respect to dz 3 which enters both this factor and Q(., ., z). 

2 °) We need to add factor [i~ 6+K — /i -3 x x /i~ K where two factors fi~ 3 are coming 
from dx and dy and [i K from oj. 



Chapter 8. Dirac energy: 3D -calculations 



115 



Taking in account all these modifications we need to replace (4.4) by 
(8.1) 1 = 

^(27r)- 3 /rV +3 J J JJJ ^(Ah^as, -/T^s, z)x 



8 

csc(t + s) csc(t — s)(t + s)~i(t — s)~? x 



3 



exp(/7i 1 --(cot(t + s) -cot(t-s))(z 2 + (z 2 + 2t/i x a) 2 + 4s 2 /i 2 ct 2 ) 
- (cot(t + s) + cot(t - s))s/i _1 a - 2s(r - /i _2 a 2 )+ 
2s(-^(t + s)- 1 (t-s)- 1 z 2 - ^V 2 (3t 2 + s 2 )]) dtds^{ f i' 1 z)dz 1 dz 2 dz 3 . 

8.1.2 Case A 2 > ///i 

Assume first that 

(8.2) X 2 := a 2 + P 2 > h 1 - 5 . 

Remark 8.1.1. (cf. remark 4.1.2). (i) We take A-admissible with respect to 
x function Q\ = Q(x/A, .). 

Then the virtue of the factor Q.{Ah~ l a r ys, —h~ 1 (3 / ys, •) under assumption 
(8.2) we need to consider only |s| < h s and therefore we can consider 
separately \t'\ < e , \ t"\ < eo and \ t'\ > eo, \ t"\ > eo; 

(ii) Note that due to section 6.3 contribution of zone {\t'\ < eo, \ t"\ < eo} 
defined by integral expressions (8.1) etc with an extra factor x eo (t') or Xe (t") 
or Xe {t) differs from the same expression for non-magnetic Schrodinger 
operator by 0(fih~ 2 ~ K x h K ) = 0(u. K+1 h~ 2 ) as k ^ 1 and by 0(n 2 h~ 2 \ log/i|) 
as K— 1. 

(iii) Furthermore, if we remove from this expression for a non-magnetic 
Schrodinger operator cut-off {|r'| > eo} then the error would not exceed the 
same expression as well. 

Let us consider contribution of zone {\t'\ > e , \ t"\ > e } defined by an 
integral expressions (8.1) etc with an extra factor (l — Xe Q {f)). Due to 
remark 8.1.1(i) we need to consider only t' , t" belonging to the same tick. 
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Let us consider first zone 

(8.3) {\s\ > h/X 2 , |sin(t)| > C\s\}. 

Then integration by parts with respect to z' delivers one of the factors (4.7), 
(4.8). Thus integrating by parts many times in the zone where both of these 
factors are less than 1 we acquire factors (4.9), (4.10) respectively. 

On the other hand, integration by parts with respect to Z3 delivers 
factor ft|s| _1 t 2 |z3| 2 ; therefore integrating by parts many times in the zone 
{| z 3| > h,2\s\ ~2|t|} delivers factor 

(8.4) (l + lsl/r^lzsl 2 )-'. 

Case 1 < k < 3. Multiplying by |z| _K and integrating we get after 
multiplication by |sin(t)|~ 2 expression (4.11) albeit with k replaced by 
(k — 1) i.e. 

(8.5) x |sin(t)r~ K |sri + Hi-K 

Integrating by t over k-th tick intersected with {t : | sin(t) | > |s|} we get 
(4.12) modified the same way i.e. 

{|s|-5 + 5 K 1<K<2, 
|s|-5(l + |log|s||) K = 2, 

|S|3~5 K 2 < K < 3. 

This expression (8.6) must be either integrated with respect to s over 
{|s| < h/\ 2 } or multiplied by (h/X 2 y\s\~' due to factor Q and integrated 
over {|s| > h/X 2 }, resulting in both cases in the same answer which is 
the value of (8.6) x |s| calculated as s = h/X 2 which is similarly modified 
expression (4.13) i.e. 

r A 1 "" 1< k < 2, 

(8.7) Chl(l + I log(n/A 2 )|)A- 1 k = 2, 

[ h 2 - K X- 3+K 2<k<3. 

In addition to zone (8.3) we need to consider zone (4.14) defined as 
{|sin(t')| x |s|, |sin(t")| < |s|}; its tween {|sin(t')| x |s|, |sin(t")| < |s|} is 
considered in the same way. 
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In zone (4.14) | cot(t') — cot(t")| x | sin(t")| 1 and in this case factors 

(4.7) , (4.8) are replaced by (4.7)', (4.8)' and (4.9), (4.10) by 4.9, (4.10)' while 
factor (8.4) is preserved. 

Then, multiplying by |z|~ K and integrating we get after multiplication 
by | s i n ( t /r ) | 1 1 sin(t / )| 1 similarly modified expression (4.11)' 

(8.5) ' ni-3«| s in(t // )|5-5«| s |- 1 ; 

then integrating by \t"\ over one tick but with | sin(t")| < |s| we get modified 
expression (4.12)' i.e. 

(8.6) ' h^ K \s\^ K . 

Finally, either integrating over {|s| < h/e\} or multiplying by |s| _ '(/j/£A)' 
and integrating over {|s| > h/e\} we get in both cases the same answer 
(8.6)' x |s|, calculated as s = h/\ 2 not exceeding (8.7). 

Therefore the total contribution of zones (8.3) and (4.14) is given by 
expression (8.7). Then multiplying by |/c|~ 3 /i K /7~ 3 A 3 we get after summation 
with respect to k : \k\ > 1 the value as k = 1 i.e. (8.7) x fi K h~ 3 \ 3 . 

Therefore we arrive to analogue of proposition 4.1.3 

(8.8) \\ l3 mbda T -\ T X '\ < /? W (A) := 

r A 1 "" 1< k < 2, 

Qi K+1 /T 2 A 3 i (1 + | log(n/A 2 )|)A- 1 k = 2, 

{ fi 2 - K h 2 - K \- 3+K 2<k<3 

where I J' is a Tauberian expression for I albeit with T = eo/i -1 . 

Case < k < 1. In this case integration over {z : |z| < C } and 
multiplication by | sin( t) | 2 results in 



(8.9) C|sin(t)|- 2 p 2 



C 1_K < k < 1, 

(l + (logC/p))+) « = 1 



with 

(8.10) C = C(s, /c) = min((^/|s|)i|/c|, l) , 

(8.11) p = min((£/|s|)U) 



21 ) Corresponding to |x — y| < Cq/i 1 . 
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(cf. expression (8.5)). 

One can see easily that 

(8.12) Integral by t over one tick intersected with {t : | sin(t) | > |s|} does 
not exceed C/i|s|~ 1 C 1 ~ K as \s\ > h~ and C(/i|s|~ 1 )2^ 1_K as |s| < h~. As k = 1 
one needs to include an extra factor | \og(/p\, p calculated as | sin(t)| = 1 
for \s\ > h and sin( t) | = (\s\/h)2 for |s| < h~. 

Meanwhile in the zone (4.14) integration over {z : |z| < Co} and mul- 
tiplication by s x | sin(t // )| 1 results in expression (8.11) and integration 
by {t" : |sin(t // )| < |s|} results in much lesser term Ch( 1 ~ K , with some 
logarithmic factor as k = 1. 

Now integral by {s : |s| < 1} does not exceed Ch\k\ 1 ~ K as < k < 1 and 
C/i(l + | log h\ + | log \k\) as k = 1. Finally, multiplication by /i K+3 /i~ 3 |/c|~ 2 
and summation with respect to k : \k\ > 1 results in Cp K+1 h~ 2 as < /t < 1 
and Qi 2 /?~ 2 | logyu/?| as /c = 1. 

Therefore we arrive to analogue of proposition 4.1.3: 



(8.13) \\iambda T 



Combining with (8.8) we arrive to 




1 < K < 1, 

(l+|l0g(/i/))|) « = 1. 



Proposition 8.1.2. (cf. corollary 4-1-3). (i) For the pilot-model operator 
with 1 < fj, < h^ 1 under assumption (8.2) estimates (8.8), (8.13) hold for 
l</t<3 ; 0</t<l respectively. 

(ii) In particular, as A = 1 

(8.14) |I T - l T, | < R w := 

{1 < K < 2,K^ 1, 

(l+|l0g(/i/))|) K=l,2, 
l/- K h 2 ~ K 2<k<3. 

Following the same proof as in proposition 4.1.4 one can prove easily 
that 

(8.15) |I T, -X W | < Cp 2 h- 1 - K 
which does not exceed /? W (A) implying 
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Corollary 8.1.3. (cf. corollary 4- 1-5). In frames of proposition 8.1.2 

(8.16) \\~[-Z™\ < C/? W (A). 

Remark 8.1.4- After summation with respect to partition we restore esti- 
mate (8.14) under some non-degeneracy assumption. Comparing with the 
Tauberian remainder estimate 0(h~ 2 ~ K ) we conclude that /? w is lesser as 
pi < /?~ K /( K+1 ) (0 < k < 2) and /i < h~3 (2 < k < 3) (with some correction 
for k= 1,2). 

In particular, /j, < br^ would suffice for k > 1 and therefore improvement 
as in subsubsection 4.1.3. "Improvement" is not needed. 

8.1.3 Case A 2 > /ih 

Assume now that /i -1 < A and A 2 < fih 1+s . Then in contrast to the previous 
we will need to compute contributions of pair of ticks (k', k") with k! ^ k". 
Then if r', t" belong to k'-th and k"-th ticks respectively we denote r = k' — k" 
and s=t'-t"- 2vrr. 

However we still need to distinguish cases 1 < k < 3 and < k < 1. 

Case 1 < k < 3. In this case we according to the previous subsection need 
just to take results for d = 2, plug k := (k — 1), a = A and multiply by 
h~ 1 \ resulting in 

Proposition 8.1.5. (cf. proposition 4- 1-10). For the pilot-model operator 
with 1 < ii < h' 1 , 1 < k < 3 

(8.17) 11^-2^)1 < (8.8) + (8.18) as A 2 > fih, 
with 

(8.18) C^-h-^-^X^fih/X 2 ) 1 
and 

(8.19) |l7 A) - Zm| < (8.20) + (8.21) + (8.22) as A 2 < fxh 
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with 



(8.20) 
(8.21) 
(8.22) 



i-^A 3 (l + (log(^/A 2 )) + ) 2 , 
l-^A 3 (l + (log(/iVA 2 )) + ) 



Proof. Easy details are left to the reader. 



□ 



Case < k < 1. In this case we need to modify arguments of the previous 
subsubsection. 

(i) Consider case k' — k" = r = first. As we did not use A 2 3> [ih estimate 
(8.13) for contributions of such terms remains valid. 

(ii) Consider now case k' — k" = r ^ 0. Note that then according to (8.4) 
one should use ( = \k' + k"\ ■ (h\k - /c'|)l 

As A 2 > fih the last expression before summation with respect to k (= k') 
acquires factor (\k f — k"\\ 2 / jih)^ 1 and we get after summation with respect 
to k' 7^ 0, k" 7^ the same right-hand expression of (8.13). 

Meanwhile, as A 2 < fih we get before summation (as < k < 1) 



(8.23) C f i K+1 h- 2 \ 3 \k'\- 1 \k"\- 1 (\k' + k"\/\k' - k"\^Y~ K x 

{l + \k' -k"\\ 2 / l rhY l 



(with a logarithmic factor as k = 1). 

It is enough to sum with as 1 < k", k! = k" + r, r > 1; then the sum is 



r>l 



so we get 



(8.24) 




< K < 1, 



K 



= 1 



where case k = 1 we left to the reader. 



(iii) Consider now k' = 0, k" ^ 0; its tween case k' ^ 0, k" = is addressed 
in the same way. To do this we need to modify our expression fo I: namely, 
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we do not divide by t' and then take Fourier transform F t /_^-i T ; instead we 
take Fourier transform F t /_^-i r /, then integrate by t' to r and divide by ft; 
modifying this way (8.1) we arrive instead of (8.1) to the expression similar 
to (4.30). 

Then repeating the similar arguments of section 4 with the above modi- 
fication, we arrive to 

C/i k /T 3 A 3 |/cT Hk (1 + \k"\X 2 /fihy' 

instead of (8.23) and summation with respect to k" : \k"\ > 1 returns 
Cfx K h- 3 \ 3 {fih/\ 2 y as A 2 > fih and C^h^X^h/X 2 ) 1 ^ as A 2 < fih. 
So we arrive to 

, o f 1 < K < 1, 

(8.25) C(/xVA 2 )V K ^ 3 A 3 , , . M . 2 

[ (1 + l0g(/i/?)j K = 1 

as A 2 > fih and 

' 1 < K < 1, 



(8.26) C/J+^/rH«A 2+K 

as A 2 < fih. 

Thus we proved 



(l + l0g(/i/7)) 2 K=l 



Proposition 8.1.6. For the pilot-model operator with 1 < /j, < h \ < 

K < 1 

(8.27) |l[ A) -Z^ } | < (8.13) + (8.25) as A 2 > ///?, 
and 

(8.28) \\J X) - 2j5[)| < (8.13) + (8.24) + (8.26). as A 2 < ji/i 

8.2 General operators 

Consider now general operators satisfying conditions (0.5)-(0.6) and either 
(0.7) or (0.8) as r = 0. 

First, using proposition 7.1.1 and propagation results we conclude that 
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Proposition 8.2.1. (cf. proposition 4-2.1). Let conditions (0.5)-(0.6) be 
fulfilled. Then 

(i) Contribution of the zone {{x,y) '■ \x — y\ > e /U -1 } (before rescaling) to 
l T does not exceed Cfi K+1 h~ 2 and 

(ii) Under condition (0.7) contribution of time {\t\ > e } (before rescaling) 
to l T does not exceed Cfi K+1 h~ 2 . 

Note that 0(/i K+1 /?~ 1 ) either coincides or is smaller than the remainder 
estimate (8.14) and as long as we are interested in Weyl approximation X w 
we should be completely happy with it. 

Under condition (0.7) we are done arriving to statement (i) of theo- 
rem 8.2.3 below. 

Under condition (0.8) let us introduce the scaling function 

(8.29) lambda(x) = eolWF -1 ! + ^A, A = Cq^ 1 

and consider covering of 6(0, 1) by A-admissible elements. In virtue of 
proposition 8.2.1 (i) we need to consider only pairs {x,y) belonging to the 
same element. Due to the rescaling arguments, the contribution of one such 
element to |/ T — X w | does not exceed Cu,h~ 1 ~ K \ 3 and therefore 

(8.30) Under condition (0.8) contribution of the zone 

{(x,y):\x-y\<C of i-\ |WF _1 | < A} 

to |/ T — X w | does not exceed C/i/?~ 2 ~ K A 3 ; 

as A = A we get 0(fi~ 2 h~ 2 ~ K ). 

Therefore we need to consider only zone 

{(x,y):\x-y\<C of i-\ IVVF" 1 ] > A}. 

We are going to prove that 

Proposition 8.2.2. (cf. proposition 4-2.2). Let conditions (0.5)-(0.6) be 
fulfilled. Further, let fi < h -1 . Then for each z contribution of zone 
{(x,y) G 6(z,A(z)), \x-y\ < C At _1 |} with A(z) > Co// -1 to |I T -X W | does 
not exceed 



Chapter 8. Dirac energy: 3D -calculations 



123 



(i) /? W (A) which is defined as the right-hand expression of (8.8) or (8.13) 
for l</t<3 ; 0</t<l provided A > (fih)^ 6 ; 

(ii) the right-hand expression of (8.17) or (8.27) for 1<k<3 ; 0<k<1 
provided A > (yu/j) 5 / 

(Hi) the right-hand expression of (8.19) or (8.28) for l</t<3 ; 0</t<l 
provided A < (yu/j) 5 - 

Consider components of /? W (A). Note that expressions (8.8), (8.13) and 
(8.20)-(8.22), (8.24), (8.26) contain A in the positive powers while expressions 
(8.18), (8.25) contain A in the negative powers. Therefore after integration 
with respect A -1 dX we get expressions (8.8) and (8.13) as A = 1 (i.e. we get 
(8.14) and we get all other expressions as A = (/u/?)2; one can check easily 
that these latter expressions are less than (8.14) + C/)~ 2 ~ K . 

Therefore we arrive to 



Theorem 8.2.3. (cf. theorem 4.2.3). 
H< h~\ Then 

(i) Under condition (0.7) estimate 



Let conditions (0.5) -(0.6). Let 



(8.31) |l-l w | < CR w + Ch- 2 - K ; 
(ii) Under condition (0.8) 

(8.32) |l - l w | < CR W + CpSh-\- K - & + Ch- 2 -\ 

Proof of proposition 8.2.2. First of all applying proposition 2.2.2 and the 
same arguments as in the analysis of section 2 we arrive to 

(8.33) Contribution of the zone {| sin(0(t'))| > Qh, | s\n(0(t"))\ > C h} to 
the error |I T — X w | does not exceed /? W (A). 

Therefore we need to explore contributions of three remaining zones 

(4.43) {| sin(^(t / ))| < 2C h, \ si n (6*( t r/ ) ) | < 2C h} 
and 

(4.44) {|sin(^(t / ))| > 2C Q h, |sin(0(t"))| < C Q h} 
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(the same analysis would cover its tween zone { | sin {9(t')) | > 2C h, | sin(0(t"))| < 

c h}). 

We would like to follow the arguments of the proof of proposition 4.2.2 
but there is a problem: the magnetic drift is controlled by a and it may be 
much less than A. However let us apply few general arguments. 

Zone (4.43) as A 2 > fih. First, if we consider zone {(x,y) : |x — y| > 
£} (in the non-rescaled coordinates) we estimate the contribution of the 
corresponding pair (k' , k") by Ch 2 £~ K \k'\~ 1 \k"\^ 1 x h~ d X d and summation 
with respect to \k'\ < fiX, \k"\ < jiX returns 

(8.34) Cfi 2 h 2 ~ d r K {l + | \ogfiX\) 2 . 

Meanwhile, if we consider zone {(x,y) : |x — y| > £} (in the non-rescaled 
coordinates as well ) then estimating U(x,y, t) by h~ d (no rescaling, d = 
2, 3) we we estimate the contribution of the corresponding pair (k' , k") by 
Ch 2 x br 2d x l d ~ K X d and again summation with respect to k', k" returns 

(8.35) CiL 2 h 2 ~ 2d i d - K X d (\ + | log/iA|) 2 . 
Comparing (8.34) and (8.35) we see that their sum is minimal 

(8.36) Cu 2 h 2 - d - K X d (l + |log/iA|) 2 
as i = h. 

Comparing with (8.8) and (8.13) we conclude that for d = 3 this is less 
than /? w as < k < 1, A > fih. To deal with 1 < k < 3 let us note first 
that there must be a factor 

(8.37) C(l + \k' - k"\X 2 /fihy' 

before summation. Really, as |^3 + 773 1 > e A this factor is due to the shift 
in (x 3 +y 3 ) and as + n 3 \ > e A but \VV/F\ x A or \V X3 V/F\ x A this 
factor is due to the shift in (£1 + r/i) or (£ 3 + r? 3 ) respectively. 
Therefore 

(8.38) As A 2 > fih 1 " 5 contribution of k' ^ k" is negligible. 

This allows us also eliminate logarithmic factors in (8.36) and (8.37). 

Let us apply cut-off Q by hD X3 and Q' by hD xy3 . Assume first that either 
I £3 1 < Cp on the support of (the symbol of) Q or 1 773 1 < Cp on the support 
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of (the symbol of) Q', in our estimate we acquire factor p and we get an 
expression below /? W (A). 

Therefore we can assume that 

(8.39) |£ 3 | > p := max(( / u/?)5, p^ 1 ) on the support of (the symbol of) Q 

(or 1 773 j > max(( j u/))i/i _1 ) on the support of (the symbol of) Q'. 

Let us consider corresponding classical dynamics. Assumption (8.39) 
implies that in k winding it travels the "observable" distance /c£ 3 and if 
there was no "turn back" the contribution of k-th winding to estimate would 
not exceed 

cp 2 h^\kr 2 ^ [ i^r^sA 3 
j\&\>p 

where integral is x 1 as < k < 1, x logp as k = 1 and x p 1 ~ K as 
1 < k < 3. Obviously multiplied by p 2 h~ 1 it does not exceed /? W (A) and 
therefore summation with respect to k returns 0(/? w (A)). 

Now we need to consider effect of the "turn back". Let (3 = |V X3 V / F\. 
Note first that as |£ 3 | > (5X there is no turn back for time \ t\ < e\. Therefore 
we need to consider only |£ 3 | < /3X. 

There is come back on k-th winding as ^ 3 is close to (*(k) ~ constfip^ 1 k 
and the distance |x — y| x /i _1 |/c||^ 3 — C(^)l an d then contribution of k-th 
winding (as |x — y| > Ch) does not exceed 

Cn 2+ «h- x \k\- 2 -" [ |6 - C(k)\- K d£ 3 X 3 < 

3 J p 2+K h- 1 \k\- 2 - K 0</t<2, 

\p% l ~ K \k\- A 2<k<3, 

with a logarithmic factor (1 + log \ k\ + \ log ph\) as k = 2. The sum of the 
expression above is the same as this expression for k = 1 (we denote it 
temporarily by RX 3 ) which obviously less than /? W (A). 

Finally, we need to consider the contribution of k-th winding as \x—y\ < h. 
However it is done by the standard way as U(x, y, t) after cut by £ 3 and ry 3 
acquires factor Ch 2 . 

Therefore, under assumption A 2 > ph 1 ^ 5 we are done. 

Case X 2 > ph. As A 2 > ph (but A 2 ^> ph) terms with k' 7^ k" are no more 
negligible but in the estimate we can ignore due to the factor (8.37) them 
except when k' = or k" = 0. 
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Applying the same arguments as in 2D-case i.e. replacing in term with 
k' = factor |/c'| _1 by factor Fr 1 we arrive to RX 3 x l/(u,h) x (/i/)/A 2 )' which 
is less than term (8.18) or (8.25) asl</t<3or0</t<l respectively. 

Zone (4.43) as A 2 < fih. In this case we apply the same arguments but 
now effectively we sum over k', k" with \k' — k"\ < fih/X 2 which effectively 
boils out to the summation over k' = 1, k" = 1, ... , u.h/\ 2 or k' = 0, 
k" = 1, ... ,fih/\ 2 ; the former results in RX 3 x fih/X 2 and the latter results 
RX 3 x u-h/X x l/(/i/)) and one can see easily that the result is less than what 
we claimed. 

Zone (4.44) as X 2 > fih. Now there is always shift with respect to (x — y) 
which is at least x | sin(t')| and therefore as 1 < k < 3 contribution of this 
zone does not exceed contribution of zone (4.43). As < k < 1 we get 
C/j, K+1 h~ 2 X 3 \k\~ 2 and summation returns Cu. K+1 h~ 2 X 3 which is the answer 
as A 2 3> let us denote it by RX 3 ; for k = 1 it includes logarithmic factor 
(l + |log(^)|). 

The same answer estimates contribution of k' ^ 0, k" ^ as A 2 > 
fiti; contributions of k' = 0, k" ^ and k" ^ 0, k' = do not exceed 
RX 3 x l/(yu/?) x (fih/X 2 ) 1 and one can see easily that the result is less than 
what we claimed. 

Zone (4.44) as X 2 < fih. Again as 1 < k < 3 contribution of (4.44) does not 
exceed contribution of (4.43) and as < k < 1 we estimate contribution 
of k' ^ 0, k" ^ by RX 3 x fih/X 2 , and contributions of k' ^ 0, k" = and 
k' = 0, k" 7^ by RX 3 x fih/X 2 x l/(yu/?) and one can see easily that the 
result is less than what we claimed. □ 

8.3 Perturbations 

Consider perturbations. However due to the presence of x 3 , hD 3 there are 
rather different components here than in subsection 4.3. 

8.3.1 Perturbations. I 

Let us start from the terms 0{\x' — y'\). During evolution these terms are 
0(/i _1 ) and therefore we need to multiply by 0(/i _2 /? _1 |/c|). 
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1 < re < 3. Consider first case k' = k" = k. Assume that A 2 > fiti. 
Recall that k-th tick' contribution to l T was not exceeding expression (8.7) 
multiplied by /i K /?~ 3 A 3 |/c| -3 i.e. 

^ K+ i h -2 X i- K 1 < re < 2, 

(8.40) CX 3 \k\- 3 { o 3 h- 2 (l+ |log(n/A 2 )|)A- 1 re = 2, 

//3/j -« A -3+ K 2 < re < 3 

(as |sin(t')| < ft or |sin(t")| < ft one should use proposition 8.2.2) and 
multiplying by Cfi~ 2 h~ 1 \k\~ 1 we get after summation with respect to k : 
\k\ > 1 

r / / k - 1 /7- 3 a 1 - k i < k < 2, 

(8.41) CA 3 J / u/?- 3 (l + |log(n/A 2 )|)A- 1 re = 2, 

[ fih^- K \- 3+K 2 < re < 3 

and summation with respect to A returns the same expression as A = 1 
which does not exceed Ch~ 2 ~ K . 

Next we need to consider ^ /t' ^ /t" ^ 0. Recall that the contribution 
of such pair to l T does not exceed 

(8.42) C\ 3 \k'\- 3 2\k"\- 3 2(u.h/\ 2 \k' - k"\)'x 

^ K+ i h - 2> i- K 1 < re < 2, 

/i 3 /T 2 (l + |log(///j/* 2 )l)* _1 K = 2. 

H 3 h- K X- 3+K 2 < re < 3. 

Multiplying by Cfi^h^lk'l (as the perturbation goes to the first "factor") 
we get 

(8.43) C\ 3 \k'\- 1 2\k"\-hvh/\ 2 \k' - k"\)'x 

^ K -i h -3 X i- K 1 < re < 2, 

/i/T 3 (l + llog^/A^DA- 1 re = 2, 

A t/7- 1 - K A- 3+K 2 < re < 3. 

Then summation with respect to k' of \k'\~i(fih/\ 2 \k'— k"\) ! returns \k"\^^(fih/X 2 )' 
as A 2 > fih and summation with respect to k" then returns the same expres- 
sion as we got considering k' = k" ^ with an extra factor (/x/j/A 2 )'. 
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On the other hand, as A 2 < fih one should replace (8.7) by C(/i/?)2 2 K ; 
multiplying by /^/?- 3 A 3 |/c'|-i|/c"|-i(l + \ 2 \k' - k"\/ u,h)-> x ^- 2 h~ l \k'\ we 
get 

(8.44) C/i-^ + 5 K /)-f'^A 3 |/c / |-5|/c // |-i(l + X 2 \k' - k"\/u,h)-' 
and one can see easily that summation with respect to k' , k" returns 

Qi-=:+H-§-2- K A 3 (/i/>/A 2 ^ x C^ K h- 2 -^ K \ 2 

and summation with respect to A returns C jA K h~ l ~^ K which again is less 
than Ch- 2 ~ K . 

3 3 

If k' ^ 0, k" = we recall that \k"\~2 should be replaced by h~? in 
(8.42)-(8.43) and perturbation brings factor /j, 2 /7 1 1 /c' | so we get instead of 
(8.43) 

C\ 3 \k'\- 1 *u.-h--*(u.h/\ 2 \k , \)'x 

^ K -I h -3 X 1-K 1 < K < 2, 

/i/T 3 (l + \\og{fi h/X^DX- 1 k = 2, 
fih' 1 -^- 3 ^ 2<k<3 

and summation with respect to \k'\ and then A > (fj,h)^ returns Coh K ~^ \T^^- K . 
Similarly, as A 2 < fih we start from (8.44) which is replaced by 

C>- 2+ H- 4 -2- K A 3 |/c'|^(l + X 2 \k'\/u,h)-' 

and summation with respect to k' and returns C ji~\ + h K h~ 7 z~^ K X 2 and 
then summation with respect to A < (/i/)) 5 returns the same answer 

As \k'\ < \k"\ the estimate of the contribution of the pair (k' , k") to an 
error is obviously less than the same estimate for the pair (k", k'); this takes 
care of k' = 0, k" ^ 0. Finally for k' = k" = we can use the standard 
arguments to estimate the same way. 

As the result we conclude that the error will not exceed Ch~ 2 ~ K . 

< k < 1. The similar arguments work for < k < 1. However one must 
reconsider zone {(x,y) : |x — y| > Cq^T 1 } where contribution of the pair 
(A 7 , k") to Weyl asymptotics does not exceed 

(8.45) C/i 2 /)- 2 C 1_K |/cr 3 A 3 x Cfi K+1 h- 2 \k\- 2 - K X 3 
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with C x /x 1 1 ^ | as k' = k" = k ^ 0. Summation with respect to k : \k\ > 1 
returns then Cfi K+1 h~ 2 \ 3 as it should. 

As /c' = /c" = /c 7^ it results in the approximation error not exceeding 

(8.46) CAt 2 /7- 2 C 1_K |/f| _3 A 3 x fi^h-^kl x C// K - 1 /7- 3 |/c|- 1 - K A 3 

and summation with respect to k, A returns Cfi K ^ 1 h^ 3 . 

As 7^ /c' 7^ /c' 7^ 0, | k'\ > \k"\ it results in the approximation error not 
exceeding 

(8.47) C f i 2 h- 2 (u,- 1 \k'\)) 1 ~ K \k / \-l\k"\- 3 2x 

(1 + X 2 \k' - k"\/fih)-'X 3 x ^h-^k'l x 

C// K - 1 /?- 3 |/f'|i-' t |/f"|-i(l + A 2 |/c' - k"\/iihy'\ 3 

and the sum with respect to (k', k") and A : A 2 > fih does not exceed 
Cfi^h- 2 . 

Meanwhile as A 2 < fih the sum with respect to (k', k") does not exceed 
C fi^ 1 h~ 3 (u.h / \ 2 )2~ K \ 3 and the sum with respect to A is less than Cn. K ~ x br 3 . 
As k' ^ 0, k" = (8.47) becomes 

(8.48) C / u K - 1 /)- 3 |/c / |5- K ( /[ x/))-i(l + \ 2 \k'\/iih)-'\ 3 

and then summation with respect to k' returns C /i K_1 h~ 3 (fj,h)~^(fih/ X 2 ) 1 A 3 
as A 2 > fih and then the sum with respect to A is less than C^' 1 ^ 3 . 

Meanwhile as A 2 < fih the sum with respect to k' does not exceed 
C/x K ^ 1 /) _3 (/i/)/A 2 )2 _K (/i/)) _ 2 A 3 and then the sum with respect to A is less 
than C/x K_1 /7~ 3 . 

As \k'\ < \k"\ the estimate of the contribution of the pair (k' , k") to an 
error is obviously less than the same estimate for the pair (k" , k'); this takes 
care of k! = 0, k" ^ 0. Finally, for k' = k" = we can use the standard 
arguments to estimate contribution by Ch~ 2 ~ K . 

For K = lwe apply the same arguments with the obvious modifications, 
concluding that the error does not exceed Ch~ 3 (l + | log fih\). 

Note that in contrast to the case 1 < k < 3 now the error is larger than 
Ch- 2 - K . 

To improve this remainder estimate let us notice that the second term 
in the approximation is and the third term also acquires factor 2 /7 1 1 A" | 
(as k' = k" = k 0); so we get 

Cfi K - 1 h- 3 \k\- 1 - K x /x 2 /? 1 1 /c| A 3 x Cfi K - 3 h- 4 \k\- K X 3 
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which we should sum as \k\ < fi 2 h; we also should sum 

C// 2 /7- 2 C 1- *|/cr 3 A 3 x Cfi K+1 h~ 2 \k\' 2 ' K \ 3 

as \k\ > fi 2 h and both sums return C / u K+1 /7~ 2 ( / u 2 /7)~ 1 ~ K A 3 x C/i~ 1 ~ K h~ 3 ~ K \ 3 . 
Summation with respect to A returns C[i~ l ~ K h~ 3 ~ K . One can estimate 
contribution of all all other pairs (k', k") in the same way. As k = 1 we get 
C>- 2 /r 4 (l + \\o gf ih\). 
So we arrive to 

Theorem 8.3.1. (cf. theorem 4-3.6). Let conditions (0.5)-(0.6) and either 
(0.7) or (0.8) be fulfilled. Leth'^K h' 1 . Then 

(i) As 1 < k < 3 

(8.49) |l - 1° | < Ch- 2 - K 
where 

(8.50) l°:= J Q 1 ^(x + y),x-y)-\el {x+y) (x,y,0)\ 2 dxdy 

and e°(x, y, r) «s the Schwartz kernel of the spectral projector of the general- 
ized pilot-model operator 

(8.51) A° z := h 2 D 2 + (hD 2 - + h 2 D\ + \/ (x), V°{x) := l/fa, z 2 , x 3 ) 

which in contrast to (6.36) V° does not contain linear terms with respect to 
(*i,x 2 ); 

(a) As o < k < i 

(8.52) |l - 1° | < Ch- 2 - K + Cfi- 1 - K h- 3 - K 
and as k = 1 

(8.53) |l - 1° | < Ch~ 3 + C>- 2 /T 4 | \ogfih\. 

Remark 8.3.2. Note that the remainder estimate does not exceed Ch~ 2 ~ K 
for /i > /j-VU+k) as < /t < 1 and for > /?~2(1 + | log pt/7|) 5 as k = 1. 
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8.3.2 Perturbations. II 

If we go from the generalized pilot-model approximation to the magnetic 
Weyl approximation then we need to counter (x 3 — y 3 )X terms and this is 
0(fi~ 1 \k'\\) due to the finite speed of propagation. Then dealing as before 
we acquire an extra factor \k'\\. 

1 < k < 3. Then as k' = k" = k ^ we have \k\~ l rather than \k\~ 2 and 
the sum with respect to k contains an extra logarithmic factor (l + log(yU 2 /7))A 
in comparison to what we got in subsubsection 8.3.1 "Perturbations. I". 

Consider case ^ k' ^ k" ^ 0, \k'\ > \k"\. Then as A 2 > fih we arrive to 
the same expression (but with an extra factor (f^h/ A 2 )'). However as A 2 < fih 
we need to consider \k'\~2 |/c"|~2 x \k'\ 2 \ with k' , k" : \k' — k"\ < fxh/\ 2 which 
in comparison to what we got before after summation returns an extra factor 
u.h/\ 2 x A and the sum with respect to A does not exceed what we got above. 

Consider ^ k', k" = 0. Again no change as A 2 > \ih and an extra factor 
(/i/?/A 2 ) x A as A 2 < fih and the sum with respect to A does not exceed what 
we got above. 

Therefore, as a result we get Cfx K ~ 1 h~ 3 (l + | \ogfi 2 h\) as 1 < k < 2 and 
Cfih~ 1 ~ K (l + | logyU 2 /?|) as 2 < k < 3 with an extra factor (1 + | \ogfxh\) as 
k = 2. In both cases it is less than Ch~ 2 ~ K unless /j, > /7 _1 | log h\~ r . 

To improve our result in this latter case we apply three term approxi- 
mation and note that the second term is while the third term contains 
summation of /? w A 3 |/c|~ 3 x (yU -2 /)|/c| 2 A) 2 for k : \k\ < fih? and summation 
of R w \ 3 \k\- 3 for k:\k\> u,hl as k' = k" = k ^ where R w = fx 3 h~ K as 

2 < k < 3, /? w = n R+1 h~ 2 as 1 < k < 2 with an extra factor (1 + | log fih\) 
as k = 2. 

Other pairs (k' , k") should be modified in the same way and we leave 
easy but tedious details to the reader. As a result we get /? w /i" 2 /) -1 without 
any logarithmic factor thus arriving to the same estimate as before but 
without logarithmic factor; as 1 < k < 3 it does not exceed Ch~ 2 ~ K . So we 
arrive to statement (i) of theorem 8.3.3 below. 

< k < 1. Now expression (8.46) acquires factor \k\\ and becomes 
C/j, K ~ 1 h~ 3 \k\~ K \ 4 and summation with respect to k : \k\ < fih2\~2 returns 
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We also should sum Cjj, K+1 h 2 \k\ 2 K A 3 for \k\ > /ih^X 2 with the same 
result. 

5 1 

Finally summation with respect to A returns Ch~z~2 K . 

Other pairs (k' , k") should be modified in the same way and we leave 
easy but tedious details to the reader. Three term approximation does not 
improve results. As k = 1 we acquire factor (1 + | log h\). So we arrive to 
statement (ii) of theorem 8.3.3 below. 

Theorem 8.3.3. (cf. theorem 4-3.4)- Let conditions (0.5)-(0.6) and either 
(0.7) or (0.8) be fulfilled. Let h~\ < u, < h' 1 . Then 

(i) As 1 < k < 3 

(8.54) ||-J MW | < Ch- 2 - K ; 

(ii) As < k < 1 

(8.55) |l -X MW | < Ch^-h 
and as k = 1 

(8.56) |l -X MW | < C/T 3 (l + \\ogfih\). 
8.3.3 Perturbations. Ill 

Now in the case < k < 1 we consider the pilot-model approximation. 
As (8.55) is larger than (8.53) we consider improvement of the previous 
subsubsection only. Now as perturbation is (x 3 — y 3 ) 2 which translates 
into fi~ 2 \k'\ 2 rather than expression (8.46) acquires factor uT 1 ^'^ and 
becomes Cfi K ~ 2 h~ 3 \k\ 1 ~ K \ 3 and we need to sum this as \k\ < /i/?3 resulting 
in 

C/j, K - 2 h- 3 {fih^) 2 - K \ 3 x C/T5-5 K A 3 . 

We should sum Cfi K+1 h~ 2 \k\~ 2 ~ K \ 3 for \k\ > jih^ with the same result. 

Finally, summation with respect to A returns Ch~^~3 K . 

Other pairs (k' , k") should be modified in the same way and we leave 
easy but tedious details to the reader. Three term approximation does not 
improve results. As k = 1 we acquire factor (1 + | \ogh\). So we arrive to 
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Theorem 8.3.4. Let conditions (0.5) -(0.6) and either (0.7) or (0.8) be 
fulfilled. Let h~2 < u. < h' 1 . Then as < k < 1 

(8.57) |l -T| < Ch- 2 - K + C/)-J-5 K + Cfi- 1 -^- 3 -* 
and as k = 1 

(8.58) |l - T| < Ch~ 3 + Cfi- 2 h-\l + | \og fih\) 
where I is defined for a simplified pilot-model operator 

(8.59) A z := h 2 Df + {hD 2 - fi Xl ) 2 + h 2 Dj + V{x), 

V(x) := V(z) + d Z3 V(z)(x 3 -z 3 ). 

8.4 Superstrong magnetic field 

Consider pilot-model Schrodinger-Pauli operator (1.72). Recall that accord- 
ing to propositions 7.2.1 (ii) , 7.2.2(H) respectively 

{h~ K as < k < 1, 

+ {\og/ih) + ) as « = 1, 

^l^-i) /,-!(«+!) as i < K < 3 

and 

(7.3) |I-I T |< 



/ /? K as < k < 1, 

/? _1 (l + (log;u/?) + ) as k=1, 



C>lH * ^l^-i) /,-§(«+!) as i < K < 3, K ^ 2, 

jj,ih~i(l + (log/i/?) + ) as k = 2. 

Let us assume that F = 1 to avoid some unpleasant correction terms as 
reduced to F = 1. 

(i) Consider the generalized pilot-model approximation (8.51). Then o x h 
and under condition 

(8.60) \V + 2mfih\ > e Vm G Z+ 
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the shift with respect to X3 is x |t| which implies that we need to consider 
only t' x t" as the total contribution of other intervals will be less, and 
that contribution of intervals \t'\ < C h, \t"\ < C h does not exceed the 
right-hand of (8.2) multiplied by Ch, and that contribution of intervals 
t' x T, t" x T (T > C h) does not exceed Cfih^ 1 x T~ K so summation 
with respect to T~ 1 dT returns Ciihr 1 ~ K . 

Then the error does not exceed the right-hand expression of (8.3). 

(ii) Consider the magnetic Weyl approximation. In this case a x T and 
again under assumption (8.60) we can consider only t' x t", and we estimate 
contribution of intervals \t'\ < C h, \t"\ < CqH as above, but contribution 
of intervals t' x T, t" x T [T > C h) does not exceed C/i/? -2 x T l ~ K 
and summation with respect to T^dT returns Cfih~ 1 ~ K as 1 < k < 3. As 

1 3 1 

< k < 1 we sum this way only for T < resulting in C/j,h~2~2 K and 
we sum Cfih' 1 T~ K for T > with the same result. As k = 1 we get 
Cfxh~ 2 {l + \\ogh\). 

(iii) Consider the pilot-model approximation as < k < 1. In this case 
a x T 2 (as we pass from the generalized pilot-model approximation) and 
again under assumption (8.60) we can consider only t' x t", and we estimate 
contribution of intervals \t'\ < C Q h, \t"\ < C h as above, but contribution 
of intervals t' x T, t" x T (T > C h) does not exceed Cfih~ 2 x T 2 ~ K 

1 4 1 

summation with respect to T dT returns Cu.h~3~3 K as we sum with 
respect to T _1 dT for T < hi; for T > /?5 we sum Cfih~ 1 T~ 1 ~ K with the 
same result. 

Then we arrive to 

Theorem 8.4.1. (i) For the Schrddinger-Pauli operator under conditions 
(0.5)-(0.6) and (8.60) 



(8.61) |I-X MW |<C^ 



' fxh~^~^ K 0<k<1, 
/i/T 2 (l + I \ogh\) K = l, 

^ K+ l) h -^ K+ 3) 1< K < 3, K ^ 2, 



{ fiih-2{l + {\ogfih) + ) K = 2. 

(ii) Further, as < k < 1 

rih~ 1 - K < k < 1, 

^ 2 (1 + (log ///?) + ) K=l 



(8.62) |I-I°|<C 
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and 



(8.63) 



I - II < c 



\ih 1 K + \ih 3 3 
/i/7- 2 (l + |log/7|) 



< K < 1, 



K = 1 



Let us get rid off condition (8.60). First, we still can assume without any 
loss of the generality that (8.60) holds for all m £ Z + save m = m. Then 
making ^-admissible partitions in x and y with 



one can replace (8.60) by (0.7) and then making ^-admissible partitions in x 
and y with 



one can replace (0.7) by (0.8) arriving to 

Theorem 8.4.2. Theorem 8.4-1 remains valid with assumption (8.60) re- 
placed by (0.8). 

Easy but tedious details are left to the reader. 

8.5 Problems and remarks 

Remark 8.5.1. The main difference between cases 2 < k < 3, 1 < k < 2 and 
< k < 1 that in the first case the main contribution to the remainder is 
delivered by (x, y) close to one another (|x — y| <C /i _1 ), in the second case 
by (x,y) with |xy| x iC 1 and in the third case by (x, y) with |xy| x 1. 

We can calculate X MW and I plugging corresponding expressions e^ w (x, x, 0) 
and e x (x, x, 0) into I. 

Problem 8.5.2. Find nice expressions for X MW and\. 

Problem 8.5.3. As fih < 1 get rid off condition (0.5); according to Chap- 
ter 13 we do not need it for estimate \ \ — l T | but we want to get rid off it in 
estimates for \\ -Z w |, |l - J MW | and |I-T|. To do this 

(i) Assume first that condition (0.7) is fulfilled and consider the scaling 
function £(x) = max(e| V(x)|, fib, C/i _1 ) . 



(8.64) 



i = e\V + 2mfih\ + hi 
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(ii) Assuming then that condition (0.8) is fulfilled and consider the scaling 
function £{x) = max(e(| V(x)\ + |W(x)| 2 )5, (fih)K C/i" 1 ). 

Problem 8.5.4. Weaken non- degeneracy condition (0.8) which seems to 
be excessive here (in contrast to 2D case). 

Problem 8.5.5. Get rid off condition (0.6) assuming instead that |F| + 
|V® F| > e. 

9 Estimates of expression (0.4) 
9.1 Basics 

In the multiparticle quantum theory we need to calculate expression (0.2) 
approximating energy of the electron-electron interaction or, more often, to 
estimate expression 

(0.4) K : = 



where u(x,y) = u> K (x,y) satisfies assumption (0.3) and M is some approxi- 
mation. We concentrate on (0.4); asymptotics of (0.2) could be proven by 
the same arguments and we leave them to the reader. 

Applying partition of unity we note that both of them depend on asymp- 
totics of expression 



with 7-admissible function vb^. 

Let us denote by 7 d R(a, 7) an upper estimate of an error occurring when 
we replace in (9.1) e(x,x, 0) by h M and let 




(e(x, x, r) - h- d Af x (r)) (e(y, y, r) - h~ d Af x (r))oj(x, y) dxdy 



(9.1) 




(9.2) 



M = Ch~ d 



+ Cjih 



l-d 



be an upper estimate for e(x,x, 0). 
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Proposition 9.1.1. Under condition (0.7) 

(9.3) |K| < C J (/?(l, 7 ))V- 1 - K of 7 
and 

(9.4) | J - J\ < CM J R(l, i)i d - 1 - K d 1 
and under condition (0.8) 

(9.5) |K| < C II (/?(«, 7 ))V" 1 ~ K tt d ~ 1 djda+ 

JJ{~i<a} 

C II R{a, a)R{p, P)a d ~ l p d ~ l ~ K dadfi 

J J{a<f3} 

and 

(9.6) | J - J\ < CM II R{a, 1 ) 1 d - 1 - K a d ~ l d 1 da+ 

J J {"/<a} 

CM II R{a,a)a d - 1 p d - 1 - K dad (3 

J J{a<(3} 

where for a < C 0/ u _1 we need to modify a~ 1 R(a,'j) (see below) 22 \ 

Proof. Due to representation (3.1) and definition of R(a, 7) and M contribu- 
tion of the zone {(x, y) : |x — z| x 7, \y — z| x 7} (note that R(a, 7) < M) 
to K does not exceed C(R(a, 7)) 2 7 2rf ~ K for fixed z and 7 where a = a(z) 
provided 7 < 3ea. 

Under assumption (0.7) summation over z results in C(/?(l, 7)) 2 7 d ~ K 
and summation over 7 results in estimate (9.3); estimate (9.4) is proven in 
the same way. 

On the other hand, under assumption (0.8) summation over z : a(z) x a 
results in C(R(a, r y)) 2 a d 'y d ~ K and summation over 7, and then over a results 
the first term in the right-hand expression of (9.5); it estimates contribution 
of zone {(x, y) : |x — y| < 3eaa(x)} (in which a(y) x a(x)). 

22 ) As d — 3 we denote by a what we used to denote by 7 in sections 5-8 (i.e. (a 2 + f3 2 )^ 
in these sections notations. 
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Further, under assumption (0.8) if |x — y| > 2ea, a := a(x), (3 := a(y) 
then |x — y| x max(a, (3). Therefore as (3 > a the contribution of zone 

{(x, y) : |x - y| x /3, a(y) x a, a{y) x /?} 

does not exceed CR(a, a)a d x R(/3, /3)/3 d ~ K and summation over a, /3 returns 
the second term in the in the right-hand expression of (9.5); it estimates 
contribution of zone {(x,y) : |x — y| > e(a(x) + a(y))}. □ 

9.2 Case d = 2 

9.2.1 Tauberian asymptotics 

(i) Consider first case /i < /7 -1 . Now according to (1.58) we need just to plug 

(9.7) /?(a, 7 ) = ^ l h^ + a- 1 R T {a, 1 ) 

as h~ 2 Af x = e T (x,x, 0). Obviously the effect of C/i -1 /? -1 term to both (9.3) 
and (9.5) will be term CiT 2 br 2 . 

Plugging into (9.3) or (9.5) R(a, , ~f) = a~ 1 R T (a, 7) with R T (a, , ~f) defined 
by (1.57) we obviously eliminate from integration the last two zones (defined 
in (1.57)) as a power of 7 is negative 23 ^ and the first one as a power of 7 is 
positive. Thus 7 snaps to its minimal value in the last two zones and to its 
maximal value in the first one, and only the second zone matters. 

Further, in virtue of (1.57) R T (a, , ~f) = Cfi2 / ~f~2 in the remaining second 
zone {h < 7 < m\n(fih/a, fi" 1 )}. 

Then the power of 7 is positive, or negative as0</c<l,/c=lor 
1 < k < 2 respectively. Thus 7 snaps to m\n(fih/a, a~ l ) as < k < 1 and 7 
snaps to h as 1 < k < 2. 

Then (9.3) becomes 

{fih 1 -" 1<k<2, 
//min(|log///7|,|log//|) k=1, 
/ umin((/i/?) 1 - K ,/i- 1+K ) 0<k<1. 

Consider the first term in the right-hand expression of estimate (9.5). 
If 1 < k < 2 as we mentioned 7 snaps to h thus resulting in the integrand 
C/i/? 1_K and in C / u/? 1_K | log h\ after integration with respect to a~ l dct. 

23 ) Here and below talking about powers of a or (3 or 7 we mean in the integral with 
respect to -f~ 1 d-f or a~ 1 da or (3^ 1 df3. 
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Let < k < 1. Then 7 snaps to 7(a) := m\n(fih/a, /i^ 1 ), thus re- 
sulting in the integrand C/i min ((/i/?/o;) 1 ~ K , /i K ~ 1 ) as < k < 1 and in 
C/xlogmin(/i/a, (/i/?)" 1 )) as /c = 1. 

Then after integration with respect to or^da we get 0(/i K | log /7| ) as 
< k < 1 or 0(/i K | log /j| • I log(/i/?)|) as k = 1 and this is less than the 
estimate of the contribution of zone (9.9) below unless /i > h~ l \ log 

However in (9.5) we need also calculate the contribution of the zone 

(9.9) {(x,y) : a{x) < C Q u.~\ a{y) < C // _1 } 

which is not covered by our arguments and also the second term. 

The former obviously does not exceed Cfi 2 h~ 2 x ur A+K = C[i~ 2+K h~ 2 
where the second factor is J \x — y\~ K dxdy over zone (9.9) and this also 
covers zone {(3 < C /i -1 } in the second term of (9.5). Obviously C[i~ 2+K h~ 2 
is larger than anything we got before. 

Consider the second term in the right-hand expression of (9.5). Here we 
have R T (a, a) = a' 1 as (fih)2 > a > fi^ 1 and fi^ 1 a^ 1 h^ 1 min [(u. 2 h/a)2 , l) x 

(fih/a) 1 as < a < 1 and therefore (3 is in the negative power and 

snaps to a; so we arrive to 

J (/? T (a,a))V- K da. 

In this integral a is in the negative power and snaps to yU -1 resulting in 
(/? T (yU _1 , yU -1 )) fx K ~ 2 = 0(/i K ) which is less than we have already. 
Finally, we need to consider contribution of zone 

(9.10) {(x,y) : a(x) < C Q ^\a{y) > Qfi' 1 }. 

So, we get C/i~ 2 x fih" 1 j R T ((3, (3)(3 1 ~ K df3; in this integral (3 is in the 
negative power for sure as 1 < k < 2 and it snaps to (3 = /i -1 resulting in 
C fi' 1 h^ 1 R T (fi' 1 , /i~ 1 ) / u K ~ 2 which is less than we got already. 

As < k < 1, (3 is also in the negative power except zone {^T 1 < (3 < 
(/i/?)^} which is possible as fi > h~3\ then as (3 snaps to (/i/)) 5 resulting in 
C(fih)~2~2 K and one can see easily that it is less than what we got already. 
The same is true for k = 1 as well. 

So we arrive to estimate 

(9.11) |K T | < Cfi- 2+K h- 2 + C|log%( / u K - 1 + /? 1 - K ). 
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(ii) Let /i > h^ 1 . Now we need just to plug R(a, r y) = a~ 1 R T (a, 7) > 1 with 
/? T («,7) delivered by (1.76). Then in (9.3) we need to snap 7 to ^"2/72 and 
set /?(7) = 0.2 h~ 2. So (9.3) becomes 

(9.12) |K T | < Cfi^h- 1 ^. 

Meanwhile, under assumption (0.8)' in the first term of the right-hand 
expression of (9.5) a is in power and logarithmic factor appears; we get 
Cu.2 K h~2 K \ log u.\. 

Further, in (9.5) we need to estimate also contribution of zones 

(9.13) {(x,y) : a(x) < C ^hKa{y) < Qfrh^}, 

(9.14) {(x,y): a(x) < Q^hK a(y) > Q^h^} 

and of the second term in the right-hand expression and one can see easily 
that these contributions do not exceed Qi2 K /?~2 K | \og/j,\; and so we get 

(9.15) |K T | < C>H-5 K | \ogn\. 
Thus we proved 

Proposition 9.2.1. (i) For ji < br 1 under condition (0.7) estimate (9.8) 
holds and under condition (0.8) estimate (9.11) holds; 

(ii) Let ji > h~ x . Then under condition (0.7)' estimate (9.12) holds and 
under condition (0.8) estimate (9.15) holds. 

Problem 9.2.2. Prove that estimates (9.8) and (9.12) hold under condi- 
tions (0.8) + and (0.8) +/ respectively. 

9.2.2 Weyl approximation 

Now we should add to what we got for K T corresponding expressions with 
/?(a,7) = /? w (a,7) calculated according to (1.80)-(1.81) and therefore 7 
must be snapped to min(o;, fiha^ 1 ). Then (9.3) becomes 

(9.16) |K W | < Cfi- 2 h- 2 + C{R w ) 2 {iih) 2 - K = 

fi 4 h~ 1 (/ih) 2 ~ K as /i < /H, 



Cfi- Z h- Z + C 



/i 2 h- 2 (fih) 2 - K as /j-i < n < h- 1 
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which obviously is larger than the right-hand expressions of (9.8) and so we 
can neglect the Tauberian error. 

Consider the first term in the right-hand expression of (9.5). Note that 
for a < (/i/?) 5 7 snaps to a and then a is in the positive power. Also 
note that as a < fi 2 h (1.80) defines R w = jih' 1 and a is in the positive 
power as well. Therefore as h~2 < fi < h' 1 we snap a to 1 resulting 
in Cu 2 h~ 2 (fih) 2 ~ K prescribed by (9.16) and /i < h~2 we need to consider 
interval max(/z _1 , fj, 2 h) < a < 1; then automatically a > (fih)2 and we arrive 
to 

(9.17) C J fi 4 h^{fih/a) 2 ' K da. 

Here as 1 < k < 2, a is in the positive power and snaps to 1 resulting in 
C/i 4 /?~ 1 (/i/?) 2 " K prescribed by (9.16). 

On the other hand, as < k < 1, a is in the negative degree and 
snaps to max( / u _1 , fi 2 h) resulting in Cfi 4 h^ 1 (fih) 2 ^ K fi 1 ^ K for fi < and 
C /x 4 h~ 1 (fih) 2 ~ K (a 2 for h~3 < fi < and while the former expression 
is less than C fi~ 2+K h~ 2 the latter is not (at their respective intervals). 

As k = 1 we get C/i 5 (l + | log/i 2 /)|). 

Consider the second term in the right-hand expression of (9.5). Then 
both a < (3 effectively are bounded from above by (fih)z which leaves interval 
empty for /i < h~3. On the other hand then (3 < jji 2 h for h~3 < u, < h^ 1 
and we arrive to 



Cfi 2 h- 2 [ [ i aP^dadfi x C^h^ifih) 2 - 1 ^ 

J J{a<P<(tih)?} 



which is less than the right-hand expression of (9.16) unless < fi < br\ 
and < k < 1 in which case it is less than C // h~ 1 (/j,h) 2 ~ K (a 2 h)^ 1 . 
Therefore we arrive to 

9.18 K w < Cu.~ 2+K h~ 2 + C I P VP ; P T + 

"~ \fi 2 h- 2 (fih) 2 ' K as h-iKfiKh- 1 

/i 4 /)- 1 ^/?) 2 -^/! 2 /?)^ 1 0<k<1, n<h~K 

C ( fi 5 (l + \\ogfi 2 h\) k = 1, n< h~K 



otherwise. 



Thus we proved 
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Proposition 9.2.3. For /j, < h 1 under condition (0.7) estimate (9.16) 
holds and under condition (0.8) estimate (9.18) holds. 



9.2.3 Magnetic Weyl approximation 

Finally, let us deal with magnetic Weyl approximation. Now we should add 
to what we got for K T corresponding expressions with R(a, , ~f) = R MW (a, , ~f) 
calculated according to (1.83) for /i < h^ 1 and (1.85) for /i > hr 1 . 

(i) Consider case fi < h^ 1 first. Then according to (1.83) 7 snaps to 
min(/i _1 , fiha^ 1 ) and R(a, 7) = jih^ 1 and (9.3) becomes 



(9.19) |k mw i < cir l h' 1 + cy/r 




h~* < /i < h 



-1 



Meanwhile in the first term of (9.5) a snaps to its largest possible value 
i.e 1. Further, in the second term a snaps to its largest value which is /3 and 
then (5 snaps to its largest value (yu/?)2 as h~3 < fi < h^ 1 and it becomes 
Cfi 2 h- 2 {fih) 2 ~^. 

Furthermore, contributions of zones (9.9) and (9.10) are estimated as 
before; we arrive in the end to 




2-K 



(9.20) |K MW | < C>- 2+K /T 2 + Cfi 2 h~ 2 { " ' 11 ~ H 2 ' + 



C 



n 2 h-\nh) 2 -* K h'l <n< h-\ 



otherwise, 
(ii) As /ih > 1 7 snaps to /i^/jf and (9.3) becomes 
(9.21) |K MW | < Cn 1+1 2 K h~ 1+1 2 K 

Consider the first term in expression (3.5); integrand contains a in the 
positive power and therefore a snaps to 1. Meanwhile in the second term 
both a > Co/i~2h2 and (3 > C fi~2h2 are in the negative powers and snap 
to yU~2/j2. Estimating as before contributions of zones (9.13) and (9.14) we 
arrive to estimate (9.21) again. 

All these above estimates are worse than for K T ; thus we arrive to 
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Proposition 9.2.4. (i) For u, < h 1 under condition (0.7) estimate (9.19) 
holds and under condition (0.8) estimate (9.20) holds; 

(ii) For a, > h~ l under either condition (0.7)', (0.8)' estimate (9.21) holds. 



9.3 Case d = 3 

We consider case d = 3. Note that 

Remark 9.3.1. (i) We use notations in (9.5), (9.6): A instead of a and A' 
instead of (3 (in the second term only); 

(ii) We cannot expect estimate better than 0(h~ 4 ) for /ih < 1. Meanwhile 
under condition (0.8) contribution of zone 

(9.22) {(x,y): A(x) < A, A(y) < A} 

with A = Cq^T 1 to any estimate does not exceed Cfi 2 h~ 4 x /i _6+K = 0(h~ 4 ). 
We can even take A = /i~ 2 /( 6 ~ K ) here. 



9.3.1 Tauberian estimates 

We will use results of subsection 5.4 combined with arguments of subsubsec- 
tion 6.8.1 always setting 73 = 7. 



Case fih < 1 As 73 = 7 we have 2/(7) x A7, £ x A and (5.25) becomes 

yu/T^As as A 2 7 < h, 

(9.23) /? T ( 7 ) = <J / //7- 1 A-i 7 -i as A 2 7 > /?, A7 < ///?, 

h~^(fih/\'y) 1 as A7 > /i/?. 

To estimate |K T | we need to calculate 

(9.24) K(A):= |(/? T (7))V^ ^7 
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and as A > /j, 1 it is equal to Cfx 2 h 2 A 1 7 2 K calculated for 7 = min(/i/j/A, A) 
as < k < 2 and for 7 = h/X 2 as 2 < k < 3 resulting in 

(9.25) K(A) < C/?- 4 A 3 - K + 

r Cfi 2 h- 2 X^{h/\ 2 ) 2 - K 2<k<3, 

Cl C/x 2 /j- 2 A- 1 (1 + j logAyu)) k = 2, 

[/i^-^-^in^VA.A) 2 "" 0<k<2. 

In particular, we conclude that 

(9.26) As fih < 1, \W/F\ x 1 

/i 2 /?~ K 2 < k < 3, 

(9.27) |K T | < C/)- 4 + C<J fi 2 h- 2 \\ogfi\ k = 2, 

< k < 2. 



On the other hand, under assumption (0.8) we need to calculate the first 
term in the right-hand expression of (9.5) which is equal to j K(X)X 2 d\ 
and one can see easily that it is (9.25) calculated as A = 1. 

Finally, one needs to calculate the second term in the right-hand expres- 
sion of (9.5) and we leave to the reader to prove that it allows the same 
estimate. Thus we arrive to statement (i) of proposition 9.3.2 below. 

On the other hand, assume that assumption (6.57) is fulfilled i.e. | Vj_f V/F\ 
1. Then (5.25) becomes 

(9.28) R T (j) = 

' fih- 1 (1 + |log(/x/j)|) as 7 < h, 

/J/7 -1 (l + (log(l///7))+) as h < 7 < m'mlfih.fi' 1 ), 

fih^ 1 ^? as fi^ 1 < 7 < fih, 

ah h~^(/j,h/^))' as 7 > fih. 

and in this case instead of 7 = h the lower pick is 7 = /i" 1 in Cu 2 br 2 ^ 2 ~ R 
which becomes Ch~ 2 ji K and we arrive to statement (ii) of proposition 9.3.2: 
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Proposition 9.3.2. Let fih < 1. Then 

(i) Under either non- degeneracy condition either (0.7) or (0.8) estimate 
(9.27) holds. 

In particular, |K T | = 0(h~ 4 ) as either < k, < 2 or k = 2 and 
fi < (h\ log/?|) _1 or 2 < k < 3 and fi < h~ 2+ 2 K (and this is true for sure as 
H < h~\). 

(ii) Under non- degeneracy condition (6.57), 2 < k < 3 

(9.29) |K T | < Ch- 4 + Ch- 2 fi K . 

In particular, |K T | = 0(h~ 4 ) as ji < hr 2 l K (and this is true for sure as 
/i < h-l). 

Case fih > 1 

Proposition 9.3.3. Let u,h > 1. Then under either non- degeneracy condi- 
tion (0.7)' or (0.8)' estimate 

{/i 2 /T K 2 < k < 3, 

/i 2 /r 2 |lo g/ u| k = 2, 
< k < 2. 

holds. 

Proof. An easy proof based on proposition 5.5.3 is left to the reader. □ 

9.3.2 Weyl asymptotics 

Proposition 9.3.4. Let fih < 1. Then 

(i) Under non- degeneracy condition (0.7) 

(9.31) |K W | < Or 4 + C> 3 /r 3 (^) 3 ~ K ; 

(ii) Under non- degeneracy condition (0.8) 

(9.32) |K W | < Ch- 4 + Cu. 3 h- 3 {nh)^ 3 - K) ; 
(Hi) As /i < under non- degeneracy condition (6.57) 

(9.33) |K W | < Ch- 4 + C/2 5 h- 2 (i2h) 3 - K . 

Proof. An easy proof based on proposition 5.4.3 is left to the reader. □ 
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9.3.3 Magnetic Weyl asymptotics 

We will apply proposition 5.6.4. In the isotropic settings R MW (j) given by 
(5.47) becomes 



(9.34) /? MW ( 7 ): 



A~^7~i (fih/X'y) 1 

A^/T^ 



as 7 > fxh/X, 

1 2 

as A~5/)3 < 7 < /i/?/A, 

12 

as 7 < A 3/ 7 3 j 



which immediately implies statements (i), (ii) of proposition 9.3.5 below. 
Meanwhile as a x 1, fih < 1 (5.47) becomes 



(9.35) /?r (7) := 



7 2(^/7/7)' 



C/i < 



7 h^^H^/i) 1 

k 7-5(/i 7 )^ 

which immediately imply statement (iii) of proposition 9.3.5 below. 



as 7 > yu/?, 1 < /x 2 /?, 

as 7 < jih, /i7 > 1, 

as 7 > /i/?, 1 > /x 2 /j, 

as 7 < jih, /i7 < 1 



Proposition 9.3.5. Let jih < 1. T/ien under either non- degeneracy 
condition (0.7) or (0.8) 



(9.36) 



IK 



MW 



< C/7- 4 + C// 2 /7-K 2+K ); 



(nj Lei jih > 1. Then under either non- degeneracy condition^). 7)' or (0.8)' 



(9.37) 



IK 



MW 



2 

(raj Let ji < h~3. Then under non-degeneracy condition (6.57) 



(9.38) 



IK 



MWi 



< Ch- A +Cfi 3 h-^ 
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9.4 Problems 

Problem 9.4.I. For d = 2, 3 estimate |J - J T |, |J - J w \ and |J - J MW \. 

Problem 9.4.2. For d = 2, 3 consider pilot-model approximation. 

Problem 9.4-3. For d = 2,3 compare estimates for |K W | and |J — J7 W | to 
those for |K MW | and |J — JT MW |. 

Problem 9.4.4. For d = 2,3 as fi < h^ 1 get rid off condition (0.5) replac- 
ing it by (0.7), or even (0.8) + ; or even (0.8). 

Problem 9.4.5. Get rid off condition (0.6) assuming instead that \F\ + 
|VF| > e for d = 2 (which will affect estimates for sure) and |F| + |V<S>F| > e 
for d = 3 (which most likely would not affect estimates) 

Problem 9.4.6. For d = 3 investigate K T , K w , K MW under non-degeneracy 
conditions (0.7) and |V ±fV / 'F\ < e ==> | det Hess_|_ F V/F\ x 1. 
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